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Abstract 

A  systematic  study  of  microwave  discharges  at  2.45  GHz  has  been  performed 
through  the  pressure  range  of  sustainable  electric  discharges  in  pure  oxygen  flows  of  2  to 
10  Torr.  A  corresponding  study  of  13.56  MHz  has  also  been  perfonned  at  pressures  of  2, 

4,  and  7  Torr.  Optical  emissions  from  02  (a  'A),  02(b  1 X) ,  and  O-atoms  have  been 

measured  from  the  center  of  a  p-wavc  discharge.  Discharge  residence  times  from  0.1  to  5 
ms  have  been  reported,  establishing  that  gas  temperatures  arrive  at  stationary  values  within 
the  first  0.3  ms  upon  entering  the  discharge  region.  The  0-,{b  X)  emissions,  with  a 

spectral  resolution  of  0.01  mn,  have  been  used  to  measure  the  temperature  of  the  gas, 
which  typically  reaches  a  steady-state  of  1 ,200  K.  A  theoretical  description  of  the  gas 
heating  is  fit  to  measured  temperatures,  which  detennines  that  the  fraction  of  discharge 
energy  coupled  into  gas  heating  is  17  ±  2%.  The  yield  of  02(a  'A)  comes  to  steady-state 
at  all  pressures  within  1  ms  of  entering  the  discharge  region.  The  interpretation  of  the 
measured  02  ( a  1  A)  yield,  using  a  streamlined,  nearly  analytic  model,  cast  new  light  on  the 
kinetics  within  the  electric  discharge.  The  pseudo-first  order  quenching  rate  of 
02{a  1  A)  ranges  from  6,000  s'1  for  p-wavc  discharges  to  600  s'1  for  radio  frequency  (RF) 

discharges,  independent  of  gas  pressure  and  flow  rate.  The  slower  decay  rate  for  the  RF 
discharge  corresponds  with  a  considerably  lower  ionization  rate.  The  observations  are 
consistent  with  a  second  order  reaction  channel  that  is  dependent  on  both  the  electron  and 
molecular  oxygen  ground  state  concentrations.  Destruction  of  the  02  (a  1  A)  state  by  direct 


IV 


impact  with  electrons  or  atomic  oxygen  does  not  adequately  describe  the  observed  behavior 
of  02(a  1  A)  .  The  role  of  vibrationally  excited  ground  state  oxygen  is  explored  and 
provides  a  plausible  destruction  mechanism. 
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KINETICS  OF  THE  ELECTRICAL  DISCHARGE  PUMPED 


OXYGEN-IODINE  LASER 

I  Introduction 

The  airborne  laser  (ABL)  is  a  modified  Boeing  747  designed  to  serve  as  a  missile 
defense  system.  The  weapon  system  carried  on  the  ABL  is  a  chemical  oxygen-iodine 
laser  (COIL).1  COIL  is  a  megawatt-class  laser  that  employs  the  interaction  of  basic 
hydrogen  peroxide  and  chlorine  gas  to  reactively  generate  the  first  excited  state  of 
molecular  oxygen,  02(a  ’A).1'2  The  02(a  1  A)  state  is  then  used  to  collisionally  excite 
the  /( 2Py2)  excited  state  of  atomic  iodine.  It  is  the  population  inversion  between  the 
/( 2Py2 )  upper  energy  state  and  I(2P3/2)  ground  state  that  provides  the  gain  in  a  COIL 

device.  Using  a  chemical  reaction  to  generate  the  02(a  A)  state  of  oxygen  introduces  a 
limit  to  the  magazine  depth  of  the  weapon  system  and  possesses  a  large  logistical  tail  and 
significant  safety  issues  associated  with  the  exotic  chemistry.  In  order  to  alleviate  the 
shortcomings  of  the  COIL  system,  an  electric  discharge  may  be  used  to  populate  the 
02{a  'A)  state. 

An  electrically  driven  oxygen-iodine  laser  (EOIL)  has  many  advantages.  First,  it 
eliminates  the  logistic  tail  associated  with  the  COIL  chemistry.  Second,  it  has  an 
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extended  magazine  depth.  Next,  it  ffers  closed-cycle  operation.  Finally,  it  maintains  the 
thermal  management  and  beam  quality  traits  of  the  COIL. 

The  yield  of  02  ( a  1  A)  produced  by  a  microwave  discharge  was  first  reported  in 
1978  by  Bernard  and  Pchelkin.-  They  obtained  yields  (the  fraction  of  gas  flow  in  the 
excited  state)  of  1 1%.2  The  02{a  'A)  yield,  Ya ,  is  defined  as: 


where  [02(a  *A)J  is  the  concentration  of  02{a  'A)  in  the  gas  flow,  and  [O2]0  is  the 

concentration  of  ground  state  molecular  oxygen  entering  the  discharge  region.  In  a  laser 
medium,  there  is  competition  between  absorption  and  stimulated  emission  of  photons. 
The  excited  state  must  reach  a  yield  (percent  of  population)  large  enough  to  produce 
more  stimulated  emissions  of  photons  than  those  absorbed.  This  yield  is  defined  as  the 
threshold  yield.  The  threshold  yield  required  for  lasing  in  an  oxygen-iodine  laser  is  17% 
at  room  temperature,  or  293  K.  Therefore,  using  an  electric  discharge  as  a  source  of 
02(a  'A)  for  an  oxygen-iodine  laser  will  require  the  yield  of  the  02(a  ‘A)  within  the 
discharge  to  increase  in  order  to  surpass  the  threshold  yield  of  the  laser. 

Recent  research  into  an  EOIL  device  has  resulted  in  output  laser  powers  up  to 
102  W.4  In  2004,  the  University  of  Illinois  at  Urbana-Champaign  (UIUC)  demonstrated 
the  first  operation  of  an  EOIL  device.5  In  2005,  the  maximum  02(a  'A)  yield  of  17% 

•3 

was  reported  without  the  inclusion  of  NO2  in  the  gas  flow.  The  maximum  laser  power 
of  this  early  laser  system  was  520  mW  with  a  400  W  radio  frequency  (RF)  discharge.  In 
2008,  two  parallel  1,000  W  discharges  were  used,  and  a  laser  power  of  12  W  was 
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measured.4  In  2010,  UIUC  presented  output  powers  of  102  W  for  an  EOIL  device,  using 
a  discharge  configuration  with  six  3/4-inch  quartz  tubes.6  For  an  EOIL  device  to  be 
useful  as  a  laser  weapon,  the  output  laser  power  must  be  scaled  to  kilowatts. 

The  maximum  laser  power  of  an  EOIL  device  is  constrained  by  the  amount  of 
02{a  1  A)  within  the  gas  flow  over  the  threshold  yield.  In  order  to  scale  an  EOIL  device 
to  higher  powers,  either  the  yield  must  be  increased,  or  the  number  density  of  the  gas 
above  threshold  must  be  increased.  The  number  density  above  threshold  yield  may  be 
improved  by  increasing  the  total  yield  or  the  total  gas  pressure,  while  maintaining  the 
same  yield.  The  overall  yield  of  the  discharge  is  detennined  by  the  nature  of  the  plasma 
and,  therefore,  is  difficult  to  change.  Increasing  the  pressure  of  the  oxygen  flowing 
through  an  electric  discharge,  while  maintaining  02  ( a  1  A)  yields  above  the  laser 
threshold  of  17%,  has  proved  problematic.  The  pressure  dependence  of  the  yield 
suggests  that  a  second-order  reaction  is  destroying  the  02 (a  'A)  excited  state.  Previous 
studies  in  the  literature  have  focused  primarily  on  increasing  production  of  02  ( a  1  A)  and 
not  on  completely  understanding  the  destruction  of  02(a  'A)  in  the  discharge.3'14 

In  order  to  gain  a  better  understanding  of  what  is  limiting  the  02  ( a  1  A)  yield  in 

the  electric  discharge,  a  systematic  study  is  performed  of  a  glow  discharge  sustained 
within  flowing  pure  oxygen  gas.  The  chemical  kinetics  of  an  EOIL  system  are  studied  in 
order  to  investigate  the  dominant  mechanism  responsible  for  limiting  the  yield  of 
02(a  1  A)  as  pressure  is  increased.  Studying  chemical  kinetics  within  an  electric 

discharge  is  complicated  because  of  the  need  to  employ  diagnostics  that  do  not  perturb 
the  plasma.  Typical  studies  used  the  optical  emissions  from  excited  states  within  the 
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discharge  to  monitor  the  excited  species  populations.  ’  ’  ’  Using  optical  diagnostics  to 
examine  the  system  is  complicated  by  the  variations  in  gas  density  due  to  changes  in  the 
gas  temperature,  as  well  as  contamination  of  the  collected  data  from  emissions  of  atomic 
species  within  the  measured  frequency  range. 

In  this  research,  there  has  been  an  in-depth  investigation  of  the  changes  in  gas 
temperature  as  gas  flow  rates  and  gas  pressures  change,  using  the  emissions  of  the 
02(b  X)  state.  The  02(b  'X)  state  is  the  second  excited  state  of  molecular  oxygen,  and 
the  relative  intensities  of  the  peaks  in  the  rotational  emission  spectrum  are  governed  by 
Boltzmann  statistics,  allowing  for  the  gas  temperature  to  be  extracted.  Gas  temperatures 
are  reported  in  a  number  of  studies  with  errors  typically  being  10%  or  higher.  ’  ’  ’  ’ 
Emission  intensities  are  directly  related  to  the  number  density  of  emitters,  which  is  then 
directly  related  to  the  gas  temperature  by  the  ideal  gas  law.  Uncertainty  in  the 
temperature  translates  directly  into  uncertainties  in  the  determination  of  the  species’ 
concentrations.  Using  a  high-resolution  spectrometer  and  developing  a  theoretical 
description  of  the  heat  transfer  within  the  gas  flow,  a  high-precision  gas  temperature 
measurement  was  pursued. 

Gas  temperature  also  plays  an  important  role  in  laser  weapon  performance.  The 
intent  of  coupling  energy  into  the  gas  using  the  discharge  is  to  populate  the  02  ( a  1  A) 

state.  If  this  energy  is  instead  coupled  into  gas  heating,  it  is  wasted  and  must  be  removed 
from  the  gas  flow.  Also,  heavy-heavy  particle  interactions  are  typically  directly 
proportional  to  the  temperature  of  the  gas.  One  such  interaction  is  the  quenching  of  the 
02  ( a  1  A)  state  by  ground  state  oxygen.  The  relationship  between  the  excited  state  of 
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oxygen  and  excited  state  of  iodine,  which  is  responsible  for  lasing  in  oxygen-iodine  laser 
devices,  is  inversely  affected  by  temperature.  This  will  be  discussed  further  in  Section 
2.1.  It  is  important  to  have  a  precise  and  complete  understanding  of  the  gas  temperatures 
within  the  electric  discharge  in  order  to  understand  the  performance  of  an  EOIL  device. 

After  temperature-driven  changes  in  concentration  are  removed  from  the  data,  the 
measured  02  ( a  1  A)  emissions  may  still  be  affected  by  optical  emissions  from  other 
excited  states  created  by  the  discharge.  Within  the  electric  discharge,  other  excited 
species  are  formed,  such  as  O-atoms,  and  these  may  emit  optically  within  the  same 
spectral  region  as  the  molecular  species,  02(a  'A).  When  observing  the  02(a  ‘A) 
emissions,  the  spectral  radiation  from  four  atomic  species  is  also  gathered.  The  atomic 
species  emit  at  1,257  nm,  1,266  nm,  1,280  nm,  and  1,299  mn  and  are  in  the  same  spectral 
region  as  the  02(a  ’A)  excited  state,  which  radiates  from  1,250  nm  to  1,290  nm.  By 
identifying  and  removing  the  optical  contributions  of  the  O-atom  states  from  the 
measured  intensities,  the  02  ( a  'A)  concentrations  may  be  studied  directly  within  the 
discharge.  After  taking  care  to  report  concentrations  of  emitting  species  accurately,  a 
kinetic  model  is  included  in  the  analysis  to  understand  the  behavior  of  non-emitting 
species,  such  as  electrons. 

In  order  to  understand  the  kinetic  reactions  within  the  discharge,  a  computer 
model  is  required.  The  electron  density  drives  the  reactions  of  all  other  species  within  the 
discharge;  however,  electrons  do  not  have  optical  signatures  [such  as  02(a  'A)  ]  that 
would  enable  a  measurement  of  their  concentration.  Therefore,  the  electron  number 
density  is  determined  by  numerically  solving  the  Boltzmann  equation  and  using  the 
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results  to  solve  a  chemical  kinetic  package  to  describe  the  plasma.  There  are  many 
computational  models  accepted  within  the  community.  The  models  include  at  least  50 
reactions  and  reproduce  measured  results.7, 16  However,  including  such  a  large  number  of 
kinetic  pathways  makes  the  determination  of  key  reactions  difficult.  Therefore,  a 
physically  intuitive  model  is  developed  within  this  study  that  adequately  describes  the 
experimental  results  and  gives  a  clear  understanding  of  the  physics  involved  in  the 
system. 

A  physically  intuitive  set  of  kinetics  is  included  in  a  computer  model  and  used  to 
extract  the  behavior  of  the  quenching  reaction  responsible  for  limiting  the  02  ( a  1  A) 
concentration  within  electric  discharges.  The  model  includes  the  electron  density,  O- 
atom  density,  negative  ion  density,  ground  state  oxygen  density,  and  density  of  the 
02(a  'A)  excited  state  of  molecular  oxygen.  This  model  is  shown  to  accurately  describe 
the  discharge  and  gives  a  clear  understanding  of  the  gas  phase  kinetics,  which  is  difficult 
to  discern  from  models  that  include  larger  numbers  of  kinetic  pathways.  The  results  of 
the  physically  intuitive  computer  model  are  used  to  describe  the  behavior  of  the 
quenching  mechanism  responsible  for  limiting  the  yield  of  singlet  oxygen  within  RF  and 
p-wave  electric  discharges  (13.56  MHz  and  2.45  GHz,  respectively). 

The  destruction  rate  of  02(a  'A)  increases  by  an  order  of  magnitude  between  the 
RF  and  p-wavc  electric  field  frequencies.  This  drastic  change  links  the  destruction 
mechanism  to  electron  concentrations.  The  02(a  ’A)  destruction  rate  extracted  from  the 
data  is  invariant  as  pressure  and  velocity  of  the  gas  flow  through  the  discharge  changes. 
The  role  played  by  electrons  as  well  as  atomic  oxygen  in  the  destruction  of  02  ( a  1  A)  will 
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be  explored.  Furthermore  destruction  of  02 ( a  'A)  by  collisions  with  vibrationally 
excited  oxygen  may  play  an  important  role  in  the  plasma  kinetics. 
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II  Background  and  Theory 


This  chapter  will  present  the  theory  and  background  relevant  to  the  study  of 
oxygen  plasmas,  starting  with  a  general  discussion  of  the  operation  of  COIL  devices  and 
then  the  operation  of  EOIL  devices  with  an  emphasis  on  the  advantages  and  technical 
challenges  of  an  EOIL  device.  The  optical  diagnostics  used  to  measure  the  change  in 
relative  magnitude  of  the  excited  species  without  perturbing  the  plasma  are  also 
discussed.  In  the  case  of  02(a)  and  02(b) ,  the  optical  signatures  involve  molecular 

spectroscopy;  its  theory  is  well  established  and  discussed  in  this  chapter.  The 
actinometry  technique  is  used  to  determine  the  O-atom  concentration  within  the 
discharge  region.  This  technique  and  its  associated  spectroscopy  are  also  discussed.  The 
temperature  of  the  gas  changes  as  it  flows  through  the  plasma  region,  so  a  theoretical 
description  of  gas  heating  and  cooling  is  discussed.  To  improve  the  understanding  of  the 
chemical  kinetics  occurring  within  the  plasma  from  changes  in  the  concentrations  of 
measured  excited  species,  a  model  of  the  chemical  kinetics  has  been  developed  and  is 
described  in  Section  4.6.  The  rate  package  used  in  this  model  is  dependent  on  the 
electron  energy  distribution  within  the  plasma.  A  preexisting  Boltzmann  solver,  called 
BolSig+,  is  used  to  determine  the  rate  coefficients  for  the  model  used  in  this  study. 
Finally,  the  current  status  of  research  in  the  area  of  EOIL  development  is  discussed 
within  the  context  of  power  scaling  and  increasing  the  yield  of  02(a )  in  oxygen  plasmas. 
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2.1  Chemical  Oxygen-Iodine  Laser 

The  COIL  system  was  first  demonstrated  in  1977. 19  The  electrically  excited  state 
of  oxygen,  02(a )  ,  is  populated  by  a  chemical  reaction  between  chlorine  gas  and  liquid 
basic  hydrogen  peroxide.  This  is  a  multiphase  process  governed  by  the  stoichiometry: 

Cl2  +  H202  +  2 KOH  kco,^->  02  ( a  1  A)  +  2  H20  +  KCl  (2. 1) 

20 

The  radiative  lifetime  of  the  02(a)  state  is  64  minutes.  This  excited  state  is  also 

resistant  to  collisional  deactivation  at  surfaces.  Therefore,  once  it  is  created  by  the 
chemical  reaction,  this  state  will  persist  for  long  periods  of  time.  A  population  inversion 
between  the  02(a)  and  02{X )  is  easily  achieved;  however,  the  cross-section  for  the 
stimulated  emission  is  small.  The  COIL  system  uses  the  near  resonant  energy  transfer 
from  the  electrically  excited  02(a)  state  and  atomic  iodine  I(2Py2)  state  to  create  a 
population  inversion  in  the  iodine. 

The  /(  1 P21)  and  02(a)  excited  states  are  separated  by  only  0.03  eV,  allowing  a 
very  efficient  transfer  of  the  energy  from  the  01  ( a )  state  of  oxygen  to  the  iodine  lasing 
state.  However,  the  iodine  injected  into  the  gas  flow  is  not  atomic;  it  is  in  the  molecular 
state.  The  molecular  iodine  must  be  dissociated  first  before  the  atomic  iodine  can  be 
excited.  The  energetics  for  the  dissociation  of  molecular  iodine  are  also  shown  in  Figure 
2.1  and  require  the  sacrifice  of  at  least  two,  but  typically  five,  02(a)  molecules  in  gas 

flows  of  pure  02(a) .  For  the  EOIL  system,  L  is  dissociated  via  many  pathways, 
including  collisional  dissociation  with  the  02(b )  state  and  O-atoms.  ’  ’  ’ 
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Figure  2.1  Energy  level  diagram  for  COIL,  illustrating  the  resonant  energy  transfer  from  metastable  reservoir 
02(a  1  A)  to  the  upper  lasing  level  of  atomic  iodine  /( 2  P]n )  •  The  energy  levels  for  several  key  states  involved  in 
the  dissociation  of  molecular  iodine  are  also  provided.24 


A  schematic  for  a  simple  COIL  device  is  shown  in  Figure  2.2.  As  is  observed  in 
the  schematic,  the  waste  water  and  salt  from  the  chemical  reaction  is  removed  from  the 
gas  flow.  The  gas  is  cooled  by  supersonic  expansion,  while  iodine  is  injected  into  the  gas 
flow.  The  chemicals  are  then  exhausted  from  the  laser  cavity,  taking  any  waste  heat  with 
them. 

The  gain  in  the  system  is  generated  through  the  population  inversion  of  /( 2/j  2 ) 
versus  the  I(2Py2)  state  of  the  iodine  atom,  which  emits  at  k  =  1.315  pm.  The 
population  inversion  is  attained  by  energy  transfer: 

02  (a 1  A)  +  /( : 2PV2 )  < - >  02  (X 3 2)  +  /* ( : 2PU2 )  (2.2) 
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Figure  2.2  Schematic  view  of  chemical  oxygen-iodine  system,  including  supersonic  expansion  of  gas  (adapted 
from  Carrol).3 


Assuming  that  this  is  the  only  important  kinetic  reaction  occurring,  the  steady- 
state  relationship  between  the  excited  states  and  their  ground  states  is: 


where 


(2.3) 


oJ(jr!z)^=[o!(x'!z)]  +  [o!(v)]  +  [o4a1A)]  +  [oJ((.12;)]  +  ^[o(Jp)]  +  i.5[oJ] 


20 

And,  from  Perram  et  al.  : 

Keq(T)*  0.75 exp 


'  4°2  Z  ' 

V  /  1gJ 


(2.4) 


(2.5) 


where  Tg  is  given  in  Kelvin. 

The  ratio  of  I  ( 2Pl/2 )/ 1  ( 2P^2 )  is  directly  related  to  the  yield  of  the  02(a)  that 
results  from  the  singlet  oxygen  generator  (SOG).  Note  that  the  equilibrium  constant  is 
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temperature-dependent.  By  cooling  the  gas,  the  yield,  Ya ,  required  to  reach  the  threshold 
is  decreased.  Note  that  the  value  of  [O2]0  (used  in  detennining  Ya)  in  a  conventional 
COIL  is  [02(X)]  +  [02(a)] ,  because  the  other  species  do  not  obtain  significant 
concentrations.  The  threshold  for  gain  (at  Tg  =  295  K)  is  achieved  when  Ya  is 

approximately  15%.  By  virtue  of  Keq,  the  threshold  yield  is  reduced  to  approximately  5% 
at  the  point  of  supersonic  expansion  by  reducing  the  gas  temperature  to  160  K.  Typical 
yields  from  the  RotoCOIL  SOG  used  in  a  COIL  are  80%.26 

2.2  Electric  Discharge  Pumped  Oxygen-Iodine  Laser 

The  primary  difference  between  an  EOIL  and  a  conventional  COIL  system  is  the 
SOG.  Changing  this  component  introduces  limits  to  yield,  temperature  issues,  and 
excited  species  that  do  not  exist  in  the  conventional  chemical  system  and  will  require 
adjustments  to  the  laser  system.  In  Figure  2.3,  a  schematic  of  a  typical  EOIL  system  is 
shown. 

In  Figure  2.3,  the  system  exhaust,  laser  cavity,  and  supersonic  gas  expansions  are 
all  very  similar  to  those  shown  for  the  conventional  COIL  system.  The  system’s  SOG  is 
an  electric  discharge. 

The  major  differences  between  a  chemistry  based  SOG  and  a  discharge-based 
SOG  are:  (1)  Ya  is  smaller,  approximately  10-15%;  (2)  gas  temperatures  within  the 

plasma  are  much  higher,  commonly  reaching  temperatures  of  700  K  and  higher;  and  (3)  a 
myriad  of  other  excited  species,  such  as  O-atoms  and  vibrationally  excited  Ot,  are 
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formed,  possibly  limiting  the  maximum  achievable  02(a)  population  and  reacting  with 

2  3 

I*,  which  will  reduce  the  achievable  gain  within  the  discharge.  ’ 


Laser 

Cavity 


Figure  2.3  Schematic  of  a  typical  EOIL  system,  including  a  discharge-driven  singlet  oxygen  generator,  a 
supersonic  flow  to  cool  the  gas  before  the  laser  cavity,  iodine  injection,  and  the  laser  cavity  itself  (adapted  from 
Carrol  et  al.).3 


2.3  Calculating  Electron  Energy  Dependent  Rates 

Gas  discharge  models  require  the  use  of  rate  coefficients,  which  are  dependent  on 
the  electron  energy  distribution  within  the  plasma,  as  shown  below: 

ki(E  /  N)  is  the  rate  coefficient  of  the  ith  reaction;  e  is  the  fundamental  charge  of  an 
electron;  m  is  the  electron  mass;  f(s;E/  N )  is  the  electron  energy  distribution  function 


f  9  \/2  00 

,(E/N)  =  —\io-,(e)e'l2f(e;E/N)de 

\m!  i 
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(EEDF);  a  is  the  electron  energy;  cr/A)  is  the  ith  electron  impact  cross-section;  E  is  the 
electric  field  strength;  and  N  is  the  number  density  of  collision  partners.  In  this  work,  N 
will  be  the  ground  state  of  molecular  oxygen,  because,  under  most  conditions,  70%  or 
more  of  the  gas  flow  remains  in  the  ground  state.  In  plasmas,  electrons  are  accelerated  by 
the  electric  field  and  slowed  by  collisions,  so  a  key  parameter  of  a  plasma  is  the  value  of 

97 

E/N  or  the  reduced  electric  field." 

The  BolSig+  solver  provides  a  steady-state  solution  of  the  Boltzmann  equation  for 
electrons  in  a  uniform  electric  field,  using  the  classic  two-term  expansion."  At  high 
values  of  E/N  (where  most  collisions  are  inelastic)  and/ (the  distribution  function 
becomes  anisotropic),  the  two-tenn  approximation  will  fail.  The  errors  in  calculated 
transport  coefficients  are  acceptable  for  fluid  discharge  modeling  in  the  range  of  most 
studies.  The  typical  reaction  rates  for  this  study  are  shown  in  Table  2.1.  In  this  table, 
the  reaction  rate  for  elastic  collisions  is  5.9  x  1CT8  cm3  s'1,  which  is  much  larger  than  the 
reaction  rates  for  inelastic  interactions  of  electrons  with  ground  state  oxygen.  Therefore 
the  conditions  of  this  study  are  well  within  the  conditions  where  the  two-tenn  expansion 
is  valid. 

In  general,  the  EEDF,  f£(a)  =  2/ \[k  (kBTe)  3  ~  a''2  exp  \-a/kBTe\,  maybe 

calculated  from  the  fundamental  collisional  cross-section  data  (of  the  gases  in  which  the 
plasma  is  being  sustained)  by  solving  the  Boltzmann  equation.  The  average  energy  of  the 
electrons,  s  ,  is  the  average  value  of  a  using  the  EEDF  : 

oo 

£avg  =  J  £fe(£)ds  (2-7) 

0 
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Table  2.1  Reaction  rates  for  this  work.  Values  calculated  at  4  Torr,  Tg  =  1,185  K,  and  Te  =  3.2  eV 
_  (given  in  units  of  cm3  s'1  with  the  exception  of  reaction  24  in  units  of  cm6  s'1). _ 


Reaction 
Number,  k 

Name 

Reaction 

Reaction  Rate 

Coefficient 

Source 

i 

Elastic 

Collisions 

e  +  [02(X)]  >[02(X)]  +  e 

kElas=  5.9  xio^8 

28 

2 

Dissociative 

Attachment 

e-  +  o2(X)  - . *->o~  +  o 

^attach (Te)  =1.9x10  ” 

28 

3-9 

Vibrational 

Excitation 

e+02  ( X )  — - — >  02  ( v )  +  e~ 

v  =  0.02  to  0.75eV 

=  5.3x10-" 

28 

10 

Excitation 

0.98  eV 

e~  +  02  ( X )  — ka+  >  02  ( a )  +  e~ 

A-a+(Te)  =4.2x1 0“'° 

28 

11 

Excitation 

1.63  eV 

e-  +  02  (X)  kb+-  >  <9,  ( b )  +  X 

A-b+(Te)  =9.9x10-" 

28 

12 

Excitation 
4.50  eV 

e+O^X)  kA+-  >02(A,3T)  +  e~ 

A-a+(Tc)  =  1.6  x  1010 

28 

13 

Excitation 

6.0  eV 

e ~  +  <9,  (X)  -°±-  >  <9( 3P)  +  0(3P)  +  e- 

*op+(Te)  =2.2xlO-10 

28 

14 

Excitation 

8.4  eV 

e~  +  02  (X)  ko+  >  0(3P)  +  0('D)  + 

£0D+(Te))  =  2.3  x  10  10 

28 

15 

Excitation 

9.97  eV 

e~  +  02 (X)  -°±-  > 0(3P)  +  0(lS )  +  e 

fcos+(Te))  =1.7xl0-13 

28 

16 

Ionization 
12.06  eV 

e~  +02(X ) — kj£g— >02  +  e~  +  e~ 

k ion  (T e)  )  =1.0xl0-12 

28 

17 

Deexcitation 
0.98  eV 

e  +  07  (a)  — — — >  02  ( X )  +  e 

i’a-(Te)  =5.4x10-'° 

Detailed 

Balance 

18 

Recombination 

e"  +  O;  (X)  *"  >  0(  3P)  +  0(  3P) 

Arecde)  =1.6  X  10  8 

20 

19 

Atomic 

Recombination 

0  +  0  +  Wall  K“u  >6>,(X) 

^wall  =  2400 

20 

20 

Detachment 

0~  +02(a ) — ^-><93  +  e“ 

kode=  3.0xl0-10(rg/300)°5 
=  6.0x10" 

20 

21 

Ionic 

Recombination 

O  +  O ,+  (X)  **■  >  O  +  O,  (X) 

k0rec  =  2-0  X  10-7  [Tg  /300)  ‘ 

=  5.1xl0-8 

20 

22 

Quenching  by 
02(X) 

[<92(a)l  +  [02(X)]  — »  2[0,(X)] 

_20Q/ 

k02X  =  3.0  x  lCT18e  ^ 

2.5  xlO"8 

29 

23 

Quenching  by 
O-atoms 

[0]  +  [02(a)]-+0  +  [02(X)] 

kaQ  =  2.0  x  10"16 

30 

24 

Three  Body 
Reaction 

[O]  +  [0,(X)]  +  [02(a)]  -+  [O]  +  2[0,(X)] 

^=l-0xl0-32 

17 

25 

Quenching  by 
Ozone 

[03]  +  [02(a)]^[03]  +  [02(X)] 

-2840/ 

kam  =2.0xl0-"e  /Tg 
=  1.8xl0-12 

29 

26 

Electron 

Interaction 

[e“]  +  [O,  (a)] -+[e"  ]  +  +!// 

k  *  f(s  )  x  10-10 

=  6.7x10-'° 

Scaled 

27 

Vibrationally 
Excited  02 

[02(v)]  +  [02(a)]  -+  [02m  +  [02(X)] 

Jfcv  «  4x10-'° 

This 

Work 

28 

Pooling 

[02(a)]  +  [02(«)]->[02(X)]  +  [02(fc)] 

kpooi  =  9  xl0-17e™' 

=  8.2  x  10  18 

30 

29 

Quenching  of 
O-  by  O-atoms 

0~  +  M  — >0  +  M  +  e~ 

fro*®  =8x10-'° 

31 

30 

0(‘D) 

quenched 

02(X) 

[0(  '£>)]  +  [0(  3P )]  -+  [0(  3P)]  +  [0(  3P)] 

K-oid  =8x10 

18 

31 

0(‘D)  excites 
02(a) 

[0(  '£>)]  +  [O,  (X)]  -*  [0( 3P )]  +  [02 (a)] 

ka.D  =  1.6  XlO"2  0-6OTg) 

=  1.5  x  1012 

29 

32 

O('D)  excites 
02(b) 

[0(  'D)]  +  [02 (X)]  ->  [0( 3P)]  +  [02 (a)] 

kblD  =  =2.56x10'"  e(-6ing) 

=  2.4x10-" 

29 

15 


Using  this  definition  of  average  energy,  the  temperature  of  the  electrons  in  electron-volts, 
Te ,  is  defined  as: 


Kinetic  rates  are  calculated  from  their  integral  form  shown  in  Equation  2.6  and  are 
dependent  on  the  EEDF.  However,  using  the  EEDF  is  cumbersome,  so  the  electron 
temperature,  as  defined  above,  is  used  to  define  the  system  for  the  sake  of  convenience 
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when  discussing  the  plasma’s  state. 

In  Table  2.1,  all  of  the  reactions  used  in  this  study  are  listed,  and  approximate 
values  are  given.  Reactions  1-16  are  used  in  the  BolSig+  model  to  detennine  the 
EEDF.-  Reaction  1  represents  the  rate  of  elastic  collisions  between  the  oxygen  in  the 
gas  flow  and  electrons.  Even  though  the  rate  coefficient  of  this  reaction  is  large,  the 
energy  exchange  between  the  very  light  electrons  and  heavy  oxygen  molecules  is  very 
small.  As  a  result,  the  reaction  does  not  play  a  significant  role  in  the  energy  balance  of 
the  plasma.  The  reaction  rates  are  detennined  for  a  plasma  maintained  at  a  pressure  of  4 
Torr,  a  flow  velocity  of  20  m/s,  and  an  E-field  operating  at  2.45  GHz.  These  particular 
plasma  conditions  are  relevant  to  the  experimental  discussions  in  Chapter  4.  Reactions 
17-32  will  be  used  throughout  this  document,  and  their  approximate  values  are 
determined  assuming  a  gas  temperature  of  1 ,200  K,  which  is  typical  for  this  study. 

For  a  given  gas  composition  and  field  frequency,  BolSig+  numerically  solves  the 
Boltzmann  equation  to  detennine  the  EEDF,  referencing  tabulated  cross-section  data. 
Using  the  calculated  EEDF,  the  model  then  numerically  integrates  to  determine  the 
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electron  energy  dependent  rate  coefficients  associated  with  reactions  1-16.  The  BolSig+ 
solver  has  been  cited  in  over  35  publications  and  used  to  detennine  the  EEDF  in  a  variety 
of  gas  mixtures,  including  oxygen,  nitrogen,  argon,  and  helium.'  ’  '  The  cross-sections 
included  in  Figure  2.4  are  inputs  to  the  BolSig+  model.  The  measurement  of  these  cross- 
sections  was  reported  by  Phelps,  1985.  The  typical  electron  temperatures  of  the 
discharges  studied  in  this  work  are  approximately  3  eV. 


Figure  2.4  Cross-sections  for  the  production  of  excited  states  of  oxygen  included  in  the  BolSig+ 
modei.  •  elastic  collisions  reaction  1;  o  dissociative  attachment  reaction  2;  ■  vibrational  excitation 
0.02  eV  reaction  3;  □  vibrational  excitation  0.19  eV  reactions  4  and  5;  ♦  vibrational  excitation  0.38  eV 
reactions  6  and  7;  0  vibrational  excitation  0.57  eV  reaction  8;  +  vibrational  excitation  0.75  eV 
reaction  9;  T  02(a)  production  reaction  10;  V  02(b)  production  reaction  11;  AO:(A)  production 
reaction  12;  A  dissociation  reaction  13;  X  dissociation  producing  O^D)  reaction  14;  |  dissociation 
producing  Of'S)  reaction  15;  —  ionization  reaction  16. 


17 


Even  though  the  average  electron  energy  is  3  eV,  reactions  with  energy  thresholds 
larger  than  3  eV,  such  as  ionization,  may  have  significant  rates  of  production,  because  the 
tails  of  the  EEDFs  extend  into  higher  energies.  Some  EEDFs  for  different  E/Ns  are 
shown  in  Figure  2.5.  Typical  E/N  values  for  this  work  range  from  80  to  200  Townsend 
(Td),  which  is  defined  as  1  x  10~17  V-cm.2 

The  average  E/N  values  for  this  study  are  55.6  Td,  which  yields  an  average 
electron  energy  of  3.2  eV.  As  seen  in  the  55.6  Td  case,  there  is  a  significant  population 
of  electrons  with  energies  above  5  eV.  Because  the  electron  populations  with  high 
energies  can  be  large,  the  probability  of  exciting  high-energy  oxygen  states  in  a  plasma 
can  also  become  large. 

It  is  important  to  note  that  this  energy  distribution  would  be  exponential  if 
Boltzmann  statistics  prevailed.  Consequentially,  the  reduced  energy  distribution  would  be 
linear,  on  a  semi-log  plot,  Figure  2.5.  This  difference  in  the  EEDF  demonstrates  the  need 
for  using  the  BolSig+  model  for  this  study.  The  rate  coefficients  in  Figures  2.6  and  2.7 
are  for  the  reactions  included  in  BolSig+.  The  shaded  areas  in  these  figures  are  the 
energy  ranges  of  interest  for  the  two  different  oscillation  frequencies  studied  in  this  work. 
The  rate  coefficients  are  calculated  by  the  integration  shown  in  Equation  2.6.  In 
producing  Figures  2.6  and  2.7,  the  BolSig+  code  was  run  over  a  range  of  E/N  with  values 
from  1  to  250  Td  and  1  -  450  Td  consecutively. 

Introducing  an  oscillating  field  changes  the  EEDF,  if  the  frequency,  a>,  of  the  field 
is  fast  enough  that  the  field  direction  changes  before  the  average  electron  experiences  a 
collision. 
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Electron  Energy,  s,  (eV) 


Figure  2.5  Comparison  of  Maxwellian  distribution  at  Te  =  3.2  eV,  — ,  to  the  electron  energy 
distribution  functions  in  oxygen  with  a  p-wave  field  with  reduced  electric  field  values  of  •  10,  o  55.6, 
and  ■  200  Td.27 


The  approach  to  describing  oscillating  systems  is  to  introduce  an  effective 
electric  field  ( Eeff  )  that  includes  the  effects  of  the  field  frequency7  as  shown  in  and 
MacDonald,  194935: 

E 

where  E  is  the  rms  value  of  the  electric  field,  vc  is  the  collision  frequency,  and  co  is  the 
radian  frequency  of  the  oscillating  field. 
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Figure  2.6  The  reaction  rates  as  a  function  of  reduced  electric  field  for  a  pure  oxygen  discharge  and 
an  RF  field,  13.56  MHz.27  •  attachment  reaction  2;  o  vibrational  excitation  0.02  eV  reaction  3;  ■ 
vibrational  excitation  0.19  eV  reactions  4  and  5;  □  vibrational  excitation  0.38  eV  reactions  6  and  7;  ♦ 
vibrational  excitation  0.57  eV  reaction  8;  0  vibrational  excitation  0.75  eV  reaction  9;  ▼  02(a) 
production  reaction  10;  V  02(b)  production  reaction  11;  A02(A)  production  reaction  12;  A 
dissociation  reaction  13;  +  dissociation  producing  O('D)  reaction  14;  X  dissociation  producing  0(*S) 
reaction  15;  —  ionization  reaction  16. 
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Figure  2.7  The  reaction  rates  as  a  function  of  reduced  electric  field  for  pure  oxygen  gas  in  a  p-wave 
frequency  oscillating  field;  2.45  GHz.27  •  attachment  reaction  2;  o  vibrational  excitation  0.02  eV 
reaction  3;  ■  vibrational  excitation  0.19  eV  reactions  4  and  5;  □  vibrational  excitation  0.38  eV 
reactions  6  and  7;  ♦  vibrational  excitation  0.57  eV  reaction  8;  0  vibrational  excitation  0.75  eV 
reaction  9;  T  02(a)  production  reaction  10;  V  02(b)  production  reaction  11;  AO:(A)  production 
reaction  12;  A  dissociation  reaction  13;  +  dissociation  producing  CH’D)  reaction  14;  X  dissociation 
producing  0{:S)  reaction  15;  —  ionization  reaction  16. 
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Note  that  when  co  is  sufficiently  small  in  comparison  to  the  collision  frequency,  the 
fraction  in  Equation  2.9  approaches  1,  and  hence,  the  effective  electric  field  is  equal  to 
the  rms  value  of  the  electric  field.  ’  At  7  Torr  and  above,  the  collision  frequency  is 
large  enough  that  Eefj  is  not  sufficiently  perturbed  from  the  rms  value  of  the  electric 
field.  However,  at  pressures  lower  than  7  Torr,  the  effects  become  increasingly  stronger. 
In  Figure  2.7,  the  reaction  rates  are  shown  for  a  p-wavc  frequency  electric  field  operating 
at  2.45  GHz  in  4  Torr  of  oxygen.  The  different  ranges  of  Ees/N  for  the  two  different  E- 
field  frequencies  have  a  large  affect  upon  the  ionization  rate  and  will  be  discussed  further 
in  Section  4.8. 

2.4  Energetic  Species  Monitored  within  the  Discharge 

Within  an  electric  discharge,  a  number  of  excited  species  are  created.  Infonnation 
on  plasma  conditions  can  be  extracted  from  observing  the  emissions  from  excited 
species.  The  energy  levels  of  most  of  the  monitored  species  are  shown  in  Figure  2.8.  It 
requires  5.2  eV  to  dissociate  the  oxygen  molecule  into  two  neutral  atoms,  1.6  eV  to  excite 
the  02{b)  state,  and  0.98  eV  to  excite  the  metastable  02(a)  state.  Figure  2.8  also  shows 

that  the  energy  of  the  02(X)  state  and  O  state  are  more  energetically  stable  than  the 
ground  state  of  oxygen. 
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I nternuclear  Spacing,  (nm) 

Figure  2.8  Morse  potentials  representing  the  energy  levels  of  ground  state,  02(a),  02(b),  and  the 
dissociation  energy  of  the  oxygen  molecule. 


2.4.1  Spectral  Emissions  from  the  02(a)  State 

The  transition  of  the  singlet  oxygen  a-state  to  the  ground  state  of  oxygen, 

02  (a,  'A, v'  =  0  ->  X,  3Sg,v"  =  0),  is  centered  at  1 ,268  nm  and  is  spin  and  angular 

momentum  disallowed;  however,  the  transition  is  magnetic  dipole  allowed.  The 
spectrum,  measured  with  a  0.33  meter  spectrometer,  is  shown  in  Figure  4.1 1,  while  the 
possible  branches  for  this  system  are  shown  in  Figure  2.9/  The  expected  transitions  are 
for  AK  =  0,  ±1,  ±2.  K  is  the  quantum  number  related  to  total  angular  momentum, 
excluding  spin,  which  becomes  a  good  quantum  number  under  Hund’s  case  b.  Hund’s 
case  b  predicts  quantum  mechanical  energy  levels  when  the  spin  vector,  S,  is  not  strongly 
coupled  to  the  intra-nuclear  axis.  Therefore,  the  possible  values  of  the  total  angular 
momentum  quantum  number,  J,  are  from  K+S  to  the  absolute  value  of  K-S,  in  integer 
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steps.  The  0-,{a)  state  has  an  angular  momentum  quantum  number  of  two,  and  the  spin 
is  equal  to  zero.  The  subscript,  g,  refers  to  the  inversion  symmetry  of  the  eigenstate  of 

'in 

the  energy  level. 


3  4  5  J(=K) 


Figure  2.9  Branches  of  Magnetic  Dipole  Branches  for  02(a)  emissions  (adapted  from  Herzberg).37 

The  02  {a  A„  — »  3f  3X',  j  transition  has  a  complex  and  dense  rotational  structure, 

with  nine  rotational  branches  (akAJ  =  °P,  pP,  PQ,  ~P,  ~Q,  ~R,  rR,  rQ,  sR).  37  The 
negative  or  broken  circle  transitions  are  asymmetric  states  and  are  absent  from  the 
emission  spectra  due  to  the  zero  nuclear  spin.  The  same  will  be  true  in  the  02  ( b ) 
emissions  in  Figure  2.10. 

2.4.2  Spectral  Emissions  from  the  0?(h)  State 

The  emission  spectrum  for  the  02  (b)  excited  state  to  the  02(X)  ground  state  is 

centered  at  762  mn.  The  total  orbital  angular  momentum  is  equal  to  zero  for  both  states 
involved.  The  spin  of  the  upper  state  is  zero,  while  the  spin  of  the  ground  state  is  one. 
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The  doublets  that  appear  in  the  02(b)  spectra  are  caused  by  the  close  spacing  between 
the  pP  and  PQ  branches  (along  the  P  branch  side  of  the  spectra)  and  the  rR  and  RQ 

'Xl 

branches  (on  the  R  branch  side  of  the  spectra).  These  transitions  are  shown  in 
Figure  2.10.  The  intensity  of  the  emission  lines,  Ib(J ')  in  the  PP  and  PQ  branches,  may 
be  used  to  derive  the  rotational  temperature  and  are  defined  as: 

Ib(J')  =  CSJrN(J' )  (2.10) 


1  2  3  4  5  J(=K) 


Figure  2.10  Possible  branches  for  the  02(b)  emissions  to  the  ground  state  (adapted  from  Herzberg).37 


where 


N(J') 

Xw) 

j' 


fhcB^ 

V  kTR  J 


(2  J '+ 1 )  exp  ( ■ -F(  J  ')hc  /  kTR) 


(2.11) 
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and 


S(j')  =  rotational  line  strength  = 


\  0'+ 1)/2 
(J'+.  25)/ 2 
(J'-.  25)/ 2 
(J'+l)/2 


V 


(2.12) 


A (Jr)  is  the  concentration  of  02  (b,v'  =  0,  J') ;  v  indicates  the  vibrationally  excited 
state;  C  is  an  arbitrary  scale  factor,  including  radiometric  and  detectivity  factors; 
Bv,=0  =  1.39138  cm~l ;  Dv,=0  =  5.486  xlCT6  cm ~l  15’38’39;  F (J1)  =  rotational  energy, 


which  is  equal  to  Bv,=0J'(J'  +  l)  -  Dv,=0  (,/'(./'  +  1))^ ;  and  T R  =  rotational  temperature. 

The  value  of  Tr  may  be  extracted  from  the  measurement  of  the  excited  states’  emissions. 
This  method  is  used  to  determine  the  gas  temperature  in  Chapter  4. 


2.5  Heat  Transport  for  a  Gas  Flow  in  a  Cylindrical  Vessel 

The  velocity  of  the  gas  flow  is  controlled  by  gas  flow  controllers.  The  gas  flow 
controllers  measure  gas  flow  in  standard  cubic  centimeters  per  minute  (SCCM). 
Calculating  the  gas  velocity  from  the  measured  SCCM  is  more  or  less  a  unit  conversion 
from  SCCM  to  cm/s.  Standard  cubic  centimeters  assume  that  the  gas  is  at  standard 
temperature  and  pressure  (STP),  meaning  that  pressure  is  760  Torr  and  the  gas  is  at 
300  K.  A  plug  flow  model  is  used,  ignoring  the  effects  of  the  wall  on  the  velocity.  The 
validity  of  this  assumption  is  based  on  the  flow  having  a  Reynolds  number  of  less  than 
500,  which  is  discussed  in  greater  detail  in  Appendix  A.  So,  the  mapping  from  SCCM  to 
cm/s  is: 


25 


V 


g 


h 

NA 


(0.0432) 


AP 


(2.13) 


where  vg  is  the  gas  velocity;  N  is  the  number  density  of  the  oxygen  in  the  gas  flow 

detennined  by  the  gas  pressure  and  temperature;  A  is  the  area  of  the  flow  tube;  and  h  is 
the  mass  flow  rate  measured  in  SCCM.  The  quantities  on  the  right-hand  side  of 
Equation  2.13  are  all  measured  quantities;  therefore,  velocity  is  a  simple  calculation. 

Flow  velocities  are  calculated  in  Chapter  3. 

Using  a  simple  heat  transfer  model,  the  fraction  of  the  discharge  power  coupled  to 
the  gas  temperature  is  detennined.  Also,  an  experimental  detennination  of  the  Nusselt 
number  is  possible.40  The  Nusselt  number  will  be  discussed  in  greater  detail  in 
Chapter  4. 

As  discussed  in  Appendix  A,  the  energy  equation  for  gas  in  a  gas  flow  may  be 
written  as40: 


pVcp^  =  fpin-shT{Tg-T,)  (2.14) 

Pin  is  the  power  coupled  into  the  gas;  /  is  the  fraction  of  the  power  going  to  gas  heating; 
S  =2  nrl  is  the  surface  area  of  the  cylindrical  tube,  where  the  tube  radius  r  —  0.005  m 
and  the  length  of  the  plasma  /  =  0.025  m;  hr  is  the  heat  transfer  coefficient,  which  varies 
less  than  15%  in  the  temperature  range  of  interest;  p  is  the  mass  density  of  the  gas;  V  is 

the  discharge  volume;  and  C p  =  7.005  kJ/(kg  K),  which  is  the  heat  capacity  of  oxygen.40 

The  pressure  in  the  flow  tube  was  measured  along  the  length  of  the  flow  tube  and 
exhibits  a  less  than  1%  change  in  value.  Therefore,  the  gas  pressure  is  assumed  to  remain 
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constant,  and  any  change  in  the  gas  temperature  will  change  the  gas  density,  but  not  the 
gas  pressure,  as  governed  by  the  ideal  gas  law.  Expressing  the  density  in  tenns  of 
pressure  and  temperature  via  the  ideal  gas  law,  the  equation  may  now  be  written  as: 


m  VCp 


_^_dT^ 
kBTg  dt 


Shr(Tg 


T0)  +  /Pm 


(2.15) 


After  introducing  Pq/P  as  a  unitless  scaling  term,  the  pressure-dependent  data  is 
collapsed  to  a  single  functional  fonn: 


1  dTs  k,ShT  |  kBfP„ 

Tedir  PnmVCp '  s  P0mVCp 

where  kB  is  the  Boltzmann  constant;  P  is  pressure;  m  is  the  mass  of  oxygen;  and 


is  the  pressure  nonnalized  time  of  the  system. 
The  solution  to  Equation  2. 16  is: 


(2.16) 


(2.17) 


Tj  ^ 

(j3  +  T0a  e^') 


(2.18) 


where 


«=  ^  ,  r  =  Toa  +  0,a nd  P  =  (2.19) 

P0mVCp  P0mVCp 

If  it  is  assumed  that  the  temperature  comes  to  a  steady-state,  then  a  steady-state  gas 
temperature  may  be  introduced  to  Equation  2.18.  So,  as  tr  — »  oo  ,  T  (tr)  — >  Tss ,  where 


Tss  is  the  steady-state  temperature  of  the  gas.  In  this  limit,  the  []  tenn  in 
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Equation  2.18  may  be  neglected,  and  Tg  may  be  written  as: 


T  = 


To/  e  r  y 

Tna  er a 


Equation  2.16  may  now  be  rewritten  as: 


_L^ 

Tg  dtr 


a{Tg-T0)  +  p 


(2.20) 


(2.21) 


And  after  multiplying  both  sides  of  the  equation  by  Tg  and  including  the  steady-state 
temperature  defined  in  Equation  2.20: 


T,  dK 


=  a(Tm 


TJ 


(2.22) 


Using  the  new  form  of  the  equation,  the  solution  may  now  be  written  as: 


TA)  = 


TJss  e 


tr 


(T.-Ta)  +  Tt 


tr 


(2.23) 


This  functional  fonn  of  the  time-dependent  temperature  applies  to  the  temperature  profile 
within  the  electric  discharge  discussed  in  this  work.  The  Nusselt,  Nu,  number  describes 
the  relationship  between  convective  heat  transfer  and  conductive  heat  transfer.  Nusselt 
values  greater  than  100  are  indicative  of  turbulent  flows.  Therefore  a  value  significantly 
less  than  100  is  suggestive  of  a  flow  operating  in  the  laminar  regime.  For  a  cylindrical 


flow  tube  24 : 


hrD  . 

Nu  =  -d—  (2.24) 

kT 

where  D  is  the  diameter  of  the  flow  tube,  kT  is  the  thennal  conductivity  of  the  fluid,  and 
hT  is  the  convective  heat  transfer  coefficient.  For  all  cylindrical  flow  tubes,  the  Nusselt 
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number  should  equal  approximately  4. 3 64. 40  This  analysis  will  play  an  important  role  in 
determining  the  time-dependent  temperature  of  the  gas  flow. 


In  Figure  2.11,  the  functional  fonn  of  the  gas  temperature  (using  Equation  2. 18)  is 
shown  for  oxygen  gas  flowing  in  a  cylindrical  tube  and  being  heated  via  a  capacitively 
coupled  100  W  discharge.  The  x-axis  in  Figure  2.11  is  the  amount  of  time  required  for 
the  gas  to  travel  through  the  discharge  region.  This  amount  of  time  is  nominally  the 
length  of  the  discharge  divided  by  the  velocity  of  the  gas.  Residence  time  will  be 
discussed  in  greater  detail  in  Section  4.2.  The  fraction  of  power  coupled  into  gas  heating 
is  arbitrarily  chosen  to  be  20%.  As  calculated  in  Appendix  B,  a  significant  amount  of 
power  is  coupled  into  the  formation  of  excited  states  within  the  gas  flow;  therefore,  all 
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the  power  does  not  need  to  be  coupled  into  gas  heating  or  lost  via  radiation.  There  are 
two  important  parameters  to  be  extracted  from  this  functional  form.  First  is  the  initial 
slope  of  the  curve,  which  is  determined  primarily  by  the  strength  of  the  heating  tenn. 
Second  is  the  steady-state  temperature  of  the  gas,  which  is  the  balance  between  the 
heating  term  and  heat-loss  tenn. 

2.6  Determination  of  O-atom  Concentrations  Using  Actinometry 

The  actinometry  technique  uses  the  emissions  from  the  Of 3.v  3  S  — »  3 p  P) 
transition  at  844  mn  and  the  Ar(4s  — »  4 p)  transition  at  750  nrn  to  detennine  the  O-atom 
number  density  within  an  oxygen  plasma.41  Energy  level  diagrams  for  these  two  states 
are  shown  in  Figure  2.12. 


Figure  2.12  Schematic  drawing  of  the  energy  levels  involved  in  determining  the  concentration  of  the 
ground  state  oxygen  atoms  in  the  actinometery  method. 


It  is  assumed  that  the  two  main  pathways  of  production  of  the  0(  P)  state  are 

38 

direct  excitation  by  electrons  and  dissociative  excitation,  as  shown  below  : 


30 


e~  +0(2  p  3P) 
e~  +  02(X,3Z) 


(2.25) 


0(3p3P)  +  e- 
0(2 p  3P)  +  0(3 p  3P)  +  e~ 

where  k°  is  the  rate  coefficient  for  excitation  of  O-atoms  from  the  ground  state  to  the 
0(  '  P)  state,  and  kde  is  the  rate  coefficient  for  dissociation  of  molecular  oxygen  to  two 
O-atoms  [one  excited  to  the  0(  P )  state]. 

It  is  assumed  that  the  mechanism  for  excitation  of  Ar(2px )  is  direct  electron 

impact: 

kAr 

e  +Ar(\p )  — - — >  Ar(2pl)  +  e  (2.26) 

where  kf  is  the  rate  of  excitation  from  Ar(\p  — »  2px)  .  As  shown  in  Figure  2.13,  the 

excitation  cross-sections  for  the  Ar( 2 px )  and  0(  P)  states  are  similar,  which  allows  this 

method  to  detennine  the  O-atom  population  from  the  argon  emissions.  The  primary  loss 
mechanism  for  these  two  excited  states  is  radiative  : 


0(3 p  3 P)  >  0(3s  3S)  +  hv  (844. 6«m) 

Ar(2p\)  — IAa,:  ->  Ar(2p  )  +  hv  (150 Anm) 


(2.27) 


where  ZAa  and  IAAr  are  the  sum  of  all  the  radiative  loss  rates  of  the  two  excited  states. 
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1  10  100 


e  I#V  1 

Figure  2.13  Cross-sections  for  electronic  excitation  of  the  0(3P)  and  Ar(2p0  states  and  reduced 
electron  energy  distribution  at  30  and  100  Td.38  Image  taken  from  Pagnon  et  al. 

Collisional  relaxation  with  the  ground  state  will  slightly  influence  the  actinometry  results, 

38 

so  these  reactions  are  accounted  for  in  this  analysis  : 

0(3p3p)+a(x)  ~  *  0(2/?  3ma(x) 

A  (2.28) 

Ar(2p\)  +  02{X)  —  >  Ar(2 p)  +  02(X) 

where  kg  is  the  rate  coefficient  for  collisional  relaxation  of  the  0(  ' P)  state  by  ground 
state  oxygen,  and  ^'  is  the  rate  coefficient  for  collisional  relaxation  of  the  Ar(2p{ )  state 
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by  collision  with  ground  state  oxygen.  With  all  of  these  reactions  considered,  the 

38 

intensity  of  the  argon  emissions  and  oxygen  atom  emissions  are  : 


lAr}kf 

S4,,  +  [<W  K 


([0}k° +K,[Oi(x)]) 

24  +  [o2(X)]^ 


(2.29) 


where  I  a,-  is  the  intensity  of  the  argon  emissions,  Io  is  the  intensity  of  the  O-atoms,  A  Ar  is 
the  Einstein  A  coefficient  for  argon,  and  Ao  is  the  Einstein  A  coefficient  for  O-atoms. 
Therefore,  the  ratio  of  the  measured  emissions  resulting  from  these  two  excited  state 
takes  the  fonn  of: 


K  hv„A0  [0;U)]^^,  +  [02W]^ 
I  a,  l>vmAAr  [ Ar]k f  ZA,  +  [0,(X)]k°Q 


This  relationship  may  also  be  written  as: 


L  1  f°1 

I  A,  Cl[02(X)] 


where 38 

c„  =  hv„„A„  ( [02(X)]<:Ay1  r ZA4r+[Q2(X)] k£' 
"  hvMA0{  [ Ar]k f  J  {T.A,+[02(X)]k°  y 


(2.30) 


(2.31) 


(2.32) 


From  Equation  2.3 1,  it  may  be  seen  that  the  ratio  of  the  measured  intensity  of  the 
844  mn  line  and  the  750  mn  line  is  a  constant  value,  C°r,  multiplied  by  the  ratio  of 
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[O]/[O2(^0]  •  The  constant,  as  determined  by  Pagnon,38  is  2.1  ±  0.15  x  10'3  with  a  7% 
uncertainty  for  a  given  ratio  of  [0]/[02(X)]  .38  The  Pagnon  results  are  measured  over  a 
pressure  range  of  0.36  to  2  Torr  and  a  discharge  current  of  5  to  80  mA.  Note  that  CAr  is 
a  function  of  the  ratio  of  argon  and  molecular  oxygen,  as  well  as  the  Eeff/N  of  the 
discharge.  The  ratio  of  Ar/O^  in  the  Pagnon  study  is  10'  .  In  this  study,  the  ratio  of 

argon  to  oxygen  is  0.1.  Therefore,  the  value  of  CAr  is  scaled  by  an  order  of  magnitude  to 
2. 1  x  10'“,  as  dictated  by  the  increase  of  the  argon  concentration.  The  determined  values 
of  CC-r  have  an  approximate  20%  error  associated  with  them.  The  actinometry  method  is 

well  understood  and  widely  applied  for  the  detection  of  O-atoms  in  the  presence  of 
electrons.42'54 

2. 7  Operational  EOIL  System 

Bernard  and  Pchelkin  explored  the  production  of  02  (a)  within  an  electric 
discharge  in  19787  Using  a  microwave  discharge  at  70  W  and  450  mTorr,  they  reported 
a  maximum  O^a)  yield  of  1 1%.  This  value  is  short  of  the  15%  threshold  for  gain  at 

1  T 

room  temperature.  Also,  work  released  by  Schmiedberger,  Takahashi,  and  Fuji  in  2001 
reports  a  maximum  02(a)  yield  of  30%  over  the  E-field  frequency  range  of  25  to  100 
MHz.  With  a  30%  yield,  there  should  be  a  significant  amount  of  extractable  power; 
however,  the  reported  maximum  laser  power  is  approximately  2  p-watts,  which  is  similar 
to  power  levels  expected  from  florescence.  The  calibrated  high-purity  gennanium 
(HPGe)  detector,  which  they  use  to  detennine  the  T  ,  does  not  spectrally  resolve  the 
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emissions.  An  unresolved  spectrum  is  susceptible  to  corruption  if  there  is  another 
emitting  species  in  the  detected  frequency  range.  Including  emission  intensities  from 
other  excited  species,  such  as  O-atoms,  may  explain  the  discrepancy  between  the  reported 
02(a)  yield  and  laser  power. 

2. 7.1  Work  at  the  University  of  Illinois  at  Urbana-Champaign 

In  2004,  UIUC  demonstrated  the  first  operation  of  an  EOIL  device. 5  A 
schematic  of  the  system  is  shown  in  Figure  2. 14.  The  gas  flow  is  made  up  of  oxygen; 
helium,  which  is  the  carrier  gas  that  comprises  83.55%  of  the  gas  flow;  and  trace  amounts 
of  NO  in  the  flow  (0.78%).  NO  and  NO2  are  parasitic  to  O-atoms  via  a  catalytic  reaction. 
Therefore,  including  NO  in  the  discharge  decreases  the  amount  of  O-atoms  that  result 
from  the  discharge,  while  also  serving  as  an  electron  donor  to  allow  the  discharge  to  be 
self-sustaining  at  an  E/N  of  10  Td.  This  E/N  is  ideal  for  the  production  of  02(a )  ,  as 
stated  by  Stafford  and  Kushner.7  Introducing  NOo  in  the  gas  flow  downstream  of  the 
discharge  serves  two  purposes  according  to  the  authors.  Primarily  NO2  will  further 
reduce  the  amount  of  O-atoms  in  the  gas  flow,  by  the  catalytic  reaction.  Also  the 
addition  of  a  room  temperature  gas  downstream  of  the  discharge  will  cool  the  hot  effluent 
from  the  discharge  region.  O-atoms  quench  the  excited  state  of  iodine  responsible  for 
lasing  in  the  system,  so  reducing  O-atom  concentrations  increases  the  steady-state  value 
of  the  excited  iodine  state.7  It  is  reported  that,  even  though  some  O-atoms  are  helpful  in 
dissociating  E,  too  many  will  deactivate  I*,  which  has  a  detrimental  effect  on  laser 

-5 

perfonnance.  In  2005,  a  maximum  02{a )  yield  of  17%  was  reported  without  the 
inclusion  of  NO2  into  the  gas  flow.  The  maximum  laser  power  of  the  system  was  520 
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mW  with  an  RF  discharge  coupling  400  W  into  the  discharge.  The  operating  conditions 
of  the  laser  were  a  total  pressure  of  12.9  Torr  with  oxygen  comprising  15.67%  of  the  gas 
flow  or  a  partial  pressure  of  approximately  2  Torr. 


Diagnostic  Duct  Lasing  Cavity 


Figure  2.14  Schematic  design  of  an  electric  discharge  pumped  oxygen-iodine  laser  demonstrated  by 
the  University  of  Illinois  at  Urbana-Champaign.  Partial  pressure  of  oxygen  was  2  Torr  in  the  flow 
and  520  mW  laser  power  coupled  out  of  the  system. 3 


In  January  2008,  UIUC  applied  a  100  W  RF  discharge  to  predissociate  the  1 2 
before  it  was  injected  into  the  gas  flow.55  Operating  at  a  total  pressure  of  300  Torr,  a 
mixture  of  He:02:NO  at  50: 10:0.05  was  passed  through  a  700  W  RF  discharge.  The  pre¬ 
dissociation  of  the  F  resulted  in  a  gain  of  0.067%  cm'1  and  a  maximum  laser  power  of  6 
W.  This  report  was  followed  closely  by  a  report  in  which  the  primary  discharge  was 
changed  from  a  single  discharge  in  a  gas  flow  tube  with  a  radius  of  4.9  cm  to  two  flow 
tubes  with  radii  of  1.6  cm.  Each  flow  tube  had  power  capacitively  coupled  into  it  using 
two  1,000  W  RF  discharges.  Once  again,  the  F  was  pre-dissociated,  using  a  100  W  RF 
discharge;  the  reported  gain  was  0.17%  cm'1  with  a  maximum  power  of  12  W.4 
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2. 7.2  Work  by  Alan  Hill  at  Kirtland  AFB 

Alan  Hill  proposed  a  method  for  producing  02{a )  with  higher  energy 
efficiency.55  The  major  problems  with  electrically  driven  SOGs  are  that  self-sustained 
discharges,  operate  at  reduced  electric  fields  that  do  not  produce  02(a )  efficiently,  20  - 
40  Td.  Therefore  a  large  portion  of  the  energy  coupled  into  the  discharge  goes  to  the 
creation  of  excited  species  which  are  not  02  ( a  )  and  instead  result  in  heating  the  gas 

flow.  Heating  the  gas  is  detrimental  to  an  EOIL  system  in  two  ways.  First,  any  power 
that  is  coupled  into  the  temperature  of  the  gas  is  not  available  for  the  production  of 
02(a) .  Second,  the  pooling  reaction, 

[02(a)\  +  [02(a)\  ->  [02(X)]  +  [02(b)\ ,  (2.33) 

increases  exponentially  with  gas  temperature.7  So,  not  only  is  the  energy  used  to  heat  the 
gas-wasted  power,  but  it  also  decreases  the  amount  of  02(a )  created  in  these  systems.  It 
is  possible  for  an  outside  source  of  electrons  to  sustain  the  discharge.  Therefore,  the  E- 
field  could  be  held  at  low  levels  that  are  ideal  for  the  production  of  02(a)  ,  and  this 

would  also  keep  the  gas  temperature  low.  Current  results,  using  pre-ionization  by  UV 
radiation,  with  voltage  pulses  of -150  Td,  -10-30  ns  pulses  at  50,000  pulses  s'1,  and  a  DC 
E/N  of -10  Td  promise:  (1)  singlet  oxygen  yields  of  30%;  (2)  an  electric  excitation 
efficiency  of  40%;  (3)  specific  energy  loading  of  150  kJ/m  O2;  and  (4)  a  temperature 
increase  of  less  than  125  K.  Lasing  from  this  system  has  not  been  reported  yet.57 
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2. 7.3  Work  at  The  Ohio  State  University 

Work  at  The  Ohio  State  University  (OSU)  is  based  on  using  the  pulser-sustainer 
method  of  maintaining  a  discharge  first  proposed  by  Hill.  ’  Typical  glow  discharges 
have  two  drawbacks,  which  are  addressed  by  using  pulser-sustainer  systems.  One  is  that 
the  supersonic  nozzles  employed  in  EOIL  systems  require  significant  gas  pressures 
upstream  of  the  supersonic  nozzle.  Typical  electric  discharges  in  oxygen  are  not  stable  at 
these  pressures.  Pulser-sustainer  methods  establish  stable  discharges  at  high  pressures. 
Secondly,  the  optimal  E/N  for  the  production  of  the  02(a)  state  is  approximately  10  Td, 

which  is  significantly  less  than  the  E/N  in  which  a  typical  oxygen  discharge  is 
sustainable.7  By  combining  high-voltage  pulses  with  a  constant  field  of  relatively  low 
voltage,  a  pulser-sustainer  is  able  to  maintain  an  average  E/N  around  10  Td.7  The 
frequency  and  power  of  the  sustainer  pulse  is  tailored  to  keep  the  discharge  stable.  In 
January  2007,  OSU  demonstrated  a  peak  laser  power  of  0.28  W,  using  a  pulser-sustainer 
SOG.57  The  gas  flow  was  15%  02/85%  He  at  a  discharge  pressure  of  60  Torr  and  a 
discharge  power  of  1.55  kW.  A  maximum  gain  of  0.04%  cm'1  was  reported  at  the  same 
conditions  but  at  an  80  Torr  discharge  pressure.  In  June  2007,  OSU  implemented 
UIUC’s  technique  and  included  NO  in  the  gas  flow  through  the  discharge.  Maximum 
laser  power  was  reported  at  1.24  W  with  a  gain  of  0.049%  cm'1.10  The  gas  flow  was 
reported  as  a  15%  02/He  mixture  with  190  ppm  I2  and  550  ppm  NO  at  a  discharge 
pressure  of  107  Torr  and  a  discharge  power  of  2.4  kW.  Both  groups  reported  an  increase 
in  laser  power  by  including  NO  in  the  gas  flow  through  the  discharge.  Because  the  E/N 
in  pulser-sustainer  methods  is  already  ideal  for  the  production  of  02(a)  ,  the  increase  in 
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laser  power  may  be  attributed  to  the  reduction  in  the  quenching  effect  that  O-atoms  have 
on  the  I*  density. 

2.8  Scaling  EOIL  System 

UIUC  has  shown  that  an  EOIL  will  operate  at  laser  powers  up  to  102.5  W.64  The 
next  step  in  the  evolution  of  the  EOIL  weapon  system  is  to  show  a  path  to  scaling  from 
100  W  to  kW,  as  well  as  higher  output  energy  systems.  In  oxygen-iodine  laser  systems 
the  extractable  power  is  detennined  by  Ya ,  as  well  as  the  number  of  02(a)  molecules  in 
the  system.  In  order  to  scale  laser  power  from  the  current  system  to  a  feasible  weapon 
system,  the  oxygen  pressure  in  the  cavity  must  be  increased,  while  Ya  is  held  constant  or 

increased.  The  historic  results  of  this  system  suggest  that  accomplishing  this  power 
scaling  may  be  problematic. 

Figure  2.15  shows  the  results  of  both  experiments,  as  well  as  computer  models 
reported  by  many  different  research  groups.  A  large  parameter  space  is  explored  in  these 
studies,  including  different  buffer  gases,  E-field  frequency  oscillations,  coupled  discharge 

2  3  12  56  59  63 

power,  and  discharge  geometries.  ’  ’  ’  ’ 

The  general  trend  of  the  measured  results  shows  a  decrease  in  the  achievable 
02(a)  yield  as  a  function  of  the  partial  pressure  of  oxygen.  There  are  three  points  shown 

in  Figure  2.14  that  break  this  trend.  They  are  the  results  of  A  Hill,57  O  Ionin,62  and  □ 
Naparitovich.61  These  three  values  are  the  result  of  computer  models,  implying  that  there 
may  be  a  discrepancy  between  simulated  and  measured  kinetics,  which  detennine  the 
value  of  02(a)  that  is  achievable  within  an  oxygen  glow  discharge.  Understanding  this 
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difference  is  important  to  determining  the  feasibility  of  using  an  electrically  driven  SOG 
within  a  weapon  system. 

The  research  group  at  UIUC  has  scaled  the  output  power  of  an  EOIL  device  by 
two  orders  of  magnitude  since  first  light  was  achieved  in  2005. 64, 65  A  number  of 
different  modifications  have  been  made  to  the  lasing  system  in  order  to  increase  output 
laser  power.  One  significant  change  made  by  the  research  group  at  UIUC  is  the  number 
of  flow  tubes  passing  through  the  discharge  region.  All  of  the  flow  tubes  are  combined 
after  the  discharge  region  before  the  iodine  is  injected  into  the  gas  flow. 

By  increasing  the  number  of  flow  tubes,  from  one  in  200855  to  six  in  20 10,64  the 
group  has  increased  the  volume  of  the  laser  resonant  cavity,  while  maintaining  a  constant 
gain  of  approximately  0.26%  cm'1.  It  is  of  particular  note  that  the  number  of  discharge 
regions  has  been  scaled,  while  holding  the  dimensions  of  the  discharge  region  constant. 
This  suggests  the  importance  of  the  surface  area  to  discharge  volume  and  may  be  linked 
to  gas  kinetics  governing  the  concentration  of  02(a)  within  the  discharge. 

Recent  publications  have  expressed  the  need  to  completely  define  the  chemical 
kinetics  that  dominate  excited  state  populations  for  EOIL  devices.65  Zimmerman  et  al. 
published  a  study  on  the  effects  of  frequency  and  discharge  geometry  upon  the  discharge 
effluent.66  It  was  found  that  the  conditions  favorable  for  02(a)  production  involve 

reducing  the  flow  tube  diameter,  while  increasing  the  residence  time  to  limit  the  O-atom 
concentrations. 
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Figure  2.15  ■  Schmiedberg(RF)12;  AHill(Controlled  Avalanch)57;  ARakimova(RF)17;  O  Ionin62;  □ 
Naparitovich61;  Pchelkin(1978  MW)2;  Itami(1999  MW)";  Savin63;  x  Carrol(RF)3;  +  Rawlins(2005 

MW)59;  •  Rakimova16;  ♦  super-elastic  theoretical  limit.16  Historical  plot  reporting  the  results  of  a  number 
of  different  groups  operating  under  various  different  pressures,  power  frequencies,  discharge  geometries.  Both 
measured  data  and  simulated  data  are  reported. 


By  reducing  the  radius  of  the  flow  tubes  the  surface  area  to  volume  ratio  of  the 
discharge  will  be  increased.  Therefore  the  concentrations  of  any  excited  species  that  is 
primarily  destroyed  by  collisions  at  the  wall  will  be  reduced.  A  paper  by  Davis  et  al. 
discusses  power  scaling  in  a  1  to  5  kW  supersonic  discharge-flow  reactor.  While 
coupling  1  to  5  kW  of  power  via  a  microwave  discharge  into  a  gas  flow,  it  was  found  that 
the  quenching  of  I*  by  O-atoms  is  not  fast  enough  to  explain  the  loss  of  excited  species 
within  the  gas  flow.  This  result  points  to  a  destruction  mechanism  which  is  other  than 
direct  collisions  with  O-atoms.  Lee  et  al.  presented  a  study  of  catalytically  produced 
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01  (a)  by  using  the  flux  of  O-atoms  that  result  from  the  discharge.  In  this  work  an 

iodine  oxide  catalyst  is  applied  to  the  walls  of  the  flow  tube  downstream  of  the  discharge 
region.  Results  suggest  that  the  O-atoms  created  within  the  discharge  region  react  at  the 
walls  producing  02{a )  .  Other  excited  species,  such  as  vibrationally  or  electrically 
excited  O2  and  O3,  may  also  be  involved  in  the  catalysis  at  the  wall.  A  2010 
publication  by  Azyazov  et  al.69  suggested  that  the  production  of  vibrationally  excited 
ozone  is  significant  in  plasmas  with  an  excess  of  ground  state  oxygen.  This  study  argues 
that  including  this  mechanism  renders  including  the  three-body  reaction  (reaction  24) 
unnecessary  to  reproduce  measured  [02(o)] .  Recent  publications  by  the  EOIL 
community  have  established  the  need  for  a  greater  understanding  of  the  governing 
kinetics  involved  in  EOIL  devices  to  allow  scaling  to  kW  laser  powers. 
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Ill  Equipment 


The  apparatus  is  designed  to  study  the  behavior  of  an  electric  discharge-driven 
SOG.  A  glow  discharge  is  created  by  flowing  pure  oxygen  contained  in  a  glass  flow  tube 
through  an  intense  electric  field.  Oscillating  electric  fields  are  used  at  a  microwave 
frequency  of  2.45  GHz  and  a  radio  frequency  of  13.56  MHz.  Emissions  from  the  excited 
states  of  02(a )  ,  02(b) ,  and  various  O-atom  emissions  are  used  to  characterize  the 
discharge.  Absolute  concentrations  of  O-atoms  are  detennined  using  the  actinometry 
technique  by  comparing  O-atom  emissions  to  atomic  argon  emissions. 


3.1  Microwave  Discharge 

A  microwave  discharge  operating  at  2.45  GHz  sustains  the  glow  discharge  of  a 
pure  oxygen  gas  flow  in  a  Pyrex  tube.  This  apparatus  is  shown  in  Figure  3.1. 
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The  flow  tube  has  an  outer  diameter  of  1  cm  and  an  inner  diameter  of  0.8  cm.  With  the 


use  of  a  bifurcated  optical  bundle,  the  emissions  from  both  singlet  states  are  collected 
simultaneously  at  a  given  spatial  location  within  the  discharge  and  in  the  discharge 
afterglow.  The  fiber  is  located  at  the  center  of  the  discharge  region  and  emissions  are 
collected  for  an  integration  time  of  two  minutes.  The  observation  region  of  the 
collimated  optics  is  approximately  one  cm  in  diameter.  The  discharge  region  in  this 
study  is  also  one  cm  in  diameter  and  averages  about  3  cm  long;  therefore,  the  optics 
gather  light  from  roughly  one  third  of  the  discharge  region.  In  order  to  measure  the 
changes  in  the  excited  species  as  a  function  of  residence  time  within  the  discharge  the 
flow  velocity  of  the  gas  is  varied.  The  emissions  are  spectrally  resolved,  and  the  intensity 
is  detennined  by  summing  the  emissions  for  each  wavelength  range.  For  the  02(b) 

emissions,  the  wavelength  range  is  761  mn  to  768  mn,  and  for  02  ( a )  is  1,240  nm  to 

1,275  nm.  Two  spectrometers  are  employed. 

A  Roper  Scientific  InGaAs  detector  array,  attached  to  a  0.33  meter  McPherson 
spectrometer  equipped  with  a  600  grooves/mm  grating  blazing  at  1,250  nm,  is  used  to 
observe  the  emissions  from  02(a)  centered  at  1,272  nm.  The  slit  width  of  the 

monochromator  is  250  pm,  resulting  in  a  spectral  resolution  of  approximately  1.7  nm. 

The  InGaAs  array  is  cryogenically  cooled  to  approximately  -100°  C.  Also,  a  Roper 
Scientific  intensified  charged-coupled  device  (CCD)  array,  attached  to  a  1.33  meter 
McPherson  spectrometer  equipped  with  a  2,400  grooves/mm  grating  blazing  at  700  nm, 
is  used  to  observe  the  02(b)  emissions  centered  at  762  nm.  The  slit  width  is  10  pm, 
which  yields  a  resolution  of  approximately  0.01  nm.  The  spectral  resolution  of  the 
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emissions  allows  the  removal  of  interference  of  other  emitting  species  within  the  spectral 
band  of  the  excited  state  emissions.  The  total  intensity  of  the  emissions  is  detennined  by 
summing  the  molecular  emissions  of  the  appropriate  frequency  range:  761  to  768  nm  for 
02(b )  and  1,240  to  1,275  nm  for  02(a)  . 


1.20  1.22  1.24  1.26  1.28  1.30 

Wavelength  (pm) 


Figure  3.2  Comparison  of  the  (— )  theoretical  black  body  D  IK  to  the  (— )  experimental  black  body. 


The  detectivity  of  both  detectors  is  a  function  of  the  wavelength  of  the  light  being 
detected.  This  dependence  is  determined  using  a  black  body  light  source  for  calibration. 
The  emissions  of  a  1,000°  K  black  body  are  measured  and  compared  to  the  theoretical 
intensities.  The  gratings  are  positioned  to  center  the  desired  wavelengths  upon  the 
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detector  arrays.  The  emissions  from  the  black  body,  as  detected  by  the  InGaAs  array,  are 
shown  in  Figure  3.2. 

The  theoretical  calculation  includes  an  uncertainty  of  plus  or  minus  1°  K.  The 
ratios  of  the  measured  and  theoretical  values  are  the  detectivity  factor,  shown  in 
Figure  3.3. 


Figure  3.3  The  detectivity  factor  of  the  InGaAs  array  used  to  measure  the  02(a)  emissions  from  the 
discharge. 


The  dark  current  of  the  InGaAs  array  has  structure;  therefore,  it  is  important  to  be 
careful  when  subtracting  the  dark  current  signal  from  the  measured  emissions.  The  dark 
current  is  measured  over  the  two-minute  integration  time  and  subtracted  from  the 
measurements,  as  shown  in  Figure  3.4.  The  pixel-by-pixel  change  in  dark  current  is  not 
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wavelength-dependent  and  is  confirmed  to  be  typical  by  the  manufacturer  of  the  detector 


array. 


Figure  3.4  Subtraction  of  dark  current  from  the  measured  intensity.  —  raw  data  data  after 
background  subtraction. 


An  isolated  atomic  oxygen  peak  is  used  to  calibrate  the  resolution  of  the 
spectrometer.  In  Figure  3.5,  the  resolution  of  the  system  is  determined  by  fitting  a  sum  of 
Gaussian  peaks  to  the  measured  emission  from  the  F-  D  atomic  emission  at  1,257  nm. 
The  asymmetry  in  the  peak  is  due  to  a  slight  misalignment  of  the  detector  and  grating. 

The  calibration  of  the  pixel  number  to  wavelength  is  established  by  the  spectrum  itself. 
The  fit  for  the  1.33  meter  spectrometer  is  0.00717  nm/pixel,  and  the  fit  for  the  0.33  meter 
spectrometer  is  1 .925  nm/pixel. 
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Figure  3.5  The  5F-3D  atomic  emission  at  1,257  nm  is  used  to  determine  the  resolution  of  the  0.33  m 
monochromator  used  to  monitor  02(a)  emissions.  Resolution  is  determined  to  be  0.15  nm.  Peak  is  a 
sum  of  two  Gaussian  peaks. 


Emissions  from  both  the  02(a,b)  states  are  observed  simultaneously,  using  a 
bifurcated  fiber  bundle.  A  1  cm  diameter  collimating  lens  is  positioned  on  the  end  of  the 
bundle.  The  acceptance  angle  of  the  fiber  is  68  degrees.  Combined  with  the  aperture 
stop,  this  determines  a  spatial  resolution  of  0.5  cm.  The  fiber  bundle  is  positioned 
perpendicular  to  the  flow  tube  and  mounted  on  a  sliding  mount,  allowing  for  parallel 
displacement  along  the  flow  tube.  The  focusing  lens  is  located  1.5  mm  from  the  edge  of 
the  Pyrex  flow  tube  to  maximize  the  spatial  resolution. 

The  physical  dimensions  of  the  discharge  are  determined  by  the  point  at  which  the 
plasma  begins  to  glow  and  where  it  stops  glowing.  The  discharge  glow  is  used  as  an 
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indication  of  electron  density,  and  the  point  at  which  it  tenninates  indicates  the  lack  of 
electrons  in  the  gas  flow.  The  length  of  the  discharge  plays  a  role  in  the  power  loading  as 
well  as  the  residence  time  of  the  gas. 

The  plasma  typically  terminates  a  couple  of  millimeters  outside  of  the  p-wave 
cavity;  however,  at  low  pressures  and  flow  velocities,  the  length  of  the  discharge  region 
extends  beyond  the  dimensions  of  the  cavity.  The  variation  in  the  plasma  dimension  as  a 
function  of  mass  flow  rate  is  shown  in  Figure  3.6  for  a  number  of  different  pressures. 


Figure  3.6  The  length  of  the  discharge  as  mass  flow  rate  changes,  o  2  Torr,  •  3  Torr,  □  4  Torr, 
5.2  Torr,  A  6.4  Torr,  ▲  7.5  Torr,  V  8.6  Torr. 


The  flow  rate  of  the  input  gases  is  controlled  using  a  Sierra  Series  840  Side-Trak  Mass 


Flow  Controller  calibrated  to  control  nitrogen  from  0  to  2,000  SCCM  with  a  resolution  of 
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8  SCCM.  For  the  calculation  of  residence  for  the  temperature  measurements  included  in 
this  document  these  pressure  and  flow  rate  lengths  are  used.  For  the  discharge  volume 
calculation  included  in  Chapter  4,  an  average  value  for  discharge  length  of  3  cm  is  used 
in  for  calculations  throughout  this  work. 

The  pressure  is  measured  using  an  MKS,  Model  690A,  10  Torr,  Capacitance 
Manometer.  The  lowest  achievable  pressure  is  less  than  1x10'  Torr,  which  is  beyond  the 
range  of  the  head.  The  vacuum  pump  is  an  Alcatel  2063 C2  rotary  vane  pump.  An 
Evenson  cavity  operating  at  2.45  GHz  with  a  range  of  0  to  100  W  controlled  by  an 
Opthos  power  supply  is  employed  to  produce  the  glow  discharge  in  the  gas  flow. 

Figure  3.7  shows  the  velocity  versus  the  mass  flow  rate  of  the  gas  at  a  pressure  of 

4  Torr. 


Figure  3.7  Velocity  as  a  function  of  mass  flow  rate  at  4  Torr. 
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Figure  3.7  shows  that  the  gas  velocity  ranges  from  800  to  15,000  cm/s,  while  the  mass 
flow  rate  ranges  from  2  to  2,000  SCCM.  The  gas  flow  is  adjusted  by  the  throttle  valve  to 
maintain  a  certain  pressure  (in  the  case  of  Figure  3.7,  4  Torr). 

3.2  RF  Discharge 

Using  an  RF  electric  field  to  maintain  the  electric  discharge  is  very  similar  to  the 
microwave  approach.  The  field  is  coupled  into  the  gas  by  two  copper  rings  positioned  10 
cm  apart  with  a  width  of  1  cm.  The  RF  generator  is  a  Dressier  Cesar  1310  operating  at 
13.56  MHz  with  a  power  ranging  from  0  to  1,000  W.  The  RF  power  generator  is 
matched  to  the  plasma  flow  tube  using  a  VarioMatch  matching  network  supplied  by 
Dressier.  The  auto-matching  network  is  calibrated  for  frequency  ranges  from  13  to  27 
MHz  and  up  to  1  to  1 .5  kW,  depending  on  the  frequency  of  the  field  oscillations.  A 
schematic  of  the  RF  discharge  experiment  is  shown  in  Figure  3.8. 


Figure  3.8  Schematic  of  the  RF  glow  discharge  apparatus. 


51 


The  1.33  and  0.33  meter  monochromators  used  in  the  RF  experimental  apparatus 
are  the  same  as  those  used  in  the  microwave  experiment.  Therefore,  their  calibration  and 
resolution  are  measured  and  perfonned  with  the  same  results  as  reported  in  Section  3.1. 
The  argon  inserted  into  the  gas  flow  is  used  to  monitor  the  O-atom  concentrations  via  the 
actinometry  technique,  which  was  explained  in  Section  2.6. 
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IV  Results 


4.1  02(b)  Spectrum  and  Gas  Temperature  Measurement 

The  rotational  spectra  of  the  (92(6‘E  — »  X3Z)  (0,0)  and  (1,1)  bands  near  762  mn 
are  shown  in  Figure  4.1.  Only  the  even  number  J  values  are  observed  due  to  the  selection 
rules  discussed  in  Section  2.4.2.  The  spectral  resolution  of  the  reported  02(b)  emissions 
is  slightly  better  than  the  results  currently  in  the  literature,  which  generally  have  a 
spectral  resolution  of  0.04  mn  or  greater.15 


Figure  4.1  The  emission  spectrum  of  the  02(b)  excited  state,  inside  the  p-wave  discharge,  4  Torr,  flow  velocity  of 
20  m/s,  measured  using  the  1.33  meter  spectrometer. 
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This  added  resolution  permits  the  observation  of  the  atomic  emissions  from  the  O-atom 
ion  transition  ( 4D-4D° )  at  766  nm;  O-atom  ion  transition  (3P-3S°)  at  111  nrn;  and  02 ( b ) 

(1-1)  band  centered  at  77 1  nm.  Clear  identification  of  these  emissions  permits  the 
extraction  of  their  contributions  from  the  measured  intensities  and  results  in  a  greater 
certainty  of  the  temperature  measurements  reported  in  this  work. 

The  rQ,  rR,  pQ,  and  PP  branches  are  readily  observed  and  assigned  in  Figure  4.1. 
The  emissions  from  ()( 5 S  — »  P)  around  777  nm  are  observed  at  the  far  end  of  this 
spectrum.  The  (1-1)  emissions  centered  on  771  nm  are  also  observed.  The  lifetime  of  the 
02(b)  state  is  about  1.1s,  allowing  the  nascent  distribution  to  be  collisionally  relaxed. 

Therefore,  it  is  possible  to  measure  the  gas  temperature  by  measuring  the  molecular 
rotational  temperature. 

Figure  4.2  illustrates  the  statistical  distribution  of  emission  intensity  versus  the 
rotational  energy  of  the  excited  state.  The  highly  linear  nature  in  Figure  4.2  exhibits  the 
statistical  nature  of  the  emissions,  arguing  that  the  reported  rotational  temperatures  will 
be  representative  of  the  gas  temperature.  As  discussed  in  Section  2.4.2,  Equation  2.4,  the 
slope  of  this  line  is  used  to  extract  the  rotational  temperature  of  the  02(b)  (0-0)  excited 

state.  The  intensities  are  obtained  from  a  fit  to  the  area  of  a  Gaussian  lineshape.  The 
reported  h  is  the  peak  intensity  of  the  Gaussian  fit,  while  the  width  of  all  the  peaks  are 
held  constant. 
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Figure  4.2  Natural  log  plot  of  emissions  from  02(b)  (0,0)  state  plotted  versus  energy  used  to  extract  gas 
temperature  from  emissions.  •  are  the  pP  branch  and  are  the  PQ  branch. 


The  error  bounds  represent  the  uncertainties  in  the  fit  areas,  as  generated  by  the 

86 

PeakFit  program,  which  is  typically  0.01%.  Several  gas  temperatures  have  been 
previously  reported  in  the  literature,  ’’’’’’  with  extracted  temperatures  of 
Tg  =  400  to  1,300  K  and  uncertainties  of  AT  =  20  to  61  K.  The  integrated  intensity  of  the 
02(b-X)  (1,1)  band  of  Figure  4.1,  relative  to  the  (0,0)  band,  is  9.2  +/-  0.3  %, 
corresponding  to  a  vibrational  temperature  of  896  ±  35  K.  Both  the  PP  and  PQ  branches 
yield  similar  temperatures;  Tg  =  720  D  19  and  77  =  72 1  □  20.  Gas  temperature  plays  an 
important  role  in  emission  intensity  and  is  discussed  further  in  Section  4.2. 
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4.2  Gas  Flow  Velocity  and  Residence  Time 

The  gas  temperatures,  Tg,  measured  at  the  midpoint  of  the  discharge  as  a  function 
of  oxygen  flow  rate  and  pressure  are  shown  in  Figure  4.3.  Note  the  oxygen  flow  rate 
may  be  changed  independent  of  pressure  by  throttling  the  pump  valve.  The 
measurements  reported  in  Figure  4.3  are  for  a  microwave  frequency  discharge  with  a 
100%  oxygen  gas  flow. 


Figure  4.3  Gas  temperature  as  a  function  of  flow  rate:  o  2  Torr,  •  3  Torr,  □  4  Torr,  5.2  Torr,  A 
6.4  Torr,  ▲  7.5  Torr,  V  8.6  Torr. 


The  gas  temperature,  Tg,  decreases  significantly  from  approximately  1,000  K  to 
600  K  as  the  flow  rate  increases.  However,  the  gas  temperature  appears  to  increase 
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slightly  as  a  function  of  pressure.  As  discussed  in  Section  2.5,  flow  rate  may  be 
converted  to  a  velocity,  given  the  flow  tube  geometry,  gas  temperature,  and  gas  pressure. 
If  it  is  assumed  that  the  velocity  of  the  gas  is  a  constant  throughout  its  travel  through  the 
discharge  region,  the  time  that  the  gas  spends  within  the  discharge  is: 


where  /  is  the  length  of  the  discharge  region,  vg  is  the  velocity  of  the  gas  flow,  and  t°r  is 
the  amount  of  time  the  gas  resides  within  the  discharge. 


Figure  4.4  Gas  temperature  as  a  function  of  residence  time  for  gas  pressures:  o  2  Torr,  •  3  Torr,  □ 
4  Torr,  5.2  Torr,  A  6.4  Torr,  ▲  7.5  Torr,  V  8.6  Torr. 
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The  validity  of  the  constant  velocity  assumption  is  discussed  further  in 
Section  4.3.1 .  The  length  of  the  discharge  changes  at  different  temperatures  and  mass 
flow  rates;  so,  the  length  is  measured  in  the  lab  and  varies  from  2  to  4  cm. 

Tg  as  a  function  of  residence  time  is  shown  in  Figure  4.4.  The  gas  temperature 
increases  dramatically  when  the  gas  enters  the  discharge,  changing  from  room 
temperature  to  1,200  K  in  one  millisecond.  It  then  comes  to  a  steady-state.  As  discussed 

in  Section  2.5,  it  is  appropriate  to  map  residence  time  to  a  pressure  nonnalized  time,  tr , 

which  collapses  the  pressure  dependence  of  the  data.  The  effect  of  this  mapping  is  shown 
in  Figure  4.5. 


Residence  Time,?/,  (ms) 


0.0  0.5  1.0  1.5  2.0  2.5 


Residence  Time,  C ,  (ms) 

Figure  4.5  The  gas  temperature  as  a  function  of  the  •  residence  time,  tr° ,  lower  axis  and  •  pressure 
normalized  time,  tr  ’,  upper  axis. 
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The  collapse  of  the  temperature  profiles  using  a  pressure  normalized  time  argues  for 
using  a  pressure-independent  temperature  profile,  as  derived  in  Section  2.5, 

Equation  2.18.  This  theory  is  used  throughout  the  remainder  of  this  work.  The  velocity 
calculations  and,  hence,  the  residence  times  have  been  calculated  assuming  the  gas 
temperature  is  constant.  Using  the  temperature  profile,  the  effects  upon  gas  velocity  and 
resonance  time  will  now  be  explored. 

4.3  Effects  of  Time-Dependent  Temperature 

The  pressure  gradient  in  the  flow  tube  is  not  large  enough  to  affect  the  number 
density  within  the  gas  flow,  <0.1  Torr  m  .  So,  the  pressure  is  assumed  to  be  constant 
along  the  length  of  the  flow  tube.  A  schematic  of  a  flow  tube  with  a  discharge  beginning 
at  position  xo  and  ending  at  position  xj  is  shown  below. 

Length  of  Discharge 

■ _ | _  Gas 

Flow 

◄ - 

P\vg\x\  Povgoxo 

Figure  4.6  Schematic  of  a  gas  flow  tube  with  positions  0  and  1  indicated. 

According  to  mass  conservation,  if  the  gas  density,  p,  changes,  then  the  velocity, 
vg,  must  change.  Therefore,  when  discussing  populations  as  a  function  of  time,  the 
temporal  axis  is  a  function  of  the  flow  velocity,  which  is  a  function  of  the  gas 
temperature.  The  relationship  between  length,  flow  velocity,  and  residence  time  is: 
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(4.2) 


~~  =  Vg(0  — ^  J Vg(f)dt  =  j*  dx  =  l 

k  xo 

where  /  is  the  length  of  the  discharge  region,  as  defined  in  Figure  3.6;  vg(t)  is  the  time- 
dependent  gas  flow  velocity;  xo  is  the  position  at  which  the  gas  enters  the  discharge;  x  /  is 
the  point  at  which  the  gas  exits  the  discharge;  to  is  the  time  at  which  the  gas  enters  the 
discharge;  and  ti  is  the  time  at  which  the  gas  exits  the  discharge.  The  velocity  of  the  gas 
is  a  function  of  the  time-dependent  pressure  based  on  the  ideal  gas  law  and  mass 
conservation.  Starting  with  the  ideal  gas  law: 

p  =  yk*T.  (4-3) 

where  P  is  pressure,  N  is  the  number  density,  V  is  the  volume,  kB  is  the  Boltzmann 
constant,  and  Tg  is  the  gas  temperature.  Mass  conservation  law  dictates: 

Avi  =  Pov  o  (4-4) 


where  p  is  the  density  of  the  gas.  Assuming  that  pressure  is  constant  along  the  axis  of  the 
flow,  manipulation  of  these  two  equations  shows  the  time-dependent  gas  velocity: 


v(0  =  v0 


T0 


(4.5) 


From  Section  2.5,  Equation  2.18,  the  time-dependent  temperature  is  modeled  as: 


Tj  e 

[P  +  T,a  e^) 


(4.6) 


So, 
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Solving  Equation  4.7  for  tr  yields: 

tr  =^-\n^j3  ~(T0a  + J3)  e  (4.8) 

Note  that  the  solution  for  the  time-dependent  temperature,  Equation  4.6,  is  in  terms  of  tr . 


Time,  using  v(t),  tr,  (ms) 

Figure  4.7  Comparison  of  residence  time,  defined  by  Equation  4.1,  assuming  a  constant  velocity 
throughout  the  discharge  region,  as  opposed  to  using  the  time-varying  velocity,  at  P  =  P0. 


Therefore,  pressure  nonnalization  is  built  into  the  solution  of  tr.  tr  is  the  amount  of  time 
the  gas  resides  within  the  discharge  region,  taking  into  consideration  the  time  variation  of 
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the  velocity  within  the  discharge.  The  effect  of  time  varying  temperature  on  the  overall 
residence  time  will  be  explored  in  the  following  discussion. 

4.3.1  Time  Varying  Velocity 

The  velocity  of  the  gas  flow  changes  as  the  temperature  of  the  gas  changes.  This 
change  in  velocity  will  affect  the  amount  of  time  the  gas  resides  within  the  discharge 
region.  However,  gas  temperature  reaches  a  steady-state  value  almost  immediately 
within  the  discharge  region,  as  shown  in  Figure  4.5.  Because  Tg  comes  to  steady-state  so 
quickly,  the  only  period  in  which  the  time  varying  velocity  calculation  differs  from  the 
constant  temperature  calculation  is  in  the  first  few  milliseconds  of  the  discharge. 
Therefore,  the  effect  of  including  the  time  varying  velocity  is  expected  to  be  small.  In 
Figure  4.7,  the  residence  time  is  reported  using  constant  and  time  varying  velocities.  The 
linear  relationship  between  these  two  values  shows  that  the  residence  time  of  the  gas  is 
virtually  unaffected  by  the  velocity  change  within  the  gas. 
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Figure  4.8  Two  plots  showing  the  difference  in  the  temperature  versus  time  plot  made  by  including  a 
time-dependent  velocity.  •  plot  using  f  and  •  plot  using  tr. 

Figure  4.8  shows  the  gas  temperature  as  a  function  of  both  forms  of  the  residence 
time,  tr  and/’ .  As  expected,  the  figure  shows  very  little  difference  between  the  two  plots. 
This  is  because  the  temperature  varying  velocity  does  not  have  a  significant  effect  upon 
the  residence  time  of  the  gas.  Even  though  a  change  in  velocity  has  a  negligible  effect  on 
the  residence  time  of  the  gas,  the  large  temperature  change  as  the  gas  enters  the  discharge 
requires  that  the  gas  temperature  be  well  defined. 
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4.4  Power  Loading  in  Gas  Flow 

When  analyzing  plasmas  sustained  in  flowing  gases,  the  temperature  of  the  gas 
plays  an  important  role.  The  temperature  of  the  gas  changes  many  important  parameters: 
the  kinetic  rates,  the  gas  flow  velocity,  and  the  number  density  of  the  emitters  in  the  gas 
flow.  The  intensity  of  the  emissions  from  the  02(  a)  state,  Ia>  is  dependent  on  the 

concentration  of  the  excited  state  [02(a]  ]  and,  therefore,  dependent  on  Tg,  as  governed 
by  the  ideal  gas  law.  Also,  the  velocity  of  the  gas  flow  determines  the  period  of  time  in 
which  the  gas  resides  within  the  discharge  region,  affecting  the  energy  deposited  into  the 
gas.  The  energy  deposited  into  the  gas  plays  a  large  role  in  determining  the  EEDF,  which 
drives  reaction  rates  of  all  reactions  involving  electrons.  So,  an  in-depth  understanding 
of  the  gas  temperature  in  the  flow  is  essential  to  interpretation  of  the  system  kinetics. 

As  discussed  in  Section  2.7,  the  functional  form  of  the  time-dependent  gas 
temperature,  Tg(t),  may  be  used  to  determine  the  value  of  the  heat  transfer  coefficient  of 
oxygen,  h  r,  and  the  fraction  of  energy  coupled  into  gas  heating,/  of  the  system. 

Fitting  this  functional  form  to  the  measured  temperature  data  yields  values  of: 

T0TSS=  6.4xl05  K2 

aTss  =  3900  5  (4.9) 

(Tss-To)  =  640K 

for  the  parameters  in  Equation  2.24.  This  system  of  equations  can  be  solved  giving 
values  of: 

T0  =542  ±18  K 

Tss=1182±10  K  (4.10) 

a  =  3.3±0.3  ^ 
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for  the  fit  parameters.  These  results  are  then  used  to  calculate  values  for  hr  and  f  as 
follows.  Using  Equation  2.14: 


PVCpd^  =  fpin-ShT(Tg-T») 


(4.11) 


a  may  be  written  as: 


pVCP 


The  values  for  this  calculation  are: 

Cp  =  1,005  J  kg1  K'1,  V=  1  n  r2  =  2.36  x  10'6, 
where  1  =  0.03  m,  r  =  0.  5  cm. 

The  density  of  oxygen  is: 

p=  (1.41x  103)  (7.5/760)  (3  00/1200)  kg  in3 


(4.12) 


(4.13) 


(4.14) 


at  7.5  Torr  and  1,200  K.  The  fraction  of  the  energy  coupled  into  the  overall  gas  heating 
is  found  from  the  detennination  of  the  value  for  //,  where: 


P  = 


IP. 

J  in 


pvcP 


(4.15) 


The  value  of  ft  is  determined  from  the  values  of  Tss,  To,  and  a  by  the  relationship: 

T0a  +  P 


T 


a 


(4.16) 


The  power  coupled  into  the  gas,  Pin,  is  85  W,  and,  from  Kays  and  Crawford,40  kj  = 
0.0717  (W  m"1  K"1).  The  theoretical  model  is  fit  to  the  measured  temperature  data,  and 
the  following  values  for  hr  and / are  extracted: 
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hT  =  29  ±  7  (W  •  m"2  •  K"1 ) 
/  =  17  ±2% 


(4.17) 


From  the  extracted  value  of  hT,  the  Nusselt  number,  Nu,  may  be  calculated: 


Nu  = 


hTD 

kT 


4.08  +  1.05 


(4.18) 


0.0  0.2  0.4  0.6  0.8  1.0  1.2  1.4 


Residence  Time,  K,  (ms) 

Figure  4.9  Temperature  as  a  function  of  time  in  the  discharge  is  governed  by  power  balance  of  heat 
coupled  into  the  discharge  and  heat  lost  to  the  wall.  Only  a  fraction  of  the  power  is  coupled  to  heat 
o  2  Torr,  •  3  Torr,  □  4  Torr,  5.2  Torr,  A  6.4  Torr,  ▲  7.5  Torr,  V  8.6  Torr,  —  Theoretical  fit. 


The  expected  Nu  is  4.364  for  cylindrical  geometries.40  This  value  is  well  within 
the  error  of  the  detennined  value.  The  steady  state  temperature  ,  Tss  =  1 182  ±  10,  is 
significantly  lower  than  anticipated  if  heat  transfer  to  the  wall  or  flow  of  excited  species 
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out  of  the  discharge  is  not  considered.  The  fit  result  /  =  0.17  ±  0.02  implies  most  of  the 
discharge  power  does  not  raise  the  gas  temperature.  The  analysis  presented  in  Appendix 
B  indicates  that  significant  power  is  transferred  to  the  reactor  walls  via  atomic  oxygen 
reactions.  The  reactor  wall  requires  active  cooling  with  compressed  air  to  prevent  glass 
softening.  Indeed,  thermal  control  and  surface  reactions  appear  important  in  most  EOIL 
demonstrations. 


Time,  t,  (ms) 

0  10  20  30  40  50  60  70 


Figure  4.10  •  Measured7,  —  Stafford7,  —  Nu  =  4.22  and  f  =  17%,  using  equation  2.18  to  describe  the 
heat  transfer  and  setting  the  input  power  outside  of  the  discharge  to  zero,  the  theoretical  result  is 
compared  to  the  measured  results  of  Carroll  and  simulated  results  of  Kushner. 
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Using  the  functional  form  of  gas  temperature  determined  in  Equation  2.18,  values 
for  hr  and / may  be  extracted  from  the  results  reported  by  Stafford  and  Kushner.7  The 
measured  data  in  Figure  4.10  are  the  results  of  an  experiment  performed  by  Carroll  et  al. 
at  UIUC  and  reported  by  Stafford  and  Kushner.7  The  flow  tube  has  a  4.83  cm  ID  with 
265  W  of  coupled  power,  via  an  RF  field,  into  a  flow  of  pure  O2  at  5  mmol/s  and  the 
pressure  being  held  at  2  Torr. 

The  fit  of  Equation  2.18  to  the  measurements  reported  by  Stafford  is  shown  in 
Figure  4.10.  From  position  z  =  0  to  20  cm,  the  gas  is  within  the  discharge  region;  at  20 
cm,  the  discharge  power  into  the  system  is  set  to  zero,  representing  the  end  of  the 
discharge  region.  The  differential  equation  that  describes  the  decay  of  temperature 
without  an  energy  source  is  a  simplification  of  Equation  2.14: 

dT 

pVCP-jt  =  -ShT(Tg-  T0)  (4.19) 


The  solution  to  this  differential  equation  is: 

f  _  hTs  .  T 

1  c 
0  1+  e 


Tg(tr)  =  T0 


CppV 


V 


J 


(4.20) 


Therefore,  outside  of  the  discharge,  the  temperature  decays  exponentially,  which  is 
clearly  the  case  in  Figure  4.10. 

Using  a  least-squares  fitting  method  inside  the  discharge  region,  the  values  of  hT 
and  / for  the  UIUC  results  are  fit  to  data  from  0  to  20  cm,  resulting  in  a  Nu  value  of  4.22, 
which  is  again  very  similar  to  the  theoretically  expected  values  of  4.364.  It  is  predicted 
that  only  17%  of  the  power  coupled  into  the  gas  results  in  heating.  The  fraction  of  power 
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coupled  into  gas  heating  in  p-wavc  discharges  and  RF  discharges  is  approximately  the 
same  for  these  studies.  This  is  an  unexpected  observation,  because  the  E-fields  are 
different,  and  the  geometry  of  the  two  experiments  is  different.  The  fraction  of  the 
energy  coupled  into  the  gas  flow  will  be  discussed  later  in  this  work. 

4.5  02(a)  Yield 

In  the  COIL,  population  inversion  of  the  iodine  atoms  is  achieved  using  the 
02(a)  state  as  an  energy  reservoir.  The  quantity  of  extractable  power  is  determined  by 

the  ratio  of  [02(a)  ]  and  [02(X)]  in  the  system.  Optical  monitoring  of  the  emissions 

from  the  excited  states  provides  a  convenient  way  to  determine  number  density  and  not 
perturb  the  performance  of  the  SOG.  The  relationship  between  number  density  of  the 
emitting  species  and  the  emission  intensity  is: 

Ia=CrAjNi  (4.21) 

where  Ia  is  the  emission  intensity;  Ay  is  the  Einstein  A  coefficient  for  the  transition;  N,  is 
the  number  density  of  emitting  species  in  the  discharge;  and  C,  is  a  factor,  including  the 
radiometry  of  the  detection  system. 

4.5.1  02(a )  Spectrum  and  Modeling 

13 

The  02(a  Ag— »X  Sg")  spectrum  centered  near  A  =  1.268  pm  has  nine  rotational 
branches  (akAJ  =  °P,  PP,  PQ,  QP,  qQ,  qR,  rR,  rQ,  sR),  as  discussed  in  Section  2.4.1. 37 
The  lifetime  of  the  02(a)  state  is  64  minutes,  making  emission  strengths  weak.  Thus,  a 
spectrometer  with  lesser  spectral  resolution  than  that  used  for  the  measurement  of  0-,  ( b ) 
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is  used  to  increase  the  throughput  and  signal  to  noise  ratio.  The  rotational  features  are 
not  resolved  as  they  were  for  the  02{b )  emissions.  Atypical  02(a)  spectrum  measured 
5  cm  downstream  of  the  discharge  is  shown  in  Figure  4.11. 

The  simulation  of  the  02  ( a )  emissions,  also  shown  in  Figure  4.1 1,  is  a  sum  of 

Lorentzian  peaks,  all  with  the  same  width  of  0.15  nm  and  relative  intensities  governed  by 
a  Boltzmann  temperature  distribution  at  600  K. 
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Figure  4.11  Typical  02(a)  spectra  resulting  from  the  discharge  and  a  simulation  of  the  spectra,  at  4 
Torr,  600  K,  and  500  SCCM,  5  cm  downstream  of  the  discharge,  o  measured  emissions,  —  simulated 
emissions. 
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The  temperature  is  a  fit  parameter  of  the  model  and  yields  values  very  similar  to 
the  results  determined  using  the  02 (b)  emissions  in  the  same  position  along  the  flow 
tube.  Peak  locations  and  Einstein  A  coefficients  are  extracted  from  the  HITRAN 
database  for  the  16Ot  (a'A-X3^)  (0-0)  and  160-180  (a*A-X3Z)  (0-0)  transitions.  The  1-1 
transition  is  centered  on  1,280  nm;  therefore,  the  contribution  from  this  emitting  state  is 
included  in  the  simulation.  The  Einstein  A  coefficient  of  the  1-1  transition  is  assumed  to 
be  the  same  value  as  the  0-0  transition  with  the  0-0  line  positions  scaled  by  the 
appropriate  energy  difference. 

Outside  of  the  discharge  region,  the  C>2(a— »X)  spectrum  is  well  isolated  from 
contaminating  emissions.  A  number  of  excited  O-atom  states  are  populated  within  the 
discharge,  resulting  in  several  atomic  oxygen  emission  lines  dominating  the  spectrum,  as 
shown  in  Figures  4.12  and  4.13.  In  Figure  4.12,  the  strong  3F  -  3D  (1,257  nm)  and  1,299 
nm  lines  bracket  the  02(a—>X)  emission.  In  Figure  4.13,  the  emissions  in  the  spectral 
range  of  1,266.8  to  1,277  nm  are  shown  to  be  largely  isolated  from  atomic  emissions. 
The  Q-branch  structure  is  still  clearly  evident.  The  X  =  1,266.8  to  1,277  nm  region  is 
integrated  to  provide  a  relative  measure  of  02(a)  concentration.  The  weaker  D-  P  and 
S-  P  lines  at  1,265  and  1,278  nm  are  outside  of  this  integration  band.  Spectral 
simulations  of  the  02{a—>X)  emission  have  been  developed  to  assess  the  dependence  of 
the  band  shape  on  rotational  temperature. 
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Figure  4.12  Emissions  from  inside  the  discharge  in  the  wavelength  region  of  02(a),  at  4  Torr,  and 
500  SCCM.  o  measured  data,  —  simulated  spectrum. 


The  gas  temperatures  inside  and  downstream  of  the  discharge  are  significantly 
different.  From  the  02(a)  simulation  in  Figures  4.1 1  and  4.13,  Tg  is  600  K  5  cm 

downstream  of  the  discharge  and  1,100  K  directly  inside  the  discharge.  The  temperature 
determined  inside  the  discharge  is  in  excellent  agreement  with  the  1 ,200  K  determined 
from  the  02(b )  emissions.  The  possible  sources  of  inaccurate  determination  of  the  Ia 

values  extracted  from  the  measured  data  are  the  wings  of  the  O-atom  emissions.  Figure 
4.14  will  take  a  closer  look  at  these  contributions. 
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Figure  4.13  Rescaled  view  of  emissions  from  inside  the  discharge  in  the  wavelength  region  of  02(a), 
at  4  Torr,  and  500  SCCM.  omeasured  data,  —  simulated  spectrum. 


Figure  4.14  shows  the  simulated  spectrum  with  the  contribution  from  the  02(a) 
(0-0)  transition  subtracted  from  the  total  spectrum.  The  remaining  signal  is  made  up  of 
the  02(a)  (1-1)  transition,  as  well  as  the  contribution  from  the  O-atom  emissions. 
Summing  the  remaining  signal  in  the  wavelength  range  of  1,266.8  to  1,277  nm  and 
dividing  by  the  sum  of  the  02(a)  emissions,  it  is  detennined  that  less  than  9%  of  the 
signal  is  made  up  of  emissions  not  directly  linked  to  the  02(a )  emissions. 
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Figure  4.14  Simulation  of  02(a)  with  contribution  from  the  (0-0)  transition  removed.  •  simulated 
spectrum  with  02(a)  (0-0)  contribution  removed,  02(a,  1-1),  —  O-atoms. 


Many  research  groups  use  02(a )  detection  techniques,  which  simply  include  a 

detector  and  a  band  pass  filter  or  very  low  resolution  spectrum  analyzers  with  a 
2  11  12 

bandwidth  of  100  nm.  ’  ’  “  These  techniques,  while  sufficient  if  the  detector  is  an 
appropriate  distance  downstream  of  the  discharge,  will  give  erroneous  results  if  the 
atomic  emissions  are  not  included  in  one’s  analysis. 
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4.5.2  02(a)  Yield  Calculation 

The  intensity  of  the  emissions  is  dependent  on  the  number  of  emitters  present  in 
the  field  of  observation,  as  well  as  the  Einstein  A  coefficient,  Ay,  of  the  excited  state. 

The  gas  is  heated  significantly  within  the  discharge,  and  the  number  density  is  changed 
accordingly.  Therefore,  the  change  of  the  emission  intensity  as  a  function  of  temperature 
must  be  accounted  for  if  the  kinetics  are  to  be  extracted  from  the  measurements.  Also, 
collision  frequencies  are  influenced  by  particle  densities. 


Mass  Flow  Rate,  h ,  (SCCM) 

Figure  4.15  Intensity  measurement  of  the  02(a)  emissions  versus  the  mass  flow  rate:  o  2  Torr,  •  3 
Torr,  □  4  Torr,  5.2  Torr,  A  6.4  Torr,  ▲  7.5  Torr,  V  8.6  Torr. 

In  order  to  determine  the  change  in  the  yield  of  02(a) ,  T  ,  as  a  function  of 
pressure  and  residence  time,  the  measured  quantities  are  the  intensity  of  the  02(a) 
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emissions,  Ia,  and  mass  flow  rate,  n  .  The  dependence  of  Ia  on  the  gas  temperature,  Tg, 
has  been  explained,  and  the  effects  have  been  removed  from  the  measured  data  in  order 
to  study  the  chemical  kinetics  of  the  system. 


Figure  4.16  Intensity  measurement  of  the  02(a )  emissions  versus  residence  time:  o  2  Torr,  •  3  Torr, 
□  4  Torr,  5.2  Torr,  A  6.4  Torr,  ▲  7.5  Torr,  V  8.6  Torr. 


The  data,  Ia  versus  the  mass  flow  rate,  h  ,  is  shown  in  Figure  4.15.  The  intensity 
of  the  02(a)  emissions,  Ia,  has  about  a  10%  error  bound  with  a  decreasing  trend  as  a 

function  of  the  gas  pressure  and  a  slight  decrease  at  higher  mass  flow  rates.  The  previous 
analysis  of  gas  temperature  is  used  to  map  from  the  mass  flow  rate,  h  ?  to  the  residence 
time,  tr,  resulting  in  Figure  4.16.  In  Figure  4.16,  the  high  mass  flow  rates  correspond  to 
the  largest  velocities,  which,  in  turn,  become  the  shortest  residence  times.  The  slight 


76 


change  in  intensity  at  the  high  mass  flow  rates  is  more  apparent,  showing  that,  before  the 
intensities  come  to  a  constant  value  within  the  first  millisecond  of  entering  the  discharge 
region,  the  concentration  of  02(a)  increases.  The  reason  for  measuring  the  02(a) 

emissions  is  to  study  the  changes  in  yield  as  a  function  of  pressure  and  velocity,  so  yield 
can  now  be  calculated. 

4.5.3  Scaling  Emission  Intensity  of  02(a)  to  Yield 

The  intent  of  this  study  is  to  explore  the  mechanism  responsible  for  limiting  the 
yield  of  02{a )  achievable  within  electric  discharges  in  oxygen  gas  flows  as  oxygen 
partial  pressure  increases.  A  mapping  of  measured  02{a )  emission  intensity,  Ia,  to  the 
absolute  concentration  of  02{a )  ,  [02(a)  ],  is  not  necessary  to  determine  plasma 
conditions  that  maximize  02(a)  production.  However,  the  majority  of  literature  that 
discusses  EOIL  systems  reports  results  in  terms  of  the  yield  of  02(a)  .  ’  ’  ’  In  order  to 
compare  the  current  results  to  those  published  in  the  literature,  a  scaling  term,  kr,  will  be 
detennined  that  will  enable  the  comparison  of  the  intensity  measurements  to  the  yields 
reported  in  the  literature. 

The  parameter  energy  per  particle,  tpart,  of  the  system  is  commonly  used  in 
literature.  ’  cpart  is  the  energy  coupled  into  the  gas  divided  by  the  number  of  particles 
that  flow  through  the  plasma  volume.  Note  that  the  number  of  particles  flowing  through 
the  discharge  volume  is  time-dependent. 

The  coupled  power,  length  of  the  discharge,  radius  of  the  flow  tube,  flow  velocity 
of  the  gas,  gas  stoichiometry,  and  gas  pressure  all  vary  between  experiments.  By  using 
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the  parameter,  energy  per  particle,  results  from  a  variety  of  different  discharge  conditions 
may  be  compared.  The  £,parl  for  the  measurements  made  in  this  study  may  be  written  as 
the  coupled  power  times  the  residence  time  of  the  gas  divided  by  the  number  density  of 
the  gas: 


^ part 


Pin  X  K 

[02(X)] 


(4.22) 


Figure  4.17  shows  the  02(a )  emissions,  measured  in  this  work,  for  a  plasma, 

sustained  by  an  electric  field  oscillating  at  the  same  frequency,  with  the  same  geometry, 
and  the  same  flow  conditions,  as  reported  by  Rakimova.  The  scaling  term,  kr,  is 
detennined  by  comparing  the  RF  results  to  those  by  Rakihmova  and  assuming  that 
concentrations  of  singlet  oxygen  will  be  the  same  under  the  same  discharge  conditions. 
Note  that  the  inner  diameter  of  the  flow  tube  used  in  this  study  is  1  cm,  which  is  0.05  cm 
smaller  than  the  flow  tube  used  in  Rakimova’s  work.  This  discrepancy  is  because  the 
wall  thickness  of  the  flow  tubes  in  the  two  studies  differs. 

The  02(a )  emissions  were  collected  for  three  pressures,  2,  4,  and  7  Torr.  When 
measuring  the  2  Torr  case,  the  flow  tube  overheated  when  the  gas  flow  within  the 
discharge  was  insufficient  to  cool  the  tube.  As  a  result,  analyzing  the  2  Torr  case  is 
difficult,  because  the  gas  does  not  reside  in  the  discharge  region  for  a  time  long  enough  to 
detennine  a  steady-state  population  for  this  pressure.  However,  a  scaling  tenn  may  be 
determined  when  using  the  4  and  7  Torr  measurements. 
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Figure  4.17  Emission  intensity  of  02(a)  as  a  function  of  energy  per  particle  in  the  current  RF 
discharge.  A  2  Torr,  o  4  Torr,  □  7  Torr. 


The  scaling  tenn,  kr,  was  calculated  for  the  4  Torr  case  and  was  found  to  be 
accurate  for  the  7  Torr  results  as  well,  showing  that  the  relationship  between  Ia  and 
[  0-,(a )  ]  is  linear  over  the  pressure  range  4-7  Torr.  Discussing  these  results  in  terms  of 


yield  requires  a  calibration  of  the  detection  system,  which  results  in  reporting 
concentrations  with  an  error  greater  than  20%.  The  majority  of  the  arguments  included  in 
this  work  will  depend  only  on  the  relative  changes  in  concentration. 

The  yield  of  02(a)  in  the  system  is  calculated  from  the  intensity  of  the  emissions, 
the  gas  temperature,  and  scaling  to  results  reported  by  Rakimova  : 


Yield,  T 


[02(X)] 


(4.23) 


79 


kr  =  Scaling  to  Rakimova 


where  Ia  is  the  measured  intensity  of  the  emissions;  \()2  (X )]  is  the  number  density  of 
the  ground  state  oxygen  within  the  flow  tube,  including  the  affect  of  gas  temperature. 
Note  that  the  temperature  dependence  of  the  emitter  is  already  built  into  the  yield 
calculation  via  the  temperature-dependent  ground  state  population  results. 

The  scaling  tenn  is  determined  by  comparing  the  RF  results  to  those  by 
Rakihmova.  Both  measurements  are  made  using  RF  discharges,  the  same  deposited 
energies  ranging  from  0.5  to  7.5  electron  volts  per  particle.  The  error  bars  are 
representative  of  the  error  associated  with  calibrated  optical  HPGe  detectors  used  in  the 
literature.  There  is  the  possibility  of  an  error  being  introduced  to  the  results  by  using  this 
scaling  system.  However,  discussions  in  chapter  5  are  based  on  the  relative  change  of 
yield  and  not  the  absolute  yield.  In  the  discussion  of  quenching  mechanism,  the  absolute 
yields  will  play  a  significant  role  only  if  they  are  orders  of  magnitude  different  from 
calculated  results.  It  is  unlikely  that  the  system  is  giving  results  below  5%  yields  or 
greater  than  15%  yields. 
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Figure  4.18  Yield  as  a  function  of  energy  per  particle  in  an  RF  discharge,  as  reported  by  Rakimova 
(solid  points)  and  this  work  (hollow  points),  o  4  Torr,  •  4  Torr  Rakimova,  □  7  Torr,  and  ■  7  Torr 
Rakimova.17 


Mapping  the  Ia  measurements  to  yields  allows  comparing  these  results  to  other 
published  works,  as  well  as  comparing  them  to  computer-simulated  concentrations. 
Besides  02(a) ,  electrons  and  atomic  oxygen  are  expected  to  have  significant 
concentrations  within  the  glow  discharge.  The  atomic  oxygen  density  is  detennined 
using  the  actinometry  technique,  which  is  discussed  in  Section  4.7.  The  electron  density 
of  the  glow  discharge  is  detennined  by  a  computer  model.  The  specifics  of  the  computer 
model  used  to  simulate  the  oxygen  discharges  in  this  study  are  discussed  in  Section  4.6. 
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4.6  Results  of  Simulations 


Current  systems  that  model  electric  discharges  in  oxygen  include  a  large  number 
of  pathways  for  creating  and  destroying  species  within  the  discharge.  These  numerical 
models  yield  results  which  fail  to  reproduce  pressure  trends  in  02{a )  concentrations. 
Detennining  means  to  correct  models  is  complicated  by  the  extensive  kinetic  packages 
implemented  in  these  models.  In  this  study,  we  seek  to  understand  why  the  02(a)  yield 
decreases  as  a  function  of  the  oxygen  pressure  within  the  gas  flow.  The  roles  of  electrons 
and  atomic  oxygen  in  the  destruction  of  02(a)  are  of  particular  interest  to  the  EOIL 

community.  In  order  to  understand  the  chemical  kinetics,  a  simulation  involving  only  the 
essential  chemical  kinetics  is  developed.  In  this  section,  the  streamlined  simulation  used 
in  this  study  is  described. 


4. 6. 1  Simplified  Analytic  Expression  of  Kinetics 

The  reactions  included  in  the  analytic  representation  of  the  system  are  shown  in 
Table  2.1.  The  electron  density  in  the  system,  at  steady-state,  is  determined  primarily  by 
direct  ionization  of  ground  state  molecular  oxygen  (reaction  3)  and  recombination  of 
electrons  with  positive  molecular  ions  (reaction  5): 


e  +  02  (  2f  ) — — — >  +  e  +  e 

e  +02  -,:'->02(X) 

e~  +02(X ) — >6T  +0 
e~  — — — >  Diffusion  from  Discharge 


(4.24) 


Electron  loss  is  dominated  by  attachment  thereby  eliminating  any  change  in  the  electron 
density  due  to  collisions  at  the  wall.  The  relative  rates  will  be  calculated  later. 
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O-atoms  are  generated  by  the  electron  impact  of  ground  state  oxygen  (reactions 
13-15)  and  destroyed  by  recombination  at  the  walls  (reaction  19): 


e+02(X)^ 

e+02(X)^ 

oCp)  +  o(3p)- 


>0{3P)  +  0(3P)  +  e~ 
■0(3P)  +  0CD)  +  e- 
»  02  (X) 


(4.25) 


The  concentration  of  02(a )  is  determined  by  direct  electron  impact  with  ground 

state  molecular  oxygen  (reaction  10)  and  the  inverse  reaction,  as  detennined  by  using 
detailed  balance  (reaction  17): 


e+02(X)- 
e~  +  02(a)~ 


>02(a)  +  e 
■02(X)  +  e 


(4.26) 


The  principal  species  involved  in  a  glow  discharge  are  ground  state  oxygen,  the 
02(a)  singlet  delta  molecular  oxygen  state,  the  0(3P)  atomic  state,  and  electrons.  The 

kinetic  energy  of  the  electrons  is  detennined  by  the  reduced  effective  electric  field.  The 
negative  ions  are  destroyed  primarily  by  ionic  recombination  (reaction  21): 


O  +  (X  — — 


■02(X) 


(4.27) 


This  streamlined  set  of  rates  is  used  in  conjunction  with  the  power  balance 
equation  (Equation  4.33)  to  describe  the  system,  assuming  that  steady-state  has  been 
reached: 


[°»]„  = 


*.-[«']  +r„  ' 

I-  -Iss 


(4.28) 


[«']  = 


Kn  [o2(X)\-katt[o2{X)\-kdifj 


(4.29) 
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(4.30) 

(4.31) 

(4.32) 


3  1 

in  units  of  (eV  m  s'  )  where  i  ranges  over  the  reactions  1-16  and  Aej  is  the  amount  of 
energy  required  to  excite  the  ith  excited  state.  Also,  As0  and  Asa  are  the  energy  required 
to  dissociate  02(X )  into  two  ground  state  O-atoms  and  the  excitation  of  02(X)  to  the 
02(a)  state  (6.0  eV  and  0.98  eV,  respectively).  The  steady-state  concentrations  of  the 
species  are  [02  (a)]  ss,  [02(X)]SS,  [0(3P)]SS,  [0]ss,  and  [e”]ss.  The  Ta  is  the  pseudo- 
first  order  decay  rate  used  to  compare  the  analytic  expression  of  the  02(a)  concentration 
to  the  measured  intensities.  Ta  is  used  in  Equation  4.28  instead  of  guessing  at  a  possible 
dominant  destruction  method  of  the  02(a)  state.  By  fitting  the  streamlined  equation  to 
the  measured  02(a)  data,  as  will  be  shown  in  section  4.8.1,  the  behavior  of  measured 
data  will  be  used  to  discuss  possible  destruction  methods.  The  rate  coefficients,  ka+  and 
ka-,  are  production  and  destruction  of  the  02(a)  state  by  electrons.  These  rate 
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coefficients  are  a  function  of  the  electron  temperature,  Te,  as  defined  in  Equation  2.8.  Te 
as  a  function  of  the  reduced  effective  E-field  is  detennined  by  solving  the  Boltzmann 
equation.  In  this  work,  BolSig+“  is  used  to  solve  the  Boltzmann  equation,  as  discussed 
previously.  All  of  the  Te  dependent  rate  coefficients  are  determined  using  the  BolSig+ 
package.  Electrons  are  produced  by  ionization  with  a  rate  coefficient  of  kj0n-  Loss 
mechanisms  for  the  electron  density  are  diffusion,  with  a  rate  coefficient  of  kditf, 
attachment,  with  a  rate  coefficient  of  katt,  and  electron  recombination  with  a  rate 
coefficient  of  krec •  The  plasma  is  assumed  to  be  quasi-neutral,  so  the  positive  ion  density 
must  be  equal  to  the  sum  of  the  negative  ion  density  and  electron  density.  O-atoms  are 
produced  by  dissociative  attachment  (reaction  1 8)  and  directly  by  electron  impact 
(reactions  14  and  15).  The  attachment  and  dissociation  rate  coefficients  are  summed  and 
represented  by  ko+.  At  steady-state  the  negative  ions  produced  by  dissociative 
attachment  are  considered  a  source  of  O-atoms.  The  destruction  of  negatively  charged 
atomic  oxygen  is  very  fast,  -  10'  cm'  s'  ,  and  results  in  negligibly  small  populations. 
This  assertion  will  be  discussed  later  in  this  section.  The  dominant  loss  mechanism  for 
O-atoms  is  recombination  at  the  walls  of  the  flow  tube,  kwau.  From  Gordiets  et  ah,  the 
radially  averaged  rates  for  wall  loss  can  be  expressed  as  : 

f  a2  2r  Y1  -i 

^Wall  ~  1  S 

{  D  Vo)  (4.34) 

where  A  is  the  characteristic  diffusion  length;  D  denotes  the  diffusion  coefficient;  r  is  the 
radius  of  the  cylindrical  flow  tube;  vth  is  the  thennal  velocity;  and  yo  is  the  probability  of 
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reaction  at  the  wall.  For  the  pressure  range  studied,  A 2/d  <^rlvthYo  .  This  assertion  is 
based  on  the  following  approximated  calculation: 

First,  from  Perram  et  al.,72’ 73  A  =  r/^ ,  where  r  =  0.005  m,  and  (f)  is  the  unitless 

2  3 

diffusion  parameter  for  a  cylinder  =  (2. 4048). ~  Therefore,  A  =  2.079  x  10”'  m. 

Also,  D  is  a  Fick  diffusion  coefficient,  which,  from  Yolles  et  al.,  is  equal  to 

2.7  x  10~5  (m2  s'1)  for  molecular  oxygen  and  02(b)  at  a  pressure  of  1  atm  and  300  K.  For 
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this  study  the  gas  in  the  discharge  region  is  at  1200  K  and  from  Wilke  and  Lee 
D  a  Tf  .  Therefore  at  4  Torr  and  1,200  K,  D  =  (760  /  4)  ( 1200/300)3/2  2.7  x  10  5 .  Using 
these  nominal  values,  the  diffusion-limited  portion  of  the  wall  rate  is 
A 2 ID  =  1.06  x  10  4  s.  So,  for  the  wall  reaction  tenn,  2r/vthy0  ;  the  velocity  of  the 

particles  in  the  gas  is  defined  as  vth  =  ^JSkBTg  jnm  =  1256  m/s  at  a  temperature  of 
1,200  K.  From  Gordiets  et  al.,18  yQ  ~  0.98  x  e~xl%^Twa"  is  the  probability  of  O-atoms 


reacting  at  the  wall  surface  and  is  equal  to  0.01877  at  450  K.  Therefore, 


2 r  _  2  x  0.005 

vthyQ  1256x0.01877 


4.24  x  lO^s. 


So,  the  statement  that  A2/d  <  2 r/vthy0  is  valid.  Therefore,  the  wall  rate  may  be  written 


as: 


k 


Wall 


VthYo 
2  XT 


(4.35) 


The  system  of  steady-state  relationships,  equations  4.28  to  4.33,  is  solved  using  a 
Mathematica  code.  The  inputs  to  the  Mathematica  code  are:  (1)  the  rates,  as  detennined 
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by  the  BolSig+  package;  (2)  the  gas  temperature;  and  (3)  the  velocity  of  the  gas  flow. 
The  measured  data  and  simulated  results  are  shown  in  Figure  4.19.  The  [02(a)  ]  comes 
to  steady-state  in  less  than  0.5  ms,  and  the  steady-state  value  is  inversely  proportional  to 
the  pressure  in  the  gas  flow. 


0.0  0.5  1.0  1.5  2.0  2.5  3.0 


Resdidence  Time,  K,  (ms) 

Figure  4.19  Analytic  model  of  the  yield  of  02(a)  gives  a  value  for  r„t  which,  on  average,  is  6,000  s'1:  o 
2  Torr,  •  3  Torr,  □  4  Torr,  A  5.2  Torr,  ▲  7.5  Torr,  and  —  theoretical  fits  from  Equation  4.28  at 
each  pressure. 


The  streamlined  set  of  equations,  4.28  through  4.33,  is  a  small  subset  of  the  rate 
packages  used  by  the  majority  of  the  groups  publishing  in  the  field.7' 11 ' 74  Note  that  Ionin 
asserts  that  a  unified  set  of  agreed  upon  rates  is  needed  in  current  rate  packages.40 
Therefore,  the  legitimacy  of  this  method  will  be  established  by  comparison  to  published 
results.  The  two  main  results  used  for  comparison  are  those  published  by  Stafford  and 
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Kusher4  and  Rakimova  et  al.17  The  Stafford7  paper  reports  the  simulated  concentrations 
of  a  capacitively  coupled  RF  discharge  sustained  in  a  3  Torr  pure  oxygen  flow.  The 
length  of  the  discharge  region  is  20  cm;  the  diameter  of  the  flow  tube  is  4.83  cm;  the  gas 

-5 

flow  velocity  is  1,000  cm/s;  and  the  deposited  power  is  0.5  W/cm  . 

The  results  of  the  streamlined  simulation  will  be  compared  to  the  simulated 
results  for  O-atoms,  02(a) ,  electron  density,  electron  average  energy,  and  the  E/N  of  the 
system.  The  results  reported  by  Rakimova  are  for  a  capacitively  coupled  RF  discharge 
in  pure  oxygen  of  pressures  ranging  from  2  to  20  Torr.  The  length  of  the  discharge  is 
approximately  3  cm,  the  radius  of  the  flow  tube  is  0.55  cm,  the  gas  flow  velocity  is 
630  cm/s,  and  the  coupled  energy  ranges  from  0  to  1,200  J/mmol  or  0  -  12  eV. 

In  Figure  4.20,  the  reduced  effective  electric  field,  E/N,  is  plotted  versus  the  flow 
velocity  of  a  gas  for  a  range  of  pressures.  The  E-field  is  oscillating  at  a  frequency  of 
2.45  GHz,  and  the  gas  flow  is  pure  oxygen.  For  each  pressure,  the  value  of  E/N  changes 
by  25%  over  the  range  of  flow  velocities,  while  the  average  value  of  E/N  for  the 
pressures  changes  by  a  factor  of  3. 

Te  drives  many  of  the  rates  included  in  the  model  and  is  detennined  by  the 
average  kinetic  energy,  £avg,  of  the  electrons.  £avg  as  a  function  of  E/N  is  an  output  of  the 

BolSig+  code.  If  these  values  match  those  of  other  reported  works,  then  the  rate  package 
included  in  BolSig+  is  adequate,  and  the  E/N  values  detennined  within  the  model  are 
representative  of  the  reported  work. 
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Figure  4.20  E/N  as  a  function  of  gas  flow  velocity  for  2.45  GHz  discharge  study  performed  in  this 
work:  •  2  Torr,  ■  3  Torr,  ▲  4  Torr ,  ▼  5.2  Torr,  ♦  7.5  Torr. 


Figure  4.21  shows  the  average  electron  energies,  £aVg,  calculated  for  three 

different  discharges.  The  simulation  presented  in  this  work  has  been  applied  to  each 
system.  Those  results  are  shown  in  black.  The  values  published  by  Stafford7  and 
Rakimova  are  the  open  points. 

The  comparison  to  the  Stafford  model  is  excellent,  with  the  values  differing  only 
by  0.05  eV.  The  simulated  Rakimova  energies  range  from  3.2  to  2.6  eV,  and  the  analytic 

IT 

simulation  ranges  from  3.36  to  3.0  eV.  Both  models  predict  & :avg  values,  which 

1  7 

decrease  as  a  function  of  pressure.  The  simulations  reported  by  Rakimova  and 
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Stafford7  use  Boltzmann  solvers  developed  by  those  research  groups,  as  opposed  to  the 
BolSig+  package  used  in  this  study. 


Pressure,  P,  (Torr) 

Figure  4.21  Simulated  results  for  average  electron  energy  by  o  Rakimova17;  •  this  work  simulation 
of  Rakimova  conditions;  A  Stafford7;  ▲  this  work  simulation  of  Stafford  conditions;  ■  this  work 
RF  conditions;  X  this  work  (x-wave. 


The  less  than  10%  error  in  the  absolute  values  of  the  electron  energies  is  most  likely  a 
consequence  of  the  different  Boltzmann  solvers  used,  as  well  as  the  different  rate 
packages.  The  calculated  reduced  electric  field,  E/N,  for  the  work  reported  by  Stafford  is 
shown  in  Figure  4. 22. 7 
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Figure  4.22  Comparison  of  steady-state  E/N  as  a  function  of  position.  (— )  this  work,  (— )  Stafford.7 

The  reduced  electric  field  as  a  function  of  position  is  shown  in  Figure  4.22  and  is 
compared  to  the  reported  results  of  Stafford.7  The  spike  in  Stafford’s  E/N  is  due  to  his 
simulation’s  prediction  that  the  discharge  needs  a  much  higher  electron  density  than  the 
results  detennined  by  the  steady-state  model.  As  a  result,  Stafford’s  model  predicts  a 
sharp  spike  in  the  E/N  at  the  start  of  the  discharge  region.  In  the  streamlined  model 
developed  for  this  work,  the  position,  x,  within  the  discharge  is  not  included.  Instead, 
change  in  distance  is  modeled  by  a  change  in  the  residence  time  or  gas  flow  velocity. 
Hence,  a  small  distance  is  a  short  residence  time  or  a  large  velocity.  From  Equation  4.33, 
at  large  velocities,  the  energy  lost  to  excited  species  is  significant.  Therefore,  this  term  in 
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the  power  balance  equation  dominants  the  values  of  E/N  and  therefore  the  electron 
temperature.  The  dominance  of  the  velocity  related  tenn  explains  the  change  in  the 
discharge  conditions  at  high  velocities,  short  dwell  times,  or  small  positions  within  the 
discharge.  By  solving  the  system  of  equations  at  steady  state  the  concentrations  vary 
smoothly  as  the  velocity  changes  and  therefore  there  is  no  spike  in  the  E/N  at  entrance 
into  the  discharge.  This  simulation  is  appropriate  for  this  work  because  the  discharge  is 
turned  on  minutes  before  the  measurements  are  made  allowing  the  discharge  to  light, 
thereby  eliminating  any  transient  energetics  which  are  included  within  time  dependent 
models. 

Figure  4.23  shows  the  comparison  of  the  simulated  concentrations  of  O-atoms  and 
02(a)  (using  the  Stafford-Kushner  model  versus  the  streamlined  model),  as  the  gas 

travels  through  the  discharge  region  (starting  at  0  cm  and  ending  at  20  cm).  The 
conditions  used  in  the  streamlined  model  are  the  same  as  those  simulated  by  Stafford  and 
Kushner.7  The  modeled  system  is  an  RF  discharge  operating  in  pure  oxygen.  The  flow 
tube  has  a  radius  of  4.83  cm,  the  discharge  length  is  20  cm,  the  gas  flow  velocity  is 
1,000  cm/s,  and  the  temperature  is  400  K,  in  the  discharge.  Note  that  the  radius  of  the 
flow  tube  in  the  Stafford  system  is  4.83  cm  as  opposed  to  the  0.5  cm  radii  used  in  this 
work  and  the  Rakimova  study.  The  large  radius  increases  the  value  of  the  volume  term 
in  Equation  4.33  thereby  decreasing  the  power  loading.  Reduced  power  loading  will 
decrease  the  rate  at  which  species  increase.  This  explains  why  the  excited  specie 
concentrations  continue  to  change  in  Figures  4.23  and  4.24  while  the  E/N  is  constant. 
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Figure  4.23  Comparison  of  O-atom  and  02(a)  concentrations  as  a  function  of  position.  (RF  discharge,  3  Torr, 
Pin  =  0.5  W/cm3,  /  =  20  cm,  r  =  4.83  cm.  vg  =  1000  cm/s.)  (— )  this  work  O-atoms,  (— )  Stafford  O-atoms,7  (-  -) 
this  work  02(a),  (-  -)  Stafford  02(a).7 


The  functional  form  of  the  simulations  is  very  similar  with  a  small  difference  in  the  rate 
of  curvature.  The  differences  in  curvature  are  related  to  the  destruction  and  creation 
pathways  included  in  the  Stafford-Kushner  model  which  are  not  included  in  the 
streamlined  model.  The  concentrations  at  the  exit  of  the  discharge  are  within  a  factor  of 
2  making  the  overall  agreement  of  these  models  acceptable.  The  electron  densities 
produced  by  the  two  models  are  also  compared. 
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Figure  4.24  Comparison  of  electron  density  as  a  function  of  position.  (RF  discharge,  3  Torr,  Pin  =  0.5  W/cnr\ 
/  =  20  cm.  r  =  4.83  cm,  vg  =  1000  cm/s.)  ( — )  this  work,  ( — )  Stafford.7 


In  Figure  4.24,  the  Stafford  results  are  reported  in  gray,  while  the  model  presented  in  this 
study  is  shown  in  black.  The  overall  shapes  of  the  two  curves  are  once  again  very 
similar;  the  only  difference  between  the  two  models  is  the  resultant  absolute 
concentrations,  which  differ  by  only  a  factor  of  3.5  or  less. 

Within  this  work  the  steady  state  value  of  the  electron  density  is  determined  using 
Equation  4.29.  Using  this  equation,  electrons  are  removed  from  discharge  region  by 
dissociative  attachment,  recombination,  and  diffusion.  If  diffusion  is  the  dominant  loss 
mechanism  of  electrons  within  the  plasma,  then  the  EEDF  could  be  affected  by  electron 
collisions  at  the  wall.  Using  Table  2.1at  4  Torr,  Tg  =  1,185  K,  Te  =  3.2  eV,  and  a  tube 
radius  of  0.5  cm,  nominal  values  for  the  electron  loss  reaction  rates  in  Equation  4.29  are: 
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K,,ch  =1-9x10  11  ans  1 

krec  =  1.6  x  10  s  cm3  s~x 
and  the  loss  rates  are: 

Kttch  x  [02(X)]  =  7.6  x  105  51 

k™  x[e"]  =  5°5“1 

k-diff  =  2-7  x  104  5-1 

where  a  nominal  ground  state  concentration  of  4  x  1016  cm  3  is  used.  The  electron 
concentration  of  5  x  I  ()  '  cm  is  used  from  Rakimova,17  shown  in  Figure  4.25.  Diffusive 
losses  for  electrons  is  kdiff  =  D/A  (1+40  Te  (293/Tg)).  The  rate  of  dissociative  attachment 
is  an  order  of  magnitude  larger  than  the  diffusive  losses  and  therefore  the  affects  of 
electrons  colliding  at  the  wall  will  not  significantly  affect  the  EEDF.  In  Figure  4.25, 
simulated  electron  densities  for  a  number  of  RF  discharges  are  reported  as  a  function  of 
pressure. 

The  electron  number  densities  simulated  by  Rakimova  decrease  as  pressure 
decreases,  while  those  calculated  in  this  work  increase  as  pressure  is  lessened.  There  is 
some  disagreement  in  the  literature  as  to  the  behavior  of  electron  densities  as  a  function 
of  pressure.  ’  ’  ’  Many  studies  that  report  an  increase  in  [e  ]  with  increasing 
pressure  are  perfonned  at  a  much  lower  pressure,  less  then  500  mTorr,  and  suggest  that 
the  direct  relationship  between  electron  density  and  pressure  stops  or  is  reversed  at  higher 
pressures.76  The  rationale  for  a  decrease  in  [  e~  ]  as  a  function  of  pressure  is  straight¬ 
forward  and  addressed  below. 
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Figure  4.25  Simulated  results  for  [e~]  by  o  Rakimova,17*  this  work  (Rakimova  conditions), 
work  (Stafford  conditions). 
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Stafford,7  A  this 


4. 6.2  Decrease  of  Electron  Number  Density  as  a  Function  of  Pressure 
The  electron  number  density,  as  presented  in  this  work  and  Rakimova’ s  work,  is 
determined  by  solving  a  number  of  coupled  differential  equations.  Neither  study  includes 
a  measurement  of  the  electron  density.  Both  models  assume  that  the  plasma  fills  the 
cross-sectional  area  of  the  flow  tube  in  the  discharge  region.  This  assumption  ignores  the 
tendency  of  the  plasma  to  concentrate  into  a  toroidal  shape,  at  a  higher  pressure.  (This 
phenomenon  is  caused  by  the  reduction  in  the  mean  free  path  of  the  electrons.)  As  such, 
the  electron  number  densities  detennined  by  both  models  are  flawed.  It  is  the  intent  of 
this  section  to  present  material  that  suggests  the  feasibility  of  an  inverse  relationship 
between  electron  density  and  gas  pressure  in  the  flow. 
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From  Figure  4.26,  it  is  apparent  that  the  ionization  rate  is  strongly  dependent  on  the 
average  energy  of  the  electrons.  From  Figure  4.21,  both  models  simulating  the  Rakimova 

geometry  show  a  10%  decrease  in  £avg  with  pressure.  From  Figure  4.24,  the  ionization 

rate  increases  from  a  value  of  2  x  1010 cm3,?-1  at  3  eV  to  1 1  x  10 l0crn  s~l  at  3.3  eV. 
Dissociative  attachment,  the  dominant  loss  mechanism  for  electrons,  does  not  change 
dramatically  as  will  be  discussed  later  in  this  section.  Therefore,  it  is  expected  that  the 
electron  density  will  decrease  as  the  pressure  of  the  gas  increases.  However,  Rakimova’s 
work  shows  an  increase  in  electron  density  as  the  gas  pressure  increases.  The  source  of 


Average  Electron  Energy,  sgvg,  (eV) 

Figure  4.26  Ionization  rate,  kion,  as  a  function  of  the  electron  temperature,  as  determined  using  the 
BolSig+  model,  for  a  p-wave  discharge  at  4  Torr. 
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1  7  • 

electrons  at  a  higher  pressure,  as  reported  by  Rakimova,  is  a  collisional  detachment  of 

electrons  from  negative  ions.  The  O  population  is  produced  by  attachment  and 

quenched  by  collisions  with  O-atoms,  02(a) ,  and  positive  ions.  Therefore,  an  increase 

in  the  number  of  collision  partners  logically  results  in  an  increase  in  detachment. 

The  mechanism  for  detachment  is  shown  below: 

O  +  Ad  — ^  O  ~\~  Ad  )■  c  , . 

(4.36) 

where  M  represents  any  collision  partner  in  the  gas  flow.  From  this  expression  it  may  be 
surmised  that  as  the  gas  pressure  in  the  flow  is  increased,  and  subsequently  the  number 
density  of  collision  partners  increase,  that  the  electron  concentration  in  the  flow  will 
increase.  However,  to  increase  [e\,  [O  ]  must  be  the  same  order  of  magnitude  as  [e\  at 
lower  pressures  that  will  serve  as  a  reservoir  of  electrons  that  may  be  depleted  as  pressure 
is  increased.  Examination  of  the  second  order  effects  of  the  negative  ion  density  on  the 
electron  density  follows. 

In  Figure  4.27,  the  rate  of  production  of  O  within  an  RF  discharge  is  shown  for 
a  range  of  average  electron  energies.  The  typical  Savg  for  an  RF  discharge  at  pressures 
over  1  Torr  is  approximately  3.2  eV.  From  the  chart  below,  the  value  is  approximately 
kattch  =  2x10'  cm  s'  ,  assuming  an  upper  limit  for  E/N  as  60  Td  or  Savg  =  4  eV.  In 

order  to  examine  the  upper  limit  of  the  negative  ion  concentration,  the  higher  E/N  will  be 
used. 
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Average  Electron  Energy,  e  (eV) 

Figure  4.27  Attachment  rate  in  cm3  s'1  plotted  versus  average  electron  energy  (eV).  •  p-wave.  o  RF. 


The  solution  to  the  concentration  of  O  within  the  discharge,  assuming  quasi-neutral 
discharge  at  steady-state,  having  a  gas  temperature  of  1 ,000  K,  is  governed  by  the 
reactions  below: 


^[°^,+tJ°2<a>L+W0L 


(4.37) 


3xl0-n[l010]  [lO16] 

1°  L  ~2.98xl0~8[l010]  +1.2xl0-9[l015]  +  8xl0~10[l015] 


(4.38) 


4  x  1015 
1.2x1  (f 


«  3  x  109 


(4.39) 
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The  approximate  value  for  the  [  O  ]  ss  value  is  roughly  an  order  of  magnitude 
smaller  than  the  electron  density.  This  approximate  value  was  calculated  assuming  a 
large  value  for  the  attachment  rate  and  that  both  the  02(a)  and  O-atom  population  are 

1%.  With  a  large  Savg,  production  is  at  a  maximum,  and,  with  small  concentrations  of  O 

and  02(a) ,  quenching  of  O  is  small  regarding  actual  plasma  conditions.  Even  after 

assuming  artificially  large  production  terms  and  small  quenching  terms,  the  O 
population  could  only  reach  10%  of  the  electron  density.  If  all  of  the  O  present  within 
the  discharge  is  collisionally  relaxed,  producing  free  electrons,  then  the  total  increase  in 

the  electron  density  would  be  10%.  As  the  pressure  increases  from  2  to  20  Torr,  the  8avg 

decreases  from  3.2  to  2.6  eV,  causing  a  reduction  in  the  ionization  rate  of  greater  than 
50%.  Therefore,  including  the  50%  decrease  in  electron  density  caused  by  reduction  in 
the  rate  of  ionization  and  the  10%  increase  due  to  increased  detachment,  the  overall 
change  to  the  electron  density  as  pressure  increases  should  be  a  reduction  of  40%. 

4. 7  O-atoms  within  the  Discharge 

One  of  the  major  differences  in  the  effluence  of  a  chemically  driven  SOG  and  a 
SOG  driven  by  an  electric  discharge  is  the  production  of  atomic  oxygen.  As  shown  in 
Figure  2.5,  the  EEDF  within  a  discharge  has  a  high-energy  tail,  hence,  there  are  a  finite 
number  of  electrons  with  enough  energy  to  dissociate  an  oxygen  molecule.  Also,  if 
significant  percentages  of  metastable  states,  such  as  02  ( a ) ,  are  maintained  in  the 
discharge,  then  these  excited  states  will  also  be  dissociated  by  electrons  with  energies  less 
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than  6.0  eV,~  which  is  the  energy  necessary  to  dissociate  an  oxygen  molecule.  O-atoms 
are  also  created  by  dissociative  attachment,  reaction  2,  which  has  a  threshold  energy  of 
approximately  3  eV.  O-atoms  introduce  reaction  pathways  that  do  not  exist  in  the 
chemistry  based  oxygen-iodine  laser  system  and  may  reduce  the  02(a)  population  and  I* 

concentration.  If  O-atoms  reach  a  significant  population,  they  could  influence  the  EEDF 
and  change  the  reaction  rates  determined  by  using  the  BolSig+  code  in  a  pure  oxygen  gas 
flow.  In  order  to  verify  that  the  O-atom  concentration  is  small  enough  not  to  affect  the 
EEDF,  O-atom  concentrations  have  been  measured.  Note  that  the  effect  of  O-atoms  on 
the  EEDF  has  been  added  to  the  BolSig+  model,  assuming  that  the  O-atom  concentration 
makes  up  less  than  10%  of  the  discharge.  (The  assumption  that  O-atom  concentration  is 
less  than  10  %  will  be  validated  by  actinometry  measurements  in  section  4.7. 1 .)  This  has 
been  accomplished  by  adding  O-atom  cross-sections,  from  Laher  et  ah,  to  the  BolSig+ 
library,  scaling  the  cross-sections  to  reflect  that  [O]  =  10%  of  [02(X)]  and  recalculating 
the  EEDFs.  The  result  is  that  there  is  less  than  a  3%  change  in  the  production  rate  of 
Ch(a)  .  Because  the  error  in  the  determination  of  yield  is  approximately  20%,  it  will  be 

assumed  that  the  change  in  the  rate  is  negligible  if  the  population  of  O-atoms  is  less  than 

10%. 

Measurement  of  the  O-atom  population  is  accomplished  by  using  the  established 
method  of  actinometry,  which  was  discussed  in  Section  2.8.  Measured  populations  for 
this  study’s  microwave  discharge  will  be  discussed  in  the  contexts  of  existing  literature, 
and  any  ramifications  will  be  examined. 
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4.7.1  Actinometry 

From  Pagnon,  the  [O]  may  be  determined  from  the  ratio  of  1 0844nnM Ar750nm  by 
the  relationship  shown  in  Section  2.8,  Equation  2.30.  The  emissions  of  IQ  and  I  a,-  that 
result  from  a  microwave  discharge  operating  at  85  W,  800  SCCM,  and  10%  argon  in  the 
flow  are  measured  and  shown  in  Figures  4.28  and  4.29.  The  emission  from  the  argon 
centered  at  750.4  mn  is  the  result  of  the  transition,  2 Pm  4p  (0.99eV)  -  2 Pm  4s  (0.87eV). 

The  peak  centered  at  75 1 .5  mn  is  a  result  of  the  transition  from  "Pm  4p  (0.97eV)  . 

2 

P3/2  4s  (0.85eV).  The  measurements  in  this  work  show  that  O-atoms  make  up  less  than 
5%  of  the  gas  flow  in  an  apparatus  with  a  radius  of  0.5  cm  for  a  p-wavc  discharge, 

Figure  4.30.  The  Rakhimova  measure  shows  that  O-atoms  make  up  about  5%  of  the  gas 
flow  in  a  RF  discharge.  The  argon  emission  line  is  at  750.4  mn,  and  the  atomic  oxygen 
line  at  844.6  mn  are  the  primary  lines  used  for  the  actinometry  detection  method.  Using 
the  actinometry  technique,  the  degree  of  dissociation,  as  measured  by  Ivonov,77  is  2%. 
The  streamlined  model  predicts  a  value  of  5%  under  the  Ivonov  conditions.  Other 
reported  yields  of  O-atoms  are:  Stafford  20%, 7  Sharma  0.2%, 78  Rakihmova  5%, 17  and  5 
±  0.5%  in  this  study. 
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Figure  4.30  Measured  percentage  of  the  gas  flow,  which  is  made  up  of  O-atoms  within  a  2.45  GHz 
discharge  at  three  pressures,  2,  4,  and  7  Torr. 

The  actinometry  measurement  agrees  with  the  simulated  results,  as  well  as  those 
reported  in  the  literature  for  similar  discharge  conditions.  As  previously  mentioned,  it 
will  be  assumed  that  a  yield  of  O-atoms,  less  than  10%,  will  have  no  effect  upon  the 
electron  energy  distribution.  Therefore,  the  BolSig+  results  are  accurate  for  the 
continued  simulation  of  this  work. 

4.8  Results  and  Discussion  of  02(a)  Data 

The  results  found  from  the  literature  (shown  in  Figure  2.14)  and  the 
measurements  in  this  work  (shown  in  Figure  4.18)  both  exhibit  a  decrease  in  the  yield  of 
02(a)  as  pressure  increases.  In  order  to  describe  the  behavior  of  the  reduction  of  the 


x 

CM 

O 


03 


104 


02{a )  concentration,  a  pseudo-first  order  decay  rate,  Ta,  has  been  introduced  to  the 
steady-state  solution  of  the  02(a)  excited  state.  The  behavior  of  Ta  will  be  used  to  asses 
a  destruction  mechanism  for  the  02(a )  state. 


4. 8. 1  Inclusion  of  Pseudo-First  Order  Decay  Rate 

The  concentration  of  02(a )  as  a  function  of  gas  pressure  and  flow  velocity  has 

been  measured  by  monitoring  the  optical  emissions  centered  on  1260  mn.  The  yields 
calculated  from  these  measurements  show  an  inverse  relationship  to  the  oxygen  pressure 
within  the  flow.  The  steady  state  equation  used  to  describe  the  Of  a)  concentration 
within  the  discharge  includes  a  pseudo-first  order  decay  rate.  Obtaining  values  of  V a 
using  the  steady  state  equation  in  combination  with  the  yields  of  Ofa )  ,  will  permit  an 
assessment  of  the  destruction  mechanisms  of  02(a )  within  the  discharge.  Calculation  of 
the  first-order  decay  rate  for  a  range  of  flow  velocities  and  gas  pressures  follows. 

In  Equation  4.28,  the  steady-state  solution  of  02{a )  includes  a  pseudo-first  order 
decay  rate  denoted  by  the  symbol  V a .  The  values  of  Ta  are  detennined  by  solving  the 
steady-state  equation  for  fn : 


r 


a 


ka+  e  ,  r  -i 

- *-[e  1 


(4.40) 


Note  that  the  fraction  [02(X)]ss/  [02(a)]  s  is  present  in  Equation  4.40.  This  value  is  the 
inverse  of  the  yield  of  02(a) ,  Ya ,  which  is  detennined  from  the  measured  intensity  of  the 
02(a )  emissions,  Ia: 
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Using  equations  4.28  through  4.33,  there  is  a  system  of  six  equations  with  seven 
unknown  values,  including  V a  as  an  unknown.  By  including  Equation  4.41  as  the 

solution  for  the  steady-state  value  of  02(a) ,  the  system  of  equations  may  be  used  to 
obtain  values  of  Ta . 


Figure  4.31  ra  as  calculated  for  each  datum  measured  with  a  p-wave  discharge  as  a  function  of  flow  rate  for 
the  pressures:  o  2  Torr,  •  3  Torr.  □  4  Torr.  A  5.2  Torr,  ▲  7.5  Torr. 


Figure  4.3 1  shows  the  values  of  U(  determined  for  pressures  ranging  from  2  to 
7.5  Torr  and  mass  flow  rates  varying  from  10  to  1,800  SCCM.  By  averaging  the  value  of 
over  the  range  of  flow  velocities,  an  average  value  of  is  calculated  for  each 
pressure.  These  average  values  vary  with  pressure,  as  shown  in  Figure  4.32.  The 
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average  values  of  gamma,  ovgra ,  are  near  6,000  s'1  for  p-wave  discharges,  with  the 
exception  of  the  values  for  2  Torr  and  7.5  Torr. 


1  2  3  4  5  6  7  8 


Pressure,  P,  (Torr) 

Figure  4.32  Average  value  of  r„  as  a  function  of  pressure.  •  p-vvave  conditions  for  2-7.5  Torr,  A  RF  conditions 
(this  work),  ■  RF  conditions  (Rakimova),17  and  ▼  RF  conditions  (Stafford).7 


Figure  4.32  shows  that  the  decay  rate  for  RF  discharges  is  an  order  of  magnitude 
less  than  those  values  calculated  for  microwave  discharges.  This  relationship  is  telling 
with  regard  to  any  suggested  mechanism.  The  behavior  of  as  a  function  of  flow 
velocity,  E-field  oscillating  frequency,  and  gas  pressure  will  be  used  to  examine 
candidate  kinetic  mechanisms  responsible  for  limiting  the  yield  of  02(a)  . 

By  changing  the  frequency  of  the  oscillating  field  from  a  microwave  frequency  to 
RF,  the  average  electron  energy,  savg,  changes  by  0.2  eV  or  less,  as  shown  in  Figure  4.21. 
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Average  Fit  Parameter  p-wave,  avgra,  (s'1) 


The  small  change  in  savg  has  a  relatively  minor  effect  upon  most  of  the  reaction  rate 
coefficients  in  this  energy  region,  as  shown  in  shaded  areas  of  Figures  2.6  and  2.7.  A 
closer  look  at  the  rates  important  to  production  and  destruction  of  electrons  is  shown  in 
Figure  4.33. 


Figure  4.33  Reaction  rate  coefficients  at  3  Torr  for  o  dissociative  attachment  p-wave,  •  dissociative 
attachment  RF,  A  dissociation  p-wave,  ▲  dissociation  RF,  □  ionization  p-wave,  ■  Ionization  RF 
(from  BolSig+). 


However,  the  ionization  rate,  kion,  varies  strongly  as  a  function  of  £avg  in  this 
energy  region.  As  shown  in  Figure  4.33,  this  rate  increases  by  an  order  of  magnitude  or 
greater  when  changing  from  RF  to  p-wave  frequencies.  Because  there  is  an  order  of 
magnitude  difference  in  Tfl  between  the  RF  and  microwave  discharges,  the  mechanism 
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responsible  for  quenching  02(a )  is  probably  related  to  the  ionization  rate  and,  most 
likely,  the  electron  density.  In  Figure  4.34,  the  kion  at  4  Torr  is  shown  for  both  RF  and 
microwave  frequency  electric  field,  as  detennined  by  using  the  BolSig+  code. 


Average  Electron  Energy,  savg ,  (eV) 

Figure  4.34  Ionization  rate  versus  average  electron  energy  for  —  |j.-wave  frequency  and  —  RF 
discharges  at  4  Torr,  Tg  =  1,100  K  (from  BolSig+). 


The  dimensions  of  the  discharge  used  to  produce  Figure  4.34  are  a  flow  tube 
diameter  of  1  cm,  the  length  of  the  discharge  region  of  3.2  cm,  and  the  average  electron 
energy  ranges  from  2.5  to  3.7  eV,  which  is  the  range  of  savg  studied  in  this  work. 

Figure  4.35  shows  the  near  linear  dependence  of  the  pseudo-first  order  decay  rate 
divided  by  the  steady-state  concentration  of  oxygen  in  the  ground  state,  ra/[02(X)]  , 
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on  the  lei.  This  argues  for  a  quenching  mechanism  that  is  dependent  on  I  e 

L  J  ss  L  Ass 

multiplied  by  [02(X)]  . 


Figure  4.35  Relation  of  r„  /[02(X) ]ss  versus  [e]  ss.  o  2  Torr,  □  3  Torr,  A  4  Torr,  5  Torr.  V  7  Torr,  • 
current  work  RF  2,  4,  and  7  Torr,  ♦  Stafford  RF7,  Rakimova  RF.  17 

The  electron  densities  are  determined  by  using  the  streamlined  computer 
simulation,  using  equations  4.28  through  4.33.  From  the  near  linear  relationship  shown 
in  Figure  4.35,  the  dominant  mechanism  for  limiting  the  must  be  of  the  form: 

r.=I.[>l,[0^>L  <4  42> 

where  is  the  rate  coefficient  required  for  this  reaction  mechanism.  As  governed  by 
the  power  balance  equation.  Equation  4.33,  the  [e ~]  and  [02(  A)]  are  inversely 
related,  making  their  product  invariant  with  residence  time. 
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Flow  Velocity,  vg,  (cm/s) 

Figure  4.36  The  ground  state  and  electron  densities  are  shown  in  the  plot  above,  including  the 
product  of  the  two  concentrations  as  a  function  of  flow  velocity  at  4  Torr.  •  \02(X)\,V  \e ],  ▲ 
\02{X) ]  •  \e\. 

The  relationship  between  the  electron  density  and  ground  state  density  is  shown  in 
Figure  4.36.  Two  key  characteristics  of  the  dominant  quenching  mechanism  of  02(a ) 

have  been  identified.  First,  ra/[02(3f  )]  ss  is  linear  with  electron  concentration,  and 
second,  Ta  in  an  RF  field  is  significantly  smaller  than  in  the  presence  of  a  microwave 
frequency  field.  These  characteristics  are  used  to  explore  possible  mechanisms 
responsible  for  destroying  02(a)  . 
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4. 8. 2  Brief  Synopsis  of  Results 


A  systematic  study  of  microwave  discharges  has  been  perfonned  through  the 
pressure  range  of  2  to  7.5  Torr  in  pure  oxygen  flows.  Both  RF  and  p-wave  frequencies 
have  been  studied.  Optical  spectroscopy  has  been  used  to  monitor  the  change  in  the 
[02(a)] ,  [02(b)] ,  and  \0~\  concentrations  as  a  function  of  pressure  and  gas  flow 

velocity.  A  technique  involving  spectroscopy,  gas  flow  velocity,  and  careful  analyses  of 
the  gas  temperature  has  been  developed  and  used  to  observe  the  excited  species  directly 
within  the  discharge  region.  In  the  literature,  experimental  measurements  of  excited 
effluent  from  discharges  are  typically  performed  many  centimeters  downstream  of  the 
discharge.  When  optical  diagnostics  are  used  immediately  outside  of  the  discharge  or 
within  the  discharge,  the  intensity  of  the  emissions  is  strongly  affected  by  the  gas 
temperature.  Because  of  this  dependence,  an  elegant  method  has  been  developed  to 
precisely  detennine  the  temperature  of  the  gas  in  the  discharge  region. 

The  gas  temperature  measurements  reported  in  this  work  have  been  extracted 
from  the  rotational  energy  distribution  of  the  excited  \02  (b)]  state.  The  resolution  of  the 
[0,(b)]  spectral  data  reported  in  this  study  has  a  resolution  of  0.01  mn,  while  the  most 

precise  resolutions  reported  in  the  literature  is  0.04  mn.  Higher  resolution  spectra  allow 
the  determination  of  the  gas  temperature  to  a  higher  degree  of  precision.  Typical  errors 
in  reported  gas  temperatures  for  this  study  are  5%  or  less,  whereas  errors  in  gas 
temperature  reported  in  the  literature  are  typically  10%  or  greater.  Because  the  gas 
temperature  influences  the  intensity  of  the  emissions  from  all  of  the  measured  species, 
higher  precision  temperature  determination  results  in  higher  precision  density 
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measurements.  The  high-resolution  gas  temperature  measurements  perfonned  in  this 
study  improve  the  quality  of  the  reported  data  by  approximately  10%  throughout  this 
work. 

The  primary  purpose  of  the  gas  temperature  measurement  is  to  obtain  an  accurate 
measurement  of  [02(o)]  using  optical  intensities.  An  understanding  of  the  gas 

temperature  also  pennits  a  brief  commentary  upon  thermal  transport  within  the  discharge 
and  the  affects  it  has  upon  EOIL  systems.  The  steady  state  temperature,  Tss  =  1 182  ±  10, 
is  significantly  lower  than  anticipated  if  heat  transfer  to  the  wall  or  flow  of  excited 
species  out  of  the  discharge  is  not  considered.  The  fit  result  /  =  0.17  ±  0.02  implies 
most  of  the  discharge  power  does  not  raise  the  gas  temperature.  The  analysis  presented 
in  Appendix  B  indicates  that  significant  power  is  transferred  to  the  reactor  walls  via 
atomic  oxygen  reactions.  The  reactor  wall  requires  active  cooling  with  compressed  air  to 
prevent  glass  softening.  Indeed,  thermal  control  and  surface  reactions  appear  important 
in  most  EOIL  demonstrations. 

[02(o)]  emission  intensity  measurements  made  within  the  discharge  region  have 
been  reported.  Four  atomic  species  present  in  the  electric  discharges  emit  in  the  same 
spectral  region,  from  1,250  to  1,300  nm,  as  the  [02(o)]  state.  These  atomic  emissions 

are  from  high-energy,  13  eV,  or  ionic  O-atoms,  which  are  not  observed  downstream  of 
the  discharge  but  have  sufficient  concentrations  within  the  discharge  region.  By  using  a 
computer  model,  the  intensity  contributions  of  these  emissions  are  separated  from  the 
reported  [02(o)]  emission  intensities.  Using  the  computer  model  allows  [02(a)] 
intensity  measurements  made  directly  inside  the  discharge  region  to  reflect  the 
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concentration  of  the  [02(o)]  state  without  introducing  errors  that  result  from  changes  in 
the  atomic  oxygen.  Using  the  model  to  reproduce  the  [02  (<:/)]  emissions  permitted  the 
identification  of  both  the  0-0  and  1-1  vibrationally  excited  transition  of  [02(o)] .  The 

intensity  of  the  1-1  transition  is  sufficient  to  determine  a  temperature  that  agrees  well 
with  the  rotational  temperatures  reported.  Identification  of  these  emissions  within  the 
discharge  has  not  been  previously  reported. 

The  interpretation  of  the  measured  [02(o)]  yield,  using  a  streamlined,  nearly 
analytic  model,  cast  new  light  on  the  determination  of  the  kinetics  within  the  electric 
discharge.  Over  the  range  of  pressures  studied,  the  values  of  Ta  are  invariant  as  a 

function  of  the  flow  velocity  of  800  to  15,000  cm/s.  The  quenching  rate  at  p-wave 
frequency  field  oscillations  ranges  from  Tfl  =  4,000  to  9,000  s'1,  while  the  values  range 
from  V a  =  400  to  600  s'1  for  RF  fields.  The  only  rate  that  changes  significantly  between 

the  two  different  field  oscillation  frequencies  is  the  rate  of  ionization. 

From  the  steady-state  simulation,  it  has  been  detennined  that  the  product  of 

x  [02(a)]  is  approximately  constant.  As  the  velocity  of  the  gas  is  increased,  the 
temperature  of  the  gas  in  the  flow  decreases.  The  concentration  of  02(X)  is  inversely 

related  to  the  temperature  of  the  gas  and,  hence,  changes  dramatically  as  the  residence 
time  of  the  gas  is  changed.  This  requires  a  subsequent  change  in  the  electron  density. 

The  primary  kinetic  pathway  that  changes  as  the  velocity  changes  is  the  ionization  rate. 

As  the  concentration  of  the  ground  state  increases,  the  ionization  rate  drops;  therefore,  the 
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electron  density  drops,  allowing  the  product  of  the  two  species  to  remain  constant  over  a 
large  range  of  discharge  conditions. 
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V  Discussion 


Through  this  comprehensive  study  of  the  behavior  of  excited  species  within  an 
electric  discharge  of  flowing  oxygen  in  both  RF  and  p-wavc  frequency  E-fields,  a  few 
truths  may  be  stated.  First,  the  intensity  of  the  Oja )  emissions,  Ia  ,  obtains  a  steady- 
state  value  after  residing  within  the  discharge  region  for  approximately  0.3  ms.  This 
residence  time  corresponds  to  3  eV/  particle,  consistent  with  the  observations  of 
Rakhimova.17  The  modeled  results  of  Stafford  and  Kushner7  suggest  that  the 
Oja)  concentration  peaks  at  6  eV/  particle  and  is  limited  by  the  dissociation  of  the 
ground  state  oxygen.  When  O-atom  concentrations  make  up  more  than  half  the  flow  and 
the  ground  state  concentrations  diminish,  then  the  production  of  02(a)  by  Reaction  10  of 

Table  2.1  is  reduced.  However,  the  measured  results  of  Oja)  in  this  study  and 
Rakhimova’s  study  reach  a  steady  value  at  a  lower  energy  deposition  than  predictions, 
suggesting  there  is  a  mechanism  for  quenching  Oja)  that  is  not  currently  included  in 
simulations.  This  destruction  pathway  will  limit  power  scaling  in  EOIL  systems, 
because,  when  increasing  deposited  energy,  higher  yields  of  02(a)  will  not  be 
obtainable,  presumably  because  the  quenching  species  will  also  be  created  at  a  higher 
concentration. 

Secondly,  it  is  observed  in  a  large  number  of  different  studies  that  the  yield  of 
Oja  ) ,  F  ,  decreases  as  the  pressure  of  oxygen  is  increased  within  the  discharge. 

However,  high  yields  of  02(a)  have  been  achieved  at  high  total  gas  pressures  when  the 
partial  pressure  of  oxygen  was  small.  This  result  further  identifies  the  destruction 
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mechanism  as  being  tied  to  the  existence  of  oxygen  within  the  gas  flow.  Hence,  the 
destruction  mechanism  is  not  related  to  collisions  with  the  carrier  gas,  argon,  or  an 
excited  species  of  argon. 

A  steady-state  description  of  the  02(a )  concentration  has  been  used  to  further 
describe  the  behavior  of  the  destruction  tenn  of  02(a)  ,  which  is  seen  in  Equation  4.28. 
The  resulting  values  of  the  quenching  tenn,  Y  a ,  do  not  change  significantly  over  the 
pressure  range,  2  to  10  Torr,  and  residence  time,  0  to  3  ms.  Therefore,  whatever 
mechanism  results  in  the  destruction  of  the  02(a )  within  the  discharge  must  also  remain 
unaffected  by  changes  in  the  gas  pressure  and  flow  velocity  of  the  gas.  The  magnitude  of 
Ta  for  the  p-wave  case  is  6,000  s'1.  Because  Ta  must  be  the  product  of  a  rate  coefficient 

and  the  concentration  of  the  quenching  species,  the  magnitude  of  Tfl  restricts  reactions 
that  may  be  responsible  for  the  observed  behavior  of  02(a)  .  For  instance,  a  species  that 
makes  up  10%  of  the  gas  flow  at  2  Torr  and  a  steady  state  gas  temperature  of  1200  K  has 

16  3 

a  concentration  of  approximately  1x10  cm'  .  In  order  for  the  product  of  the 
destructive  species  and  the  reaction  rate  coefficient  to  equal  6,000  s'1,  the  rate  coefficient 

13  3  1 

of  this  reaction  must  be  6  x  10'  cm  s'  .  Furthermore,  the  Ta  values  derived  from  the 

RF  measurements  are  500  s'1.  In  summation,  the  quenching  mechanism  must  be 
invariant  to  changes  in  pressure  and  flow  velocity,  while  changing  by  an  order  of 
magnitude  when  the  frequency  of  the  E-field  is  changed  from  2.45  GHz  to  13.56  MHz. 

The  criteria  for  the  quenching  mechanism  of  02(a)  is  very  specific,  many  species 

commonly  thought  to  limit  the  production  of  02(a )  within  the  discharge  do  not  meet 
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these  criteria.  There  is  one  relationship  that  does  meet  the  criteria,  ye  Jx[Ch(X)]  is 

invariant  with  respect  to  gas  pressure  and  flow  velocity.  Also,  the  ionization  rate  changes 
by  an  order  of  magnitude  between  RF  and  p-wave  frequencies.  Furthermore,  the 
reactions  responsible  for  reducing  the  electron  density,  namely,  recombination,  diffusion, 
and  dissociative  attachment,  do  not  change  significantly  when  the  frequency  of  the  E- 
field  changes.  Note  that  the  dominant  loss  mechanism  of  electrons  is  dissociative 
attachment.  Dissociative  attachment  has  a  rate  that  is  an  order  of  magnitude  larger  than 
diffusive  losses  under  the  conditions  studied  in  this  work.  Under  these  conditions  the 
electron  densities  will  not  be  affected  by  interactions  at  the  wall.  Because  the  electron 
production  rate  decreases  by  an  order  of  magnitude  and  the  loss  rates  do  not,  it  is 
reasonable  to  expect  the  electron  density  to  decrease  by  an  order  of  magnitude  when  the 
frequency  of  the  E-field  changes.  So,  if  the  concentration  of  the  excited  species,  [M], 

responsible  for  quenching  02(a) ,  is  dependent  on  the  product  of  [e “]  x  [02(X)] ,  then  it 

should  be  invariant  with  pressure  and  gas  flow  velocity,  making  it  a  prime  suspect  for 
quenching  02(a)  ,  providing  the  reaction  rate  and  concentration  are  large  enough. 

5.1  Evaluation  of  02(a  1  A)  Destruction  Mechanisms 

In  Table  5.1,  six  candidate  reactions  are  listed  that  may  be  the  dominant 
destruction  pathway  of  the  02(a)  state.  Reactions  22-25  are  mentioned  by  Stafford  as 

possible  dominant  quenching  reactions.  Reaction  26  is  deactivation  by  direct  electron 
impact,  which  is  included  for  completeness.  The  final  reaction  (reaction  27)  is  a  new 
candidate  for  quenching  of  02(a)  introduced  in  this  work. 
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Table  5.1  Possible  quenching  reactions  responsible  for  limiting  Ya  . 


Reaction 

Number 

Reaction 

Reaction  Rate 
Coefficient 
cm3  s'1 

Ref 

22 

[02  (a)]  +  [02  (X)]  — »  2[02  (X)] 

_200/ 

k02x  ~3.0xl  0  18 1?  /r* 

29 

23 

[0(3P)]  +  [02(a)]-^0  +  [02(X)] 

Ko  =  2.0  xlO16 

30 

24 

[0(  3P)\  +  [02  (X)]  +  \02  (a)]  [0(  3P)]  + : 

klbody  =  1-0  X  1CT32 

cm6  s'1 

17 

25 

[03]  +  [02(a)]^\03]  +  [02(X)] 

-2840/ 

Km  =  2.0xl0-“e  /,s 

29 

26 

[e“]  +  [02(a)]-»[e"]  +  ^// 

*„*/(*■»)><  10"10 

Scaled 

27 

[<92(v)]  +  \02(a)]  [02m  +  [02(X)] 

*4xlO-10 

This  Work 

Table  5.2  Concentrations  for  relevant  species  determined  by  streamlined  code. 


Maj  or 
Species 

p-wave  (cm'3) 

@  2  -  7.5  Torr,  and  20  (cm  s'1) 

RF  (cm'3) 

@  4  Torr  and  20  (cm  s'1) 

[02(X)] 

1.1  x  1016  —  5.8  x  1016 

4.6  xlO16 

[02(a)] 

1.6  x  1015  -  7.5  x  1014 

5.4  xlO15 

[0(3P)] 

1.3  x  1016  -  3.1  x  1015 

1.8  x  1015 

[o3] 

1  x  1012 

lxlO12 

[el 

2.4  x  1012  -2.1x10" 

1.2x10" 

[02(v)] 

2  x  1014 

4.6  xlO13 

Table  5.3  Destruction  rates  for  relevant  reactions  determined  by  streamlined  code. 


Reaction 

Number 

Collision 

Partner 

Expected  Ta  s'1 
RF 

Expected  Ta  s' 

p-wavc 

22 

02(X) 

0.1 

0.1 

23 

O-atoms 

0.16-0.  54 

1.6  —  2.6 

24 

0(3P)  +  02(X) 

1.6 -5.4 

16-26 

25 

03 

~3 

~30 

26 

e' 

0.7-2 

4-25 

27 

02(v) 

100  -  700 

1000  -  9000 
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The  experimental  data  argue  for  a  quenching  tenn  that  has  a  value  of  ~  6  x  103  s'1 
for  p-wave  conditions  and  a  value  of  ~3  x  102  s'1  for  RF  conditions.  The  difference 
between  the  RF  and  p-wavc  cases  involves  different  gas  temperatures,  Tg;  electron 
temperatures,  Te;  and  the  expected  concentration  of  the  quenching  species.  The  expected 
values  of  quenching  rates  for  the  reactions  in  Table  5.1  are  reported  for  both  p-wavc  and 
RF  conditions  in  Table  5.2. 

Table  5.1  includes  approximate  values  for  reactions  22-27,  which  all  destroy  the 
[02(o)]  state.  The  values  shown  in  Table  5.2  are  approximated  for  a  typical  4  Torr 

discharge,  where  the  gas  temperature  is  1,000  K,  flow  velocity  is  20  cm/s,  and  the  flow 
tube  diameter  is  1  cm.  From  Table  5.1,  reaction  rates  22-26  are  too  slow  by  at  least  an 
order  of  magnitude  to  explain  the  observed  yields  of  02(a) .  A  more  in-depth  discussion 
of  each  reaction  follows. 

5.1.1  02(X)  Collisionally  Quenching  02(a, 1  A) 

The  collision  partner  with  the  highest  concentration  is  ground  state  oxygen.  If 
collision  relaxation  of  02(a)  with  02(X)  (reaction  22)  is  the  dominant  quenching 
mechanism  of  the  02(a)  concentration,  then  the  decay  rate  would  have  the  form: 

r ,=k0A°^X)\  (5-1) 

The  reaction  rate  k()  x ,  shown  in  Table  4. 1  and  quoted  from  Atkinson,”  has  a  value  of 

3  x  1CT18  cm3  s'1  or  less  for  all  gas  temperatures.  As  shown  in  Table  4.2,  combining  this 
with  the  average  value  for  [02(X)] ,  approximately  10  cm",  it  is  calculated  that  the 
value  of  Y  a  for  this  reaction  would  be  approximately  0.1  s'1.  This  value  is  much  too 
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small  to  explain  the  large  amount  of  quenching  observed  in  this  study.  Furthermore,  the 
amount  of  02(X)  in  the  gas  flow  does  not  change  significantly  from  the  RF  case  to  the 

p-wavc  case.  Therefore,  this  reaction  does  not  explain  the  difference  in  the  decay  rate, 

,  between  the  RF  and  p-wavc  discharges.  Ground  state  molecular  oxygen  does  not 
strongly  quench  the  02(a)  state  and  does  not  exhibit  the  needed  behavior  as  a  function  of 
discharge  frequency.  Because  there  is  not  a  significant  change  in  the  destruction  rate 
between  the  p-wavc  and  RF  case  using  this  mechanism,  this  reaction  cannot  be  the 
dominant  reaction  limiting  the  concentration  of  02(a )  within  electric  discharges. 


5.1.2  O-atom  Collisionally  Quenching  02(a, 1  A) 

The  possibility  that  O-atoms  are  responsible  for  limiting  02(a)  is  feasible  for  a 

number  of  reasons.  O-atoms  are  not  created  in  a  conventional  COIL  SOG;  they  are 
resistant  to  recombination;  they  are  predicted  to  make  up  at  least  10%  of  the  gas  flow  in 
the  discharge;7, 17  and  they  look  to  have  the  correct  functional  dependence  on  the  electron 
and  ground  state  concentrations,  as  shown  in  Equation  4.31: 

[<9( 3P)]  =  k°+^e  (5.2) 

k\Vall 

The  rate  for  O-atom  collisionally  quenching  the  02(a )  state  (reaction  23)  is  reported  in 
Table  4.1  and  is: 


3  -1 


kaQ  =2.0x10  cm  s 


(5.3) 


In  order  for  this  reaction  to  obtain  a  value  of  Ta  of  approximately  1,000  s'1,  the 

20  3  1 

concentration  of  O-atoms  must  be  on  the  order  of  10“  cm'  .  Because  there  are  only  10 
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cm"  oxygen  molecules  in  the  gas  flow,  this  explanation  is  impossible  at  the  pressures 
studied  in  this  work.  Although  O-atoms  meet  many  of  the  criteria  for  being  the  dominant 
quenching  mechanism  within  the  electric  discharges,  the  rate  of  the  reaction  is  not  fast 
enough  to  describe  the  data. 

5.1.3  O3  Collisionally  Destroying  C>2(a  ' A) 

Ozone  is  a  very  reactive  species  that  is  generated  within  the  discharge  and  is  not 
present  in  conventional  COIL  systems.  Because  ozone  reacts  so  quickly,  small  amounts 
of  it  within  the  gas  flow  may  have  significant  effects  on  the  makeup  of  the  gas.  It  is  also 
this  reactivity  that  limits  ozone  to  making  up  approximately  0.001%  of  the  gas  flow.  The 
mechanisms  likely  to  dominate  the  production  of  ozone  within  the  discharge  are  shown 
below: 


[0\  +  [02(X)]  +  [M\  k°i+  >[0,]  +  [M]  k03+  =2.0x  10  35  cm6  s"1 

[0  }  +  [02{a)\  k°-  >[03]  +  [e~]  =3-0x10  10  cm3  s'1  (5.4) 

[03\  +  [02(a)\  -  >2[02(X)]  +  0  ^03fl  =3.0x10  12  cmV 

Ozone  is  offered  as  the  mechanism  responsible  for  the  limitation  of  [02(a)]  within  the 

electric  discharge  studied  by  Stafford.7  According  to  Stafford: 

[O3]«1012  cm"3 

(5.5) 

\03]kma  «  1012  (3.0  x  10^12)  =  3.0C1 

7  12  3 

Using  the  ozone  concentrations  reported  by  Stafford  of  10  “  cm"  and  then  using  the  rate 
of  ko3a  provided  by  Stafford,7  the  value  of  Tfl  would  be  approximately  3  s'1.  By  using 
the  Stafford  model,  it  has  been  predicted  that  up  to  30%  yields  are  achievable  using  an 
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RF  discharge.7  However,  the  measured  yields  reported  by  Rakimova,17  as  well  as  this 
work,  are  less  than  10%  yields.  It  would  appear  that  the  Stafford  model  does  not 
compare  well  to  measured  yields.  Therefore,  it  is  likely  that  another  mechanism  besides 
quenching  by  ozone  is  responsible  for  limiting  Ya  in  the  experimental  apparatus. 

The  production  pathway  of  ozone  is  rich.  So,  making  a  definitive  statement 
concerning  the  pressure  dependence  of  ozone  is  very  complicated.  Therefore, 
eliminating  ozone  as  a  possible  candidate  for  the  quenching  mechanism  based  on 
pressure  dependence,  changes  with  flow  velocity,  or  dependence  on  gas  temperature  is 
not  possible. 

The  production  tenn  of  ozone  is  dependent  on  O-atom  concentration,  so  it  is 
expected  that,  as  electron  density  increases  and,  subsequently,  the  O-atom  population 
increases,  ozone  population  should  increase.  Therefore,  an  ozone-dependent  decay  term, 
Ta ,  is  expected  to  scale  properly  as  the  discharge  is  changed  from  RF  to  p-wave.  It  is 

unlikely  that  a  large  enough  quantity  of  ozone  is  realized  within  the  gas  flow  in  order  for 
it  to  dominate  the  quenching  of  02(a)  within  the  discharge. 

5.1.4  ()(' D)  as  Destruction  Mechanism  02(a  1  A) 

The  0(  ]D)  state,  like  ozone,  is  a  very  reactive  species  that  does  not  occur  in  a 
conventional  COIL  system  but  is  produced  in  significant  quantities  within  an  electric 
discharge  sustained  in  oxygen  gas.  Also,  like  ozone,  the  fast  reaction  rates  of  this  excited 
state  of  atomic  oxygen  with  both  the  wall  and  other  species  within  the  discharge  keep  its 
concentration  small.  The  role  played  by  0{ ' D)  is  contentious  within  the  literature. 
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Despite  the  uncertainty  in  the  kinetics  mechanism,  the  following  discussion  assesses  the 
potential  for  0{  lD )  as  a  dominant  02(a )  deactivation  partner. 

In  the  literature,  Kushner’s  research  groups  refer  to  0{  lD)  as  a  source  of 
02(a)  .79  Kushner’s  model  does  not  include  a  reaction  between  0:(a)  and  0(  lD )  , 
which  quenches  or  pools  the  02{a )  state.  However,  the  research  group  that  includes 
Ionin74  is  currently  using  the  reaction  rate  given  below: 

[02(fl)]  +  [0(IZ>)]->[02(jr)]  +  [0(3P)],  k_aW  =lxl011cm-V1  (5.6) 

80 

This  reaction  rate  is  quoted  from  the  paper  by  Doroshenko  et  al.  However,  this  reaction 
rate  is  not  used  by  Stafford,  Rakimova,  Hicks,  or  Hill.  Hence,  for  the  remainder  of 
this  analysis,  it  will  be  assumed  that  this  reaction  may  be  neglected. 

The  approximate  concentration  of  0( ' D)  is  calculated  assuming  that  the  source 
of  0( '  D)  is  direct  electron  impact  with  ground  state  oxygen  (reaction  14)  and  is 
destroyed  collisionally  by  O-atoms  (reaction  30): 


[0(3/5)]  +  [0(1D)]^[0(3/,)]  +  [6>(3R)],  kow  =8.0x1  O^cnrV  (5.7) 

Therefore,  the  change  to  the  concentration  of  0(  P)  as  a  function  of  time  is 
described  by  the  differential  equation: 


d[0{'D) 


dt 


=  t„.[e][02(X)]  -  kou  [0(  5P)][0(  'D)l  -  k'2ul02(X)][0('D) 


(5.8) 


Assuming  steady-state  and  values  for  the  rates  and  electron  concentration,  from  a  4  Torr 
gas  flow  at  20  cm/s,  the  concentration  of  0( ' D)  is  approximately: 
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lo  ('£>)]  =^M^W]=jHiM  =  (7xlr")(23xl0'2)^i,2xiO'W. 

L  J“  KuAO+X)]  k,ml  (2.6x10-") 


The  calculation  of  the  source  terms  for  *A  and  (reactions  32  and  3 1)  are  shown 

below: 

+  +  ^1Z)  =1-6x10  12e'67/?i cm  V  (5.10) 


[ 02(X )]  +  [0(‘Z>)]  [ 02(b )]  +  [<9(3P)] ,  kbW  =  2.56  x  10'V67/7;cm  V'(5.11) 

using  this  steady-state  value  for  [0( ' D)  ].  The  resultant  Ta  is  calculated  with  the 
reaction  for  0{lD)  as  a  source  of  02(a )  subtracted  from  the  quenching  reaction. 

So,  [0(*D)  ]ss  =  1.2  x  1019,  k  =  (1  X  10'17  -  1.5  x  10'18  )  =  8.5  x  10'18,  which  implies  that 
Ta=  102. 

This  value  for  gamma  is  an  order  of  magnitude  that  is  too  small.  Furthennore,  the 
functional  form  of  the  0( '  D)  steady-state  concentration  is  also  linearly  dependent  on  the 


electron  density: 


\o('D)\ 


k\d+  [e][P2(X)\ 

Kd+  [f] 

^o2\d\Pl^)\  ^02  Id 


(5.12) 


Hence,  this  mechanism  for  gamma  will  linearly  depend  on  the  gas  pressure.  Because  the 
product  of  the  reaction  rate  and  steady-state  concentration  are  not  sufficient  enough  to 
obtain  a  Tfl  of  1,000,  and  because  this  reaction  would  yield  a  gamma  that  depends  on  the 

gas  pressure,  0(  lD )  is  not  the  primary  quenching  mechanism  for  02(a) .  It  should  be 


noted  that  the  wall  rate  for  relaxation  of  0(  D)  is  approximately  1  x  105  s'1.  If  the 
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quenching  of  0{ '  D)  is  dominated  by  the  wall  and  not  the  reaction  with  the  ground  state, 
then  the  concentration  of  this  reactant  would  be  independent  of  the  background  gas 
concentration,  thereby  meeting  this  observed  criteria  for  the  quenching  mechanism  of 
Ck (a)  .  However,  the  reaction  rates  in  Equations  5.10  and  5.11  are  well  known  and  when 
combined  with  the  a  nominal  value  of  the  ground  state  concentration  of  ~10  cm"  result 
in  a  destruction  rate  of  ~105  s'1.  A  nominal  wall  rate  is  on  the  order  of  103  s'1  which  is 
two  orders  of  magnitude  smaller  than  the  collisional  destruction  rate  of  ()(  D) .  Thus  the 
dominant  destruction  method  of  0( ' D)  is  by  collisions  with  ground  state  and  therefore 
will  be  dependent  upon  the  electron  concentration  as  shown  in  Equation  5.12.  Hence  the 
destruction  of  02(a ) by  collision  with  0( 1  ZD)  does  not  fit  the  measured  data. 


5.1.5  Removing  Power  Loading  Effects 

Because  of  the  different  plasma  geometries  and  production  methods  studied 
herein,  the  power  loading  of  the  plasmas  varies  between  the  different  measurements.  In 
order  to  remove  any  power  loading  effects  from  the  results  in  Figure  5.1,  is  divided 
by  the  power  loading  term,  which  is  the  left-hand  side  of  Equation  4.33.  At  sufficiently 
slow  flow  velocities,  Equation  4.33  may  be  written  as: 


^M,,[°2(X)L 


V 


(5.13) 


This  expression  for  power  loading  is  used  in  Figure  5.1  to  put  the  decay  rates  detennined 
from  Stafford,7  Rakimova,17  and  this  work  on  the  same  footing.  In  Figure  5. 1,  it  is 
observed  that  all  the  values  of  Tfl  appear  to  fall  on  a  single  line  once  power  loading  is 
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removed  from  their  values.  Figure  5.1  also  shows  the  expected  functional  forms  of 
reactions  24,  26,  and  27.  A  discussion  of  these  three  reactions  follows. 


Figure  5.1  Relation  of  ra  /  [02(X)]  versus  [02(X)].  o  2  Torr,  □  3  Torr,  •  4  Torr,  ■  5  Torr,  T  7  Torr,  V 
current  work  RF  4  Torr,  •  Stafford  RF7,  A  Rakimova  RF17.  Expected  ra  for  —  02(v)  Reaction  27,  electrons 
Reaction  26,  and  -  three-body  Reaction  24. 


5. 1.5.1  e  as  the  Dominate  Destruction  Mechanism  of  C>2(a  ' A) 

The  rates  of  electrons  reacting  with  the  02(a)  state,  such  as  reaction  26,  are  not 
studied  exhaustively  in  the  literature.  Therefore,  values  for  these  reactions  are  scaled 
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from  the  values  of  the  electron  interactions  with  the  ground  state.7  In  this  study,  the 
absolute  value  of  this  rate  will  not  be  important  to  the  conclusions.  The  absolute  value 
will  be  assumed  to  be  of  the  same  order  of  magnitude  as  other  reaction  rates  within  the 
plasma,  and  its  trend  as  a  function  of  pressure,  gas  flow  velocity,  and  E-field  frequency 
will  be  used  as  detennining  factors  for  the  conclusion  of  this  paper. 

As  seen  in  Figure  5.1,  the  electron-dependent,  power  loading  corrected,  pseudo- 
first  order  decay  term,  shown  in  gray,  decreases  as  the  ground  state  concentration  is 
increased.  This  does  not  match  the  measured  results,  which  are  not  dependent  on  the 
ground  state  concentration.  Therefore,  this  quenching  mechanism  may  be  eliminated 
from  consideration  based  on  the  trend  of  this  line,  as  well  as  having  too  small  of  an 
absolute  value,  as  reported  in  Table  5.2. 

5. 1.5. 2  Three-Body  Recombination 

17 

The  three-body  reaction  is  included  in  the  results  of  Rakimova  and  is  highly 
defended  as  the  current  solution  to  resolving  the  difference  between  current  models  and 
measurements.74  The  rate  of  this  reaction  is  not  well  understood  and  has  been  determined 
on  the  basis  of  matching  simulated  [02(a)]  to  experimentally  measured  [02(o)] .  The 

reported  value  for  the  three-body  reaction  is  on  the  order  of  10'  cm  s'.  Using  this 
value  for  the  three-body  reaction  at  a  gas  pressure  of  4  Torr,  and  assuming  that  O-atoms 
make  up  approximately  10%  of  the  gas  flow,  the  following  value  of  Tff  is  obtained: 

\0]  +  \02(a)\  +  \02(X)\  ***  >[0  ]  +  2 [02(X)] 

*3**  =3.0x1 KT32  cm6  s'1 
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(5.15) 


r3^=W0][02(x)]  s-1 
Y3Body  «(3.°xl0_32)  (l.OxlO16)  (l.0xl017)*10  s'1 

1  7 

The  work  reported  by  Rakimova  is  based  on  the  results  of  an  RF  discharge.  Therefore, 
the  expected  value  of  gamma  is  approximately  100  s’1.  The  rate  of  the  reaction,  shown  in 
Equation  5. 15,  is  dependent  on  the  O-atom  concentration.  As  shown  in 
Equation  4.31,  O-atoms  scale  directly  with  electron  density;  therefore,  this  reaction  will 
scale  by  an  order  of  magnitude  when  changing  from  p-wavc  to  RF  discharges.  The 
scaling  by  an  order  of  magnitude  from  p-wavc  to  RF  discharges  makes  the  three-body 
reaction  rate  a  prime  candidate  for  explaining  the  observed  behavior  of  Ya  .  Not  only  is 
the  reaction  rate  too  small;  this  mechanism  is  not  independent  of  the  ground  state 
concentration,  as  shown  in  Figure  5.1.  The  power  loading  corrected  value  of  for  the 
three-body  reaction  is  shown  as  a  dashed  black  line  in  Figure  5.1.  This  term  predicts 
correct  values  at  higher  pressures,  but  it  diverges  from  the  required  values  at  lower 
pressures.  Therefore,  another  mechanism  is  sought  after  to  describe  the  measured  data. 

5. 1.5. 3  Vibrationally  Excited  Oxygen 

Vibrationally  excited  molecular  oxygen  is  produced  within  electric  discharges  by 
direct  electron  impact,  as  seen  in  Table  2.1,  reactions  3-9: 

k 

[e~]  +  \02  (X)]  — ]  +  \02  ( v)]  (5.16) 

The  sum  of  reaction  rates  3-9,  as  detennined  for  a  4  Torr  microwave  discharge,  using  the 
BolSig+  code,  is: 

k  *4xlC r10c/M3j_1  (5.17) 

v 
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The  vibrationally  excited  states  of  oxygen  are  quenched  by  both  atomic  oxygen  and 
ground  state  oxygen  at  the  rate  of  1.4  x  10'14  cm3  s'1,  as  reported  by  Atkinson.29  The  rate 
of  destruction  of  vibrationally  excited  oxygen  via  these  methods  is  approximately: 


k  x 
v  - 


O AX)  1  =  (l.4  x  10  14 ) (l  x  1017)  =  (2.8  x  103); 


(5.18) 


The  reaction  of  vibrationally  excited  oxygen  at  the  wall  is  fast.  The  wall  rate  is 
calculated  using  the  same  fonn  as  discussed  in  Section  4.6.1,  Equation  4.31: 


kWall 


V,J  0 

2  r 


(5.19) 


For  vibrationally  excited  oxygen,  y0  =  0.2 ,  as  reported  by  Stafford.7  Therefore,  the 
expected  rate  of  02(v )  destruction  at  the  wall  is  24,000  s'1.  Destruction  at  the  wall  is 

larger  than  destruction  by  O-atoms  and  molecular  oxygen  by  an  order  of  magnitude. 

Using  the  production  rate  by  direct  electron  impact  and  destruction  at  the  wall,  the 
expected  steady-state  value  of  vibrationally  excited  molecular  oxygen  is: 


02(v) 


k  .  e 
Vl  —  L  As 


02(X)_ 


(5.20) 


Wall 


Using  average  values  for  a  microwave  discharge  at  4  Torr, 
(4xl<rlu)(3xl0“)(4xl0“) 


02(v) 


14 _ -3 


=  2  x  10  cm 


2.4  xlO4 


(5.21) 
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If  Tfl  is  determined  by  reaction  27,  then  V a  would  have  the  form: 


r,=-Uo2(v)] 


(5.22) 


Because  [0_2(v)].„  is  expected  to  have  the  value  of  2  x  1014  cm'3,  for  Ta  to  have  a  value  of 
6,000  s'1  for  a  microwave  discharge,  the  value  of  kva  must  be  3  x  10'11  cm3  s'1.  The 
required  value  of  kva  is  within  the  limits  of  the  gas  kinetic  rate.  However,  this  reaction 
rate  is  five  orders  of  magnitude  faster  than  quenching  by  atomic  oxygen  (reaction  46) 
and,  therefore,  not  likely.  But,  if  this  quenching  rate  was  believable,  then  the 
concentration  of  vibrationally  excited  oxygen  does  fit  the  other  required  conditions 
mentioned  in  this  study. 

If  it  is  assumed  that  all  excited  species  within  the  discharge  eventually  relax  back 
to  the  ground  state  and  that,  while  relaxing  to  the  ground  state,  these  species  populate  a 
vibrationally  excited  state,  then  there  could  be  significant  secondary  sources  of  Oi(v).  If 
it  could  be  presumed  that  the  [02(v)]ss  could  obtain  a  value  approaching  10%  of  the  gas 
flow,  then  the  reaction  rate  would  be  sufficiently  large  enough  to  explain  the  values  of 
r  a  at  5  x  10'  cm  s'  .  While  this  rate  is  still  faster  than  most  reactions  between  two 
neutrally  charged  species,  it  is  on  the  same  order  of  magnitude  as  ozone  reactions. 

Vibrationally  excited  oxygen  is  a  plausible  candidate  for  quenching  the  02(a)  in 
the  EOIL  system.  The  02(v)  state  is  well  populated  in  the  EOIL  SOG  and  does  not  exist 

in  the  COIL  SOG.  It  trends  correctly  with  the  change  in  the  SCCM  in  the  measured  data. 
Most  importantly,  because  of  the  dependence  on  the  electron  concentration,  as  shown  in 
Equation  4.42,  this  reaction  explains  the  large  difference  in  gamma  between  the  RF  and 
p-wave  E-fields. 
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5.2  Dominant  Destruction  Pathway 

The  concentration  of  singlet  delta  oxygen  in  microwave  discharges  of  pure 
oxygen  reaches  a  steady  value  within  a  1  ms  discharge  residence  time,  independent  of 
discharge  pressure.  Also,  the  corresponding  yields  decrease  inversely  with  pressure.  A 
steady-state  analysis  of  the  kinetics  is  adequate  to  explain  the  observations,  assuming 
only  electron  excitation  of  02(a) ,  the  inverse  super-elastic  process,  and  an  unspecified 
third  channel  for  02(a)  destruction.  The  pseudo-first  order  decay  rate  for  the 

unspecified  process  is  about  6,000  s'1  in  the  microwave  discharge,  independent  of 
pressure  and  flow  rate.  The  same  rate  in  the  RF  discharge  is  much  slower,  300  to  500  s’1. 
Modeled  electron  number  density  is  also  considerably  less  in  the  RF  discharge.  As  a 
whole,  the  results  suggest  a  02(a)  removal  process  that  is:  (1)  rapid  compared  to  the 

super-elastic  rate;  (2)  induced  by  a  species  generated  by  electron  kinetics;  and  (3)  a 
second  order  reaction. 

A  steady-state  analysis  has  been  presented  of  many  species  potentially  responsible 
for  the  destruction  of  02(a)  within  the  discharge.  It  has  been  argued  that  the  collision 

partner  must  be  dependent  on  the  product  of  the  electron  and  oxygen  ground  state 
concentrations.  A  brief  look  at  the  steady-state  concentrations  of  likely  collision  partners 
is  shown  below: 


_  Ku\e] 

Knd[02{X)\  ko2U 
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|  0(  P)  = 


k0t[e-]  [02(X)] 


vWall 


°2m 


4v1-M,[°2W 


Wall 


All  of  the  steady-state  concentrations  are  directly  related  to  the  electron  density  and, 
therefore,  are  expected  to  exhibit  the  order  of  magnitude  change  as  the  E-field  frequency 
is  changed.  However,  many  collision  partners,  such  as  0( ] D) ,  are  not  directly 

dependent  upon  the  product,  J  x  [02(  A)] ,  and,  therefore,  would  not  exhibit  the  flat 

response  as  a  function  of  gas  pressure  and  flow  rate.  In  these  steady-state 
representations,  it  may  be  seen  that  both  the  0(}P )  and  02(v)  states  have  the  correct 
dependence. 

Even  though  O-atoms  have  the  correct  behavior  versus  E-field,  pressure,  and  flow 
velocity,  the  expected  value  of  Y  a  is  not  sufficient  enough  to  explain  the  behavior  of 

02(a) .  As  seen  in  Table  5.3,  the  expected  Ta  for  the  O-atoms  is  two  orders  of 
magnitude  too  small.  The  reaction  rate  between  O  and  02(a)  is  well  known;  so,  it  is 
highly  unlikely  that  quenching  of  07{a )  by  O-atoms  is  responsible  for  limiting  the 
population  within  the  discharge.  However,  other  excited  species  that  are  dependent  on 
O-atoms  but  have  much  faster  reaction  rates  with  02(a )  could  maintain  the  correct 
independence  of  pressure  and  flow  rate  and  provide  plausible  reaction  pathways  for  the 
destruction  of  02(a )  . 
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The  reaction  rates  governing  the  behavior  of  02  (v)  are  not  as  well  understood  as 


those  for  O-atoms.  The  source  terms  for  02(v)  within  the  discharge  are  not  well 
determined.  Direct  electron  impact  with  02(X)  creating  02(v)  is  included  within  the 
BolSig+  package;  however,  there  are  many  different  pathways  in  which  02(v)  could  be 
created  within  the  discharge.  O-atoms  recombining  at  the  wall  of  the  flow  tube  could 
return  to  the  gas  flow  in  the  form  of  02(v)  .  Any  excited  molecular  species,  02(b)  or 

02(A ) ,  that  is  collisionally  quenched  within  the  discharge  could  also  fonn  02(v)  while 

losing  its  energy  and  relaxing  back  to  the  ground  state.  O3  and  higher  energy  atomic 
oxygen  species  could  also  contribute.  If  these  pathways  contribute  significantly  to  the 
concentration  of  02  ( v )  within  the  discharge,  then  the  proposed  reaction  rate  would  not 
need  to  be  so  large,  making  this  destruction  mechanism  more  palatable. 
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VI  Conclusions 


The  directed  energy  community  has  been  pursuing  a  laser-based  weapon  system 
since  the  lasing  concept  was  first  demonstrated.  After  the  Airborne  Laser  Lab  (ALL) 
successfully  destroyed  multiple  airborne  targets  in  1983,  the  concept  of  an  ABL  weapon 
appeared  to  be  on  the  horizon.  The  current  ABL  system  employs  a  COIL  that  suffers 
from  an  extensive  logistics  trail.  Using  an  EOIL  eliminates  the  logistical  limitations. 
Existing  EOIL  systems  have  achieved  laser  powers  of  102  W,  but  scaling  to  kilowatts 
seems  to  be  limited  by  the  pressure  dependence  of  the  02(a)  yield.  Detennining  the 
chemical  kinetics  responsible  for  limiting  the  yield  is  paramount  to  finding  a  means  to 
increase  the  yield  as  a  function  of  pressure  or  detennine  a  different  means  to  scale  laser 
power. 

A  systematic  study  of  microwave  discharges  has  been  perfonned  through  the 
pressure  range  of  2  to  10  Torr  in  pure  oxygen  flows.  Gas  temperatures  of  300  to  1,200  K 
have  been  determined  to  a  higher  resolution  than  previously  reported  in  the  literature, 
subsequently  increasing  the  accuracy  of  the  reported  concentrations,  by  10%.  A 
theoretical  description  of  the  gas  temperature  allows  the  detennination  of  the  fraction  of 
discharge  energy  coupled  into  gas  heating  (15  to  30%).  The  simulation  of  02(a) 

molecular  emissions  and  O-atom  species,  which  emit  in  the  same  spectral  region,  allows 
concentration  measurements  to  be  made  directly  inside  the  discharge  region.  From  the 
measured  changes  in  02{a )  concentrations,  a  pseudo-first  order  decay  rate  has  been 
detennined.  The  change  in  the  decay  rate  as  a  function  of  gas  flow  velocity,  gas  pressure, 
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and  electric  field  frequency  has  detennined  that  a  second  order  tenn  directly  related  to 
electron  concentration  is  responsible  for  quenching  the  02(a )  state.  Many  reaction 

pathways  are  considered;  however,  no  clear  quenching  mechanism  sufficiently  describes 
the  observations. 

The  gas  temperature  measurements  reported  in  this  work  have  been  extracted 
from  the  rotational  energy  distribution  of  the  excited  02(b )  state.  The  resolution  of  the 
02(b)  spectral  data  reported  in  this  study  is  0.01  nm,  while  the  most  precise  resolution 

reported  in  literature  is  0.04  nm.  Higher  resolution  spectra  allow  the  determination  of  the 
gas  temperature  to  a  higher  degree  of  precision.  Typical  errors  in  reported  gas 
temperatures  for  this  study  are  5%  or  less,  whereas  errors  in  gas  temperature  reported  in 
the  literature  are  typically  10%  or  greater.  Because  the  gas  temperature  influences  the 
intensity  of  the  emissions  from  all  of  the  measured  species,  higher  precision  temperature 
detennination  results  in  higher  precision  density  measurements.  The  high-resolution  gas 
temperature  measurements  perfonned  in  this  study  improve  the  quality  of  the  reported 
data  throughout  this  work  by  approximately  10%. 

This  study  also  includes  a  theoretical  description  of  the  gas  heating  used  to  extract 
the  fraction  of  discharge  power  coupled  into  gas  heating.  Using  the  presented  method  to 
detennine  fractional  power  coupled  to  heating  was  applied  not  only  to  the  measurements 
made  by  this  group  but  also  those  reported  by  UIUC.5  The  importance  of  these  results  is 
related  to  the  efficiency  of  a  proposed  weapon  system  using  an  electrically  driven  SOG. 
The  17%  of  discharge  power  that  is  coupled  into  gas  heating  will  not  be  available  to 
pumping  energy  useful  to  the  laser  weapon  system.  From  the  calculations  in  Appendix 
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B,  in  a  typical  microwave  discharge  at  2  Torr  of  oxygen,  66%  of  the  power  is  removed 
from  the  discharge  by  atomic  oxygen,  and  only  1  %  of  the  power  is  removed  by  the 
02  ( a )  excited  state  of  oxygen.  This  effect  may  ultimately  limit  the  usefulness  of  an 
EOIL  system  for  weapons  applications. 

Measured  emissions  intensity  of  02(a)  are  also  affected  by  variations  in  the  gas 

temperature.  An  increase  in  the  gas  temperature  decreases  the  emitting  species 
concentration.  If  the  effect  of  gas  temperature  is  not  included  in  the  analysis  of  these 
emissions,  then  one  cannot  determine  whether  the  change  in  concentration  was  caused  by 
chemical  kinetics. 

02(a)  emission  intensity  measurements  made  within  the  discharge  region  have 
been  reported.  Four  atomic  species  present  in  the  electric  discharges  emit  in  the  same 
spectral  region,  from  1,250  to  1,300  nm,  as  the  02(a)  state.  These  atomic  emissions  are 

from  high-energy  or  ionic  O-atoms,  which  are  not  observed  downstream  of  the  discharge 
but  have  sufficient  concentrations  within  the  discharge  region.  By  using  a  computer 
model,  the  intensity  contributions  of  these  emissions  are  separated  from  the  reported 
02(a)  emission  intensities.  Using  the  computer  model  allows  02  ( a )  intensity 
measurements  made  directly  inside  the  discharge  region  to  reflect  the  concentration  of  the 
02(a)  state  without  introducing  errors  that  result  from  changes  in  the  atomic  oxygen. 
Other  studies  in  the  literature  report  yields  above  30%  but  laser  powers  of  p- watts.  This 
may  be  the  result  of  extracting  02(a )  concentrations  from  spectral  emissions,  including 
emissions  from  O-atoms.  Identification  of  these  emissions  within  the  discharge  has  not 
been  previously  reported. 
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The  interpretation  of  the  measured  02(a)  yield,  using  a  streamlined,  nearly 
analytic  model,  cast  new  light  on  the  detennination  of  the  kinetics  within  the  electric 
discharge.  The  quenching  rate  at  p-wavc  frequency  field  oscillations  ranges  from 
Ta  =  4,000  to  9,000  s'1,  while  the  values  range  from  ra  =  400  to  600  s'1  for  RF  fields. 

The  only  rate  that  changes  significantly  between  the  two  different  field  oscillation 
frequencies  is  the  rate  of  ionization.  This  result  leads  directly  to  the  detennination  that 
the  02(a)  quenching  mechanism  must  be  related  to  the  ionization  rate.  It  is  a  small  step 

from  this  result  to  the  assertion  that  electron  number  density  must  be  important  to  the 
quenching  term. 

Quenching  of  the  02(  a)  state  by  direct  impact  with  electrons  does  not  explain  the 
acquired  data.  As  shown  in  Table  4.2,  the  expected  rate  of  quenching  the  02{a )  state  by 

collisions  with  electrons  is  two  orders  of  magnitude  smaller  then  required  to  explain  the 
data.  Also,  the  quenching  rate  is  invariant  as  the  flow  velocity  changes,  whereas  electron 
density  changes  linearly  with  gas  flow  velocity.  Therefore,  direct  collisions  with  O^a) 

and  electrons  is  not  responsible  for  destroying  the  6k (a) .  The  product  of  electron 
number  density  and  ground  state  oxygen  concentration  is  constant  as  a  function  of  the  gas 
flow,  as  determined  by  the  power  balance  equation.  There  are  two  reactions  that  are 
dependent  on  the  product  of  electron  number  density  and  the  ground  state  concentration, 
which  will  limit  yield  in  the  correct  manner  to  explain  the  measured  results. 

The  two  reactions  proposed  to  explain  the  pressure  dependence  of  the  02(a) 
quenching  are  dependent  on  the  electron  density  in  different  ways.  The  three-body 
reaction  reported  by  Rakimova  is  appropriate  if  it  is  assumed  that  the  electron  density 
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increases  as  the  pressure  of  the  oxygen  is  increased.  If  the  electron  density  follows  the 
trend  in  average  electron  energy  and  decreases  with  pressure,  then  the  C>2(v)  quenching 
mechanism  proposed  in  this  work  is  the  remaining  explanation  to  the  limit  in  0-,  ( a ) 

yield.  This  new  quenching  mechanism  may  present  the  community  with  the  new  focus  in 
which  scaling  of  an  EOIL  system  to  powers  required  for  weapons  applications  is  realized. 

This  work  points  to  the  C>2(v)  as  a  possible  quenching  mechanism  of  O^a)  .  In 
order  to  determine  the  validity  of  this  mechanism,  the  concentration  of  C>2(v)  within  the 
gas  flow  should  be  verified.  Optical  detection  of  02(v)  is  complicated,  because  the 
energy  separation  of  the  excited  states  does  not  produce  photons  that  are  easily 
distinguishable  from  thermal  emissions.  The  primary  quenching  rate  of  02(v)  is  the  wall 
rate,  as  discussed  in  Section  4. 9. 4. 3.  Therefore,  the  term,  ra,  should  be  dependent  on 
changes  to  the  surface  area  of  the  discharge  region. 

In  conclusion,  a  large  step  was  made  towards  a  field-ready  system  in  switching  to 
an  oxygen-iodine  laser.  However,  in  choosing  the  chemistry  based  laser  system,  a 
logistics  issue  was  introduced.  Using  an  EOIL  on  the  airframe  may  be  the  solution  if 
scaling  to  kW  laser  powers  is  realized.  In  this  work,  it  is  argued  that  power  scaling  is  a 
complicated  endeavor  that  involves  complicated  chemical  kinetics  involving  a  second 
order  quenching  of  the  02(a)  state  used  as  the  energy  reservoir  for  oxygen-iodine  laser 
systems. 
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Appendix  A 

Section  A.  1  Determination  of  Laminar  Flow 

The  oxygen  gas  flow  is  assumed  to  be  a  laminar  flow  throughout  the  calculations 
in  the  document.  This  assumption  is  based  on  the  calculation  of  the  Reynolds  number  of 
the  experimental  apparatus.  Reynolds  numbers  of  less  then  500  are  generally  considered 
to  be  of  a  laminar  type.40  The  Reynolds  number  is  defined  as: 


Re  =  ( pvgD)/ju 


(A.l) 


where  p  is  the  density  of  the  gas;  vg  is  the  velocity  of  the  gas  flow;  D  is  the  diameter  of 
the  flow  tube;  and  //  is  the  dynamic  viscosity  of  the  fluid.  The  dynamic  viscosity  of  the 
fluid  is  a  function  of  the  gas  temperature,  Tg,  as  shown  below: 


M  =  A) 


T  +C  1 

1o  T  '-'0 

T  +C 

g  J 


(A.2) 


where  pn  is  the  reference  viscosity  at  the  reference  temperature  Tq\  Tg  is  the  temperature 
of  the  gas  flow;  and  Co  is  Sutherland’s  constant  for  the  gaseous  material  in  question. 

For  the  specific  case  of  oxygen  flowing  in  a  1  cm  diameter  flow  tube,  with  a  gas 
temperature  of  1,000  K  (as  used  in  this  study),  the  value  of  the  Reynolds  number  is: 

Re  =  0.0082  (A.3) 

where  the  density  of  oxygen,  p,  is  0.3899  kg/nr1,40  and  velocity  is  100  m/s.  The  value  of 
the  dynamic  viscosity  is  calculated  using  values  of:  reference  viscosity,  p0,  is 
20.18  Ns/  m2;40  reference  temperature,  To,  is  292.25  K;  Sutherland’s  constant,  Co,  is  127 
K40;  and  the  gas  temperature,  Tg,  is  1,000  K.  This  yields  a  value  for  the  dynamic 
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viscosity  of  47.5  Ns/  m  .  The  value  of  the  Reynolds  number  is  well  below  500,  where  a 
gas  flow  transitions  from  a  laminar.  Therefore,  it  is  a  good  assumption  that  this  flow  is 
laminar. 

Section  A.2  Thermodynamic  Evaluation  of  the  System 

The  derivation  of  the  energy  equation  shown  in  Section  2.5  will  be  discussed  in 
this  section.  The  expression  will  be  derived  from  the  first  and  second  laws  of 
thennodynamics.  The  form  of  the  solution  is  similar  to  the  energy  balance  equations 
found  in  Stafford  and  Kushner,7  Rakimova,17  and  Kays  and  Crawford.40 

Starting  with  the  first  law  of  thermodynamics,  assuming  that  work  perfonned  by 
the  system  is  dW,  defined  as  Pd  V: 

dQ  =  dE  +  PdV  (A.  4) 

where  Q  is  the  heat  added  to  the  system;  E  is  the  total  energy  of  the  system;  P  is  the 
pressure  in  the  system;  and  V  is  the  volume  of  the  system.  This  law  states  that  the  heat 
either  removed  from  or  deposited  into  the  system  must  go  to  mechanical  work,  expanding 
the  volume  of  a  container,  or  be  deposited  in  the  internal  energy  of  the  system.  The  work 
done  by  the  expansion  of  the  gas  is  obtained  by  assuming  that  the  change  in  the  volume 
may  be  calculated  using  the  equation: 

V  =  A  z  (A. 5) 

where  the  change  in  the  value  of  z  is  detennined  by  the  change  in  the  velocity  of  the  gas 
flow,  and,  therefore,  the  expression  is: 

z  =  t0vg  (A.  6) 
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2 

where  A  =  nr  is  the  area  of  the  flow  tube;  to  is  the  time  it  takes  the  gas  to  travel  an 
arbitrary  distance  (in  this  case,  the  length  of  the  discharge  assuming  constant  flow 
velocity  of  100  m/s);  and  vg  is  the  velocity  of  the  gas.  Combining  equations  A. 5  and  A. 6 
with  the  ideal  gas  law  and  mass  conservation  law,  the  change  in  energy  and  work  is 
calculated,  as  the  temperature  changes  from  300  to  1,200  K.  The  volume  change,  dV,  is 
2.4  x  10'  m  ,  and  the  pressure  is  equal  to  533  N/m"  at  4  Torr,  which  results  in  a  change 
of  1.2  x  10'  J  or  1.5%  of  the  total  power.  For  the  remainder  of  the  calculation,  this  small 
amount  of  power  loss  will  be  ignored. 

So,  the  first  law  of  thennodynamics  is  reduced  to  the  expression: 

dQ  =  dE.  (A.7) 


If  it  is  assumed  that  the  volume  is  filled  with  pure  oxygen,  then  the  total  energy  of  the  gas 
flow  takes  the  form  of  either  kinetic  energy  of  the  gas  or  into  the  internal  energy  of  the 
oxygen  molecule.  Therefore,  the  first  law  may  be  written  as: 


dQ  =  dE  =  N 


1  /  X2 

C.dT  +  —  m(dv) 

2 


(A. 8) 


If  the  gas  within  the  discharge  is  assumed  to  be  of  uniform  temperature,  then  a 
simple  convective  heat  transfer  coefficient  is  defined  as: 


ql  =  -hr(Tg  -  T()  )  (A. 9) 

where  ql  is  the  heat  flux  in  (J  s'1  m'2);  h  r  is  the  heat  transfer  coefficient  in  joules  per 

second  per  meter  squared  per  Kelvin;  To  is  the  temperature  of  the  wall;  and  Tg  is  the 
temperature  of  the  gas.  The  wall  of  the  gas  flow  tube  maintains  a  temperature  just  above 
room  temperature  of  approximately  400  K.  The  energy  lost  from  the  gas  to  the  flow  tube 
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wall  occurs  over  the  surface  area  of  the  discharge,  S.  This  is  also  known  as  Newton’s 
Law  of  Cooling.  The  incremental  change  in  the  internal  energy  of  the  gas,  dE,  is: 

dE  =  dU  =  fPm  -  ShT(Tg  -  T0)  (A.  10) 

where//5,,,  is  the  fractional  power  coupled  into  heat.  From  Equation  A. 7,  the  incremental 
change  in  the  energy  of  the  system  must  equal  a  change  in  the  total  heat  of  the  system, 
dQ,  and  the  following  relationship  may  be  written: 


N 


f  1 


PS 

v  2 


CpdTg  +  -m  (dv)  =  JP  -  Shr(Tg  -  T0) 


(A.  11) 


This  relationship  shows  that  a  fraction  of  the  power  coupled  into  the  gas  which  is 
thermally  conducted  to  the  walls  of  the  flow  tube.  The  remaining  energy  is  coupled  into 
the  excited  molecular  states  of  oxygen,  as  well  as  an  overall  change  in  the  flow  velocity 
of  the  gas. 

The  mass  of  molecular  oxygen  is  32  amu.  The  change  in  the  gas  velocity  due  to 
measured  changes  in  the  gas  temperature  is  determined  to  be  10  meters  per  second  or 
less.  Using  this  result  for  the  change  in  velocity,  the  kinetic  energy  associated  with 
velocity  of  the  gas  is: 

N-m(dv)2  «  0.03J  (A.  12) 

2 

Therefore,  the  contribution  of  change  in  kinetic  energy  to  the  internal  energy  of  the 
system  will  be  ignored,  and  Equation  A.l  1  becomes: 

NCpdTg  =  fPin  -  ShT(Tg  -T0)  (A.  13) 
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This  is  the  form  of  the  energy  conservation  equation,  which  is  used  in  Section  2.5. 
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Appendix  B 


The  theoretical  treatment  of  the  plasma  within  this  study  dictates  that  the  power 
coupled  from  the  electric  field  into  the  plasma  is  divided  among  elastic  and  inelastic 
collisions  of  the  electrons  with  ground  state  oxygen.  After  steady-state  is  achieved  within 
the  discharge,  taking  about  1  ms,  17%  of  the  power  has  been  coupled  into  heating  the  gas 
in  the  discharge;  4%  has  been  coupled  into  the  ionization  of  the  gas;  17.4%  has  been 
coupled  into  excited  species  of  oxygen;  and  the  remainder  has  been  lost  to  the  walls  and 
dissipated  by  forced  convective  cooling  into  the  surrounding  air.  Determination  of  these 
percentages  follows. 

The  percentage  of  power  coupled  into  excited  species  is  determined  by  the 
concentrations  calculated  from  the  measured  concentrations  of  O-atoms,  using 
actinometry,  and  02(a) ,  by  scaling  to  the  Rakimova  data,  and  in  combination  with  the 

simulated  results  presented  in  this  work. 

The  amount  of  power  convectively  removed  from  the  discharge  region,  PM,  is 
calculated  as: 

Pm=Av^[M]  (B.l) 

where  A  is  the  cross-sectional  area  of  the  flow  tube;  v  is  the  velocity  of  the  flowing  gas; 

A E  is  the  energy  used  to  form  the  excited  state;  and  \M]  is  the  concentration  of  the 
excited  species,  M,  at  steady-state,  in  the  gas  flow. 

Nominal  values  for  the  excited  species  are  taken  from  Table  5.2:  Oxygen  atoms 
5%  =  1.3  x  1016  cm  3 ,  02(a)  2%  =  1.6  x  1015cm  3 .  The  energy  required  to  form  the 
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excited  state,  02(a )  ,  is  1  eV.  Atomic  oxygen  is  primarily  produced  within  the  discharge 

via  two  pathways  (reactions  13  and  14): 

e  +02(X)  0(3P)  +  0(  3P)  +  e 

e  +  02  (X)  koD+  >  0(  3P)  +  0(  lD )  +  e 

Reaction  13  results  in  the  production  of  two  oxygen  atoms  in  the  ground  state  and 

requires  6  eV.  The  result  of  Reaction  14  is  one  O-atom  in  the  ground  state  and  one  O- 

atom  in  the  1 D  excited  state  which  requires  8.4  eV.  At  the  electron  temperature  in  which 

this  study  occurs  (approximately  3  eV),  the  reaction  rate  coefficients  of  these  two 

reaction  pathways  are  about  the  same.  Therefore,  the  energy  required  to  generate  a  pair 

of  oxygen  atoms  from  an  oxygen  molecule  in  the  ground  state  is  an  average  of  these  two 

energy  pathways  and  is  equal  to  7.2  eV.  For  energy  conservation  to  hold  when  two 

oxygen  atoms  recombine  at  the  wall  of  the  flow  tube,  fonning  a  ground  state  oxygen 

molecule,  the  energy  deposited  to  the  wall  will  be  7.2  eV.  Also,  the  amount  of  energy 

carried  out  of  the  discharge  region  by  atomic  oxygen  is  also  7.2  eV  per  pair. 

Using  Equation  B.l  and  the  transitional  energies  previously  discussed,  the  amount 

of  energy  removed  from  the  discharge  by  energetic  species  in  the  gas  flow  will  be 

approximated.  If  v  =  2,000  cm/s,  A  =  nr 2  =  ;r  (0.5)  =  0.78  cm 2 ,  and  Vol  =  2.4  cm 3 , 

then  the  convective  power  loss  by  atomic  oxygen  is: 

AvAE  [0\  =  [0.78  cm2 ]  [2, 000  cm  /  s]  [l  .3  x  1016  cm'3  ]  3.6  eV 
PQ=  12  W 

and  by  02(a)  is  approximately: 
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AvAE\jD2  (a)]  =  0.78  c/?z2^][2,000  cm  /  5]  [^1.6  x  101"  cm  3  X.tieV 


1  02a 


=  1  W 


(B.4) 


Energy  is  also  deposited  collisionally  at  the  walls  by  excited  species.  The  majority 
of  the  energy  deposited  at  the  wall  by  excited  species  is  carried  by  atomic  oxygen.  The 
thennal  transfer  due  to  O-atoms  recombining  at  the  surface  of  the  wall  may  be 
detennined  using  the  wall  rate  in  reaction  19  and  the  approximation  of  the  concentration 
of  O-atoms.  Therefore,  using  the  wall  rate,  the  power  transferred  to  the  flow  tube  walls 
by  this  reaction  may  be  calculated  as: 

kWall  [6>]Eo/A£  =  [1600  ^‘][l.3  x  1016  cm  3][2.4  cm3] 3.6  eV 

(B.5) 

Pow  =  43  W 


Where  Vol  is  the  volume  of  the  discharge  region.  The  sum  of  the  power  being  swept  out 
of  the  discharge  region  in  the  form  of  the  excited  species,  being  collisionally  deposited  by 
O-atoms  at  the  walls,  and  the  17%  via  heat  transfers,  accounts  for  70  W.  85  W  are  being 
coupled  into  the  discharge  from  the  electric  field.  Therefore,  all  of  the  power  is  not 
accounted  for.  However,  all  of  the  concentrations  within  this  study  are  measured  by 
methods  that  include  at  least  a  20%  error.  By  adding  20%  to  the  concentrations  of  02(a) 
and  O-atoms,  the  dissipated  power  becomes  73  W.  Furthermore,  there  are  many  excited 
species  not  accounted  for  in  this  power  analysis,  including  the  02(A)  state  created  by 
reaction  12  and  the  vibrationally  excited  species;  so,  it  is  not  surprising  that  the  totality  of 
the  deposited  power  is  not  accounted  for. 
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